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To  Tracy 


ABSTRACT 


Four  continuous  models  of  population  genetics  with  nonlinear 
birth  and  death  processes  are  considered.  The  first  model  is  concerned 
with  a  one-locus,  two-allele  problem  in  which  fertil ity  /  viabil ity  is 
incorporated.  Global  dynamics  of  this  three  population  system  is 
obtained.  The  second  model  is  a  predator-prey  interaction  model  in 
which  the  prey  population  is  subdivided  into  three  genotypes  and  its 
dynamics  is  given  by  the  first  model.  The  trade-offs  between  the  fertil¬ 
ity/viability  of  the  prey  genotypes  and  their  predation  functional 
responses  relating  to  questions  of  convergence  to  pure  strains  and  per¬ 
sistence  are  discussed.  The  third  and  fourth  models  are  the  analogues 
of  the  first  two  models  but  for  a  two-locus,  two-allele  problem  and  with 
no  fertility  /  viability  differences.  Global  convergence  to  the  Hardy- 
Weinberg  manifold  is  shown  for  the  third  model.  Some  results  relating 
to  the  question  of  convergence  to  pure  strains  as  well  as  to  the  question 
of  persistence  are  obtained  for  the  fourth  model. 
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CHAPTER  1 


INTRODUCTION 


Population  genetics  is  the  study  of  how  genetic  principles 
apply  to  a  population  and  is  principally  concerned  with  the  genetic  and 
phenotypic  properties  of  the  members  of  a  population.  Ecology,  on  the 
other  hand,  is  more  concerned  with  population  size  and  distribution  as 
well  as  population  interactions  such  as  competition,  predation  and 
mutualism.  Even  though  these  two  areas  of  population  biology  are  closely 
related,  most  mathematical  studies  tend  to  specialize  more  in  one  of 
these  areas  than  another.  The  purpose  of  this  dissertation  is  to  study 
a  number  of  mathematical  models  for  which  both  the  genetic  and  the 
ecological  components  are  present. 

The  usual  setting  for  problems  in  population  genetics  is  with 
difference  equations  rather  than  differential  equations.  Difference  equa¬ 
tions  in  population  ecology  are  appropriate  for  problems  of  non-overlap¬ 
ping  generations.  When  organisms  reproduce  continuously  the  differential 
equation  approach  seems  more  appropriate.  Since  environmental  factors 
act  continuously,  it  seems  more  appropriate  to  integrate  the  genetic  and 
ecological  components  in  a  continuous  formulation.  This  is  the  approach 
undertaken  in  this  thesis. 

Even  among  continuous  models,  it  is  convenient,  for  our  purpose, 
to  classify  them  into  three  types.  First,  there  are  the  pure  genetic  models, 
representatives  of  which  are  the  ones  studied  in  Hadeler-Liberman  (1975), 
But! er- Freedman-Wal tman  (1982),  and,  Hadeler-Glas  (1983).  In  these  models, 
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the  birth  and  death  rates  are  taken  to  be  constants.  In  some  cases  the 
death  rate  is  taken  to  be  zero,  and  it  is  more  natural  to  talk  about 
frequency  of  a  genotype  rather  than  the  number  of  organisms  of  that  geno¬ 
type,  since  the  number  of  organisms  grows  unboundedly.  It  should  be 
pointed  out  that  many  of  the  continuous  models  in  population  genetics  are 
of  this  type.  Among  them  two  famous  ones  are  the  Fisher-Wright-Haldane 
model  for  one  locus  problems  (see  Crow-Kimura  (1970),  Akin-Hofbauer 
(1982),  and,  Losert-Akin  (1983))  and  the  Kimura  model  for  two  loci  prob¬ 
lems  (see  Crow-Kimura  (1970),  Kimura  (1958),  and.  Akin  (1979,  1983)).  The 
second  class  of  models  is  the  ones  that  incorporate  nonlinear  birth  and 
death  processes.  Representatives  of  this  class  are:  Freedman-Wal tman 
(1978,  1982),  and.  Beck  (1982).  In  this  class  of  models  the  birth  and 
the  death  functions  are  singled  out,  and  the  coefficient  associated  with 
the  birth  function  determines  the  mating  structure  and  selection.  The 
third  class  of  models  is  the  ones  that  incorporate  environmental  effects, 
as  well  as  nonlinear  birth  and  death  processes.  Representatives  of  this 
class  are:  Freedman-Wal tman  (1978,  1982),  Beck-Keener-Ricciardi  (1982, 
1984),  and,  Beck  (1984). 

The  models  (3.1)  and  (5.1)  to  be  studied  in  this  thesis  belong 
to  the  second  class  and  models  (4.1)  and  (6.1)  belong  to  the  third  class. 
These  models  attempt  to  extend  the  model  studied  in  Freedman-Wal tman 
(1978,  1982)  to  include  fertility  /  viability  differences  and  to  two  loci 
problems.  It  turns  out  that  even  for  the  case  studied  in  Freedman-Wal tman 
(1978,  1982),  we  are  able  to  generalize  and  unify  their  results  as  one 
single  theorem  (Theorem  4.6). 

The  genetic  component  of  all  the  models  studied  in  the  papers 
mentioned  under  class  2  and  3  can  be  derived  from  the  models  of  Nagylaki- 
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Crow  (1974).  This  derivation  is  carried  out  in  Waltman  (1984).  The 
model  (3.1)  (one-locus,  two-allele  with  fertility  /  viability  differences) 
and  (5.1)  (two-locus,  two-allele  with  equal  fertility  /  viability)  studied 
in  this  thesis  are  extensions  of  these  models.  A  different  deri¬ 
vation  can  also  be  found  in  Appendix  1  and  2.  The  main  theme  for  models 

(3.1)  and  (5.1)  is  global  convergence  to  equilibrium.  The  recent  work 

of  Akin  mentioned  previously  shows  that  even  though  global  convergence  to 
equilibrium  is  to  be  expected  for  one  locus  problems,  periodic  solutions 
(through  Hopf  bifurcation)  are  possible  in  the  two  loci  problems  if  selec¬ 
tion  and  recombination  are  incorporated. 

In  model  (4.1)  (resp.  (6.1)),  an  ecological  component  in  the 
form  of  a  predator  is  added  to  model  (3.1)  (resp.  (5.1)).  Thus  in  model 

(4.1)  (resp.  (6.1))  the  prey  population  is  modelled  by  (3.1)  (resp. 

(5.1) ).  To  model  the  predator-prey  interaction,  the  standard  generalized 
model  of  Gause  (1934)  for  predator-prey  interactions  is  employed.  There 
is  a  vast  literature  on  the  generalized  Gause  model.  Of  particular 
importance  to  the  investigations  we  are  undertaking  in  this  thesis  are: 
the  local  stability  criteria  of  Rosenzweig-MacArthur  (1963)  and  Gause- 
Smaragdova-Wi tt  (1936),  the  phase-plane  analysis  of  Freedman  (1976),  and 
the  global  stability  criteria  of  Hsu  (1978)  and  Cheng-Hsu-Lin  (1981). 

For  model  (4.1),  there  are  two  types  of  selection  forces  at  work.  On 
the  one  hand,  there  is  the  selection  due  to  fertility  /  viability  differ¬ 
ences,  and,  on  the  other  hand,  there  is  the  selection  due  to  differential 
predation  functional  responses.  The  main  theme  here  is  to  determine  the 
interplay  of  these  two  selection  forces  in  relation  to  the  questions  of 
persistence  and  non-persistence.  It  should  be  pointed  out  that  the  recent 
work  of  Freedman-So-Wal tman  (1984)  shows  that  even  in  the  absence  of  fertil- 
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i ty  /  viability  differences,  it  is  possible  for  model  (4.1)  to  possess 
periodic  solutions.  These  periodic  solutions  are,  of  course,  predator 
mediated.  In  model  (6.1),  we  attempt  to  discover  some  similarities  and 
differences  between  the  one  locus  and  the  two  loci  theories  in  relation  to 
the  question  of  persistence  and  non-persistence. 

The  dissertation  is  organized  as  follows.  In 
Chapter  2,  a  number  of  mathematical  results  used  in  the  main  body  of  this 
thesis  will  be  described.  Proofs  are  provided  for  those  results  that 
appear  to  be  new.  In  the  next  four  chapters,  four  continuous  models  in 
population  genetics  are  discussed.  Chapter  3  is  concerned  with  a  one-locus, 
two-allele  model  with  nonlinear  birth  and  death  processes  as  well  as  fer¬ 
tility/viability  differences.  It  is  shown  that  the  dynamics  of  this 
model  is  trivial  (that  is  all  solutions  converge  to  some  equilibrium). 

In  Chapter  4,  the  model  studied  in  Chapter  3  is  extended  to  include  a 

predator.  We  obtain  conditions  under  which  only  one  of  the  prey  gamete 

types  survives  as  well  as  conditions  under  which  the  system  persists.  Of 
particular  interest  and  importance  is  the  interplay  between  the  fertility/ 
viability  and  the  predator  functional  responses  of  the  prey  genotypes. 

In  Chapter  5,  the  model  studied  in  Chapter  3  is  extended  to  the  two-locus, 

two-allele  case  but  with  no  fertility  /  viability  differences.  Again  the 
dynamics  of  this  model  is  shown  to  be  trivial  and  the  analogue  of  the 
Hardy-Wei nberg  equilibrium  relation  for  discrete  models  is  obtained.  In 
Chapter  6',  the  model  studied  in  Chapter  5  is  extended  to  include  a  preda¬ 
tor.  Same  questions  as  in  Chapter  4  are  raised.  It  is  shown  that  the 
predator  functional  response  of  the  double  heterozygote  can  change  a  non- 
persistent  system  (when  considered  as  one-locus  problems)  to  that  of  per¬ 
sistence.  Finally,  a  concluding  discussion  is  given  in  Chapter  7. 


CHAPTER  2 


MATHEMATICAL  PRELIMINARIES 


In  this  chapter,  a  number  of  concepts  and  results  from  the 
theory  of  autonomous  differential  equations  and  dynamical  systems  will  be 
discussed.  Some  of  the  following  sections  also  contain  a  number  of  simple 
but  useful  technical  results  which  are  new. 

2.1.  Standard  Results. 

The  portion  of  the  mathematical  theory  which  is  considered  as 
standard  or  well-known  will  be  listed  below  along  with  references  where 
they  can  be  found. 

(i)  Existence,  uniqueness  and  continuous  dependence  on  initial 
conditions  and  parameters  theorem:  Chapter  1  of  Coddington- 
Levinson  (1955). 

(ii)  Differential  inequalities:  Section  1.6  of  Hale  (1969)  and 
Section  1.4  of  Coppel  (1965). 

(iii)  Flows,  local  flows  and  semi-flows:  Chapter  2  of  LaSalle  (1976), 
Chapter  II  of  Sell  (1971),  and,  Chapter  1  of  Bhatia-Szego 
.(1967). 

(iv)  Deriving  a  flow  from  an  autonomous  differential  equation: 

Chapter  2  of  LaSalle  (1976)  and  Section  1.7  of  Hale  (1969). 
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(v)  Poincare  -  Bendixson  theory  and  the  theory  of  index  for  planar 
systems:  Chapter  16  of  Coddi ngton  -  Levi nson  (1955)  and  Chapter 
VII  of  Hartman  (1982). 

(vi)  Omega -limit  sets,  alpha- limit  sets  and  their  well-known  prop¬ 
erties:  Section  1.8  of  Hale  (1969)  and  Chapter  2  of  LaSalle 
(1976). 

(vii)  Liapunov  and  global  stability  for  sets:  Chapter  2  of  LaSalle 
(1976)  and  Chapters  1  and  2  of  Bhatia-Szego  (1967). 

b 

(viii)  Flow  equivalence  and  C  -equivalence:  Chapter  2  of  Irwin  (1980). 


2.2.  Conventions  and  Notations. 

In  this  section,  some  conventions  and  notations  that  are  used 
in  the  main  body  of  the  thesis  will  be  discussed. 

A.  Conventi ons. 

(i)  In  all  our  applications,  the  flows  are  defined  by  autonomous 
differential  equations  on  subsets  of  Rn  . 

(ii)  By  boundedness  of  an  orbit  (or  solution)  is  meant  boundedness 
for  all  positive  time. 
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B.  Notations. 


(i) 

R+  =  the  positive  cone  in  Rn 

(ii) 

cl(R^)  =  the  closure  of  R+  in  Rn 

=  the  non-negative  cone  in  Rn 

(Hi) 

b(R^)  =  the  boundary  of  R+  in  Rn 

(iv) 

For  1  £  i i  <  •  •  •  <  i^  <  n  ,  let  I  =  {i •  si (<}  and 

I'  =  { 1 » • • • »n)  \  I  . 

H  ...  Y  =  e  d(R+)  :  x.  >  0  for  i  in  I  , 

A.  ,  ,A-  T  | 

1  1  \ 

=  0  for  i  in  I'} 

(v) 

cl(H  .  )  =  closure  of  H  ...  in  Rn 

A-  ,  S  A  •  A-  >  s  A  . 

H  \  M  \ 

=  (x  e  cl(R^)  :  x^  =0  for  i  in  I'}  . 

For  x  in  Rn  and  1  <  i  <  n  ,  x^  or  (x)^  denotes  the 
i*^  component  of  x  . 

(vi ) 

d  denotes  the  Euclidean  metric  on  Rn  or  any  given  metric 

of  a  metri c  space. 

(vii) 

'Solutions  of  a  differential  equation  are  sometimes  denoted  by 

x(t)  with  x(0)  denoted  by  xq  . 
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(viii)  If  a  flow  (XS<J>)  comes  from  a  differential  equation 

x  =  f(x)  (x  £  Rn) 

and  there  is  a  need  to  emphasize  f  ,  then  $  is  also  denoted 
by  <J>f  . 

(ix)  Lip  denotes  the  class  of  locally  Lipschitz  continuous  functions 

(x)  w(x)  denotes  the  omega -limit  set  of  x  . 

(xi)  0(x)  (resp.  0+(x))  denotes  the  orbit  through  x  (resp.  the 

non-negative  semi-orbit  through  x). 


2.3.  Positive  Invariance. 

Consider  the  differential  equation 

X  =  f ( X )  (x  £  Rn) 

where  f  e  Lip.  If  this  system  is  to  model  interactions  of  populations 
we  should  expect  cl(R^)  to  be  positively  invariant.  A  necessary  and 
sufficient  condition  for  this  is  given  by  the  following  proposition. 


Proposition  2.1:  cl(R^)  is  positively  invariant  under  the  flow  defined 
by  f  if  and  only  if  for  all  1  <  i  <  n  ,  we  have, 


f . (xn  ,x.  ,  ,0,x. 

i  1  *i-l*i 


+1 


,x  ) 
n 


>  0 
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for  all  (x-j  ,•••  ,xi_1 ,0,xi+1  ,•••  ,xn)  j_n  c!(r")  . 

Proof:  First  suppose  there  exists  1  <  iQ  <  n  and 

x  =  (x,  »•  •  •  ,x.  _,,0,x.  ,,•••, x  )  in  cl(R^)  such  that  f.  (x  )  <  0  . 

1o”  no  n  1o 

★ 

Let  x(t)  be  the  solution  such  that  x(0)  =  x  .  Then  x.  (t)  is  strict- 

o 

ly  decreasing  for  t  small  enough.  Hence,  x.  (t)  <  0  for  t  >  0  small 

o 

enough,  contradicting  the  positive  invariance  of  cl(R^)  . 

Now  consider  the  sequence  of  differential  equations 

x  =  f ( x )  +  1  =  fN(x)  ,  for  N  =  1 ,2 ,  —  . 

N. 

Since  f.  (x1 ,  •  •  •  ,xi_1 ,0  ,xi  +1  ,•  •  •  ,xn)  >  0  for  all  (x-j , —  ,xi_1 ,0,xi+1 , 

•••,x  )  in  c1(r")  and  1  <_  i  £  n  ,  by  looking  at  the  flow  on  b(R^) 

and  using  an  argument  similar  to  that  above,  we  can  show  that  cl(R^)  is 

N 

positively  invariant  under  the  flow  4>  N  defined  by  f  for  all 

f 

N  =  1 ,2 , - • •  .  Next,  we  suppose  that  cl(R^)  is  not  positively  invariant 

under  the  flow  4>-  defined  by  f  .  Then  there  exists  an  1  <  iQ  <  n  , 

a  t  >0  and  an  x  =  (x^,*«*,x^  _-j,0,x^  +-|,***,xn)  in  cKr")  such 

0  0 

that 

d(<J)f  (x  ,t  )  ,  cl  (r")  )  =  d  >  0 

Now,  by  continuous  dependence  on  parameters,  (see  Theorem  7.4  on  p.  29  of 
Coddi ngton-Levi nson  (1955)), 

★  *  *  * 

4)  N(x  ,t  )  converges  to  4>f(x  ,t  )  as  N  -*  +  °°  . 
fN 
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But  this  leads  to  a  contradiction  because 

dU  N(x  ,t  ),cl(R"))  =  0  for  all  N  =  1,2, •••  , 
completing  the  proof  of  the  proposition. 

q.e.d. 

2.4.  Boundedness . 

Let  (X,<j>)  be  a  semi-flow  and  let  A  <=  X  . 

Defi nition:  (X,cj>)  is  said  to  be  A-dissipative  if 

(i)  A  is  compact  and  positively  invariant,  and, 

(ii)  for  all  x  e  X  ,  f(x,t)  -+  A  as  t  -»■  +®  . 

Remark:  ( X ,4>)  is  A-dissipative  implies  that  all  solutions  are  bounded. 


The  following  is  a  situation  which  is  applicable  to  all  the 
models  we  study  here  and  for  which  we  can  show  the  flow  is  A-dissipative. 

Consi der 


*i  =  fi ^X1 ^  =  1 »•** 
y  =  g(x-j  ,*••  »xn»y) 

where  f^  ,  g  €  Lip  and  c1(r"+1)  is  positively  invariant.  Let 
x  =  x-j  +  •••  +  xn  and  suppose  that 


(2.1) 


' 
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X  £  f  ( X ) 

kx  +  y  <_  kf(x)  -  sy 


(2.2) 


for  some  k,s  >  0  ,  where  f  £  Lip  satisfies:  f(0)  =  0  ,  f'(0)  >  0  , 
there  exists  a  unique  K  >  0  such  that  f(K)  =  0  ,  and,  f 1 ( K)  <  0  . 

Proposition  2.2:  Let  A  =  {(x-j  ,•••  ,xn»y)  £  cl(R^+"')  :  x  £  K  and 
kx  +  y  ^  K  (M  +  sK)}  where  M  =  max  (f(x)  :  x  £  [0,K]}  .  Then  system 
(2.1)  when  considered  as  a  flow  on  cl(R^+^  is  A  -  dissipative. 


Proof:  Consider  the  scalar  differential  equation  z  =  f(z)  .  By  the 
properties  of  f  ,  we  have: 

(i )  if  zo  £  ( 0 ,  K]  ,  then  z(t)^K  as  t  -►  +°°  ,  and, 

(ii)  if  zQ  £  [K,+°°)  ,  then  z(t)'»K  as  t  +00  . 

Using  a  standard  differential  inequality  argument  (see  Theorem  6.1  on 
p.  31  of  Hale  (1969))  and  the  fact  that  cl(p"+1)  is  positively  invari¬ 
ant,  we  can  show  that 

(iii)  if  xq  £  [0 , K]  ,  then  x(t)  £  [0 , K]  for  all  t  £  0  ,  and, 

(iv)  if  x  £  cl(pj)  ,  then  lim  x(t)  £  K  . 

0  t-*+°° 

From  (2.2) ,  we  have, 

(v)  kx  +  y  £  -s(kx+y)  +  k(f(x)  +  sx)  . 
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We  first  show  that  A  is  positively  invariant.  For  if  we  let 

(xlos*  *  *  ,xno,yo)  eA’  then  by  x(t)  £  K  for  all  t  _>  0  .  On  the 

other  hand,  by  (v),  we  have. 


kx( t)  +  y(t)  <  e"st(kxo  +  yo)  +  ke"bL(M  +  sK) 


-st, 


st  , 
e  di 


0 


=  j  (M  +  sK)  -  e'st(j  (M  +  sK)  -  (kxo+yQ)) 


Therefore,  if  kx  +  y  <  —  (M  +  sK)  ,  then 

o  o  —  s 


kx(t)  +  y(t)  (M  +  sK)  for  all  t  £  0 


We  now  show  that  for  all  (xn  ,***,x  ,y  )  e  cl(R^+^)  , 

lo  no  o  + 


(x1  (t)  ,•••  ,xn(t)  ,y ( t) )  +  A  as  t  +  . 


By  ( i v ) ,  it  remains  to  prove  that  for  all  e  >  0  ,  there  exists  T  >  0 
such  that 


kx(t)  +  y(t) 


£  j  (M  +  sK)  +  e 


for  all  t  >  T  . 


By  ( i v ) ,  there  exists  >  0  such 
Let  K'  =  max  {x(tl  :  t  e  [0 ,T-j  3 }  . 


that  x(t)  <  K  +  ^-  for  all  t  £  T-j  . 
By  (v),  for  all  t  £  T-j  we  have 


kx(t)  +  y(t)  £  e”st(kxo+yQ)  +  ke  st 


rT 


0 


^  [f (x(t) )  +  sx(t)]  eST  dx 


+ 


ke 


-st 


[f (x(x) )  + sx(t)]  eST  dx  . 
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1  .  - s ( t-T i ) 

Since  f(x)  £M  for  all  x  e  cl(R+)  ,  second  term  <  (K+-S  K' )  , 

and  third  term  <  j  (M  +  sK)  +  |  .  We  can  choose  T  >  T  so  that  for  all 
t  £  T  ,  the  first  and  second  terms  are  <  .  Then  for  all  t  >.  T  , 


kx(t)  +  y(t)  <_  j  (M  +  sK)  +  e 


as  desired. 


q.e.d. 


Remark: 
for  all 

cURf1 


The  above  proposition  also  shows  that  solutions  can  be  continued 
positive  time  and  hence  system  (2.1)  defines  a  semi-flow  on 

)  . 


2.5.  Persistence. 

Consider 


x  =  f(x)  (x  e  Rn)  (2.3) 

where  x.  (i  =  l,«‘*,n)  can  be  thought  of  as  the  number  /  density  /  biomass 
of  the  organisms  in  the  iLr  population  and  that  the  interaction  within  and 
between  the  ponulatinns  is  described  bv  system  (2.3).  Assume  that 
cI(rJJ)  is  positively  invariant  and  all  solutions  initiating  from  cl(R^) 
can  be  continued  for  all  positive  time. 
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Definition:  (i)  System  (2.3)  is  said  to  be  persistent  if  for  all  x 

in  r"  ,  1 im  x.(t)  >  0  for  all  i  =  !,•••, n  . 

t-H-00 

(ii)  System  (2.3)  is  said  to  be  non-persistent  if  it  is  not  persis¬ 
tent.  That  is,  there  exists  an  x  e  R*J  and  i  (1  <  i  <  n)  such  that 
1  im  x.  (t)  =  0  .  Equivalently,  the  omega-limit  set  of  x  has  non- 

f-H-oo  1  o  ”  0 

empty  intersection  with  the  boundary  b(R^)  . 

Remarks :  (i)  It  should  be  pointed  out  that  there  are  many  mathematical 

definitions  of  the  concept  of  persistence.  For  example,  in  McGehee  - 
Armstrong  (1977)  the  authors  defined  persistence  to  mean  the  existence  of 
an  attractor  in  r"  .  The  one  we  choose  to  use  here  is  the  one  used  in 
Freedman-Wal tman  (1984). 

(ii)  According  to  the  above  definition  of  non-persistence,  system 
(2.3)  will  be  non-persistent  if  there  exists  an  asymptotically  stable 
rest  point  on  b(R^)  . 


2.6.  Wazewski  Sets. 

In  this  section  the  basic  properties  of  a  Wazewski  set  will  be 
discussed.  More  information  can  be  found  on  pp.  24-25  of  Conley  (1978) 
and  on  Sections  X.2  and  X.3  of  Hartman  (1982). 

Definition:  Let  (X,<J>)  be  a  continuous  semi-flow  and  let  A  c  X  . 
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(i)  The  eventual  exit  set.  A0  ,  of  A  is  defined  as 

A0  =  {x  e  A  :  <f>(x,t)  i-  A  for  some  t  >  0} 

(ii)  The  immediate  exit  set,  A~  ,  of  A  is  defined  as 

A"  =  {x  e  A  :  <J>(x,[0,t))  i-  A  for  all  t  >  0}  . 

(iii )  A  is  called  a  Wazewski  set  if  it  satisfies: 

(a)  A  is  flow  closed,  i.e.,  if  x  e  A  and  <J>(x,[0,t))  c  A  , 
then  <J)(x, [0,tl )  c  A  ,  and , 

(b)  A  is  closed  relative  to  A0  . 

(i v )  The  time  of  exit  map  T  :  A0  ->  [0,+°°)  is  defined  by 
T(x)  =  sup  (t  _>  0  :  <J>(x,[0,tl  )  c  A)  , 
for  all  x  e  A0  . 


Theorem  2.3.  (Wazewski):  If_  A  is  a  Wazewski  set,  then  A”  is  a 
strong  deformation  retract  of  A0  and  A0  is  open  relative  to  A  . 

Proof:  The  statement  of  the  theorem  as  given  on  p.  24  of  Conley  (1977) 
is  for  flows  but  the  proof  remains  valid  without  modification  in  the  case 
of  semi -flows. 


q.e.d. 
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Another  useful  property  of  Wazweski  sets  is  the  following. 

Corollary  2.4:  If  A  is  a  Wazewski  set,  then  the  time  of  exit  map  is 
conti nuous . 

Proof:  This  is  contained  in  the  proof  of  the  theorem.  In  fact,  condi¬ 
tions  (a)  and  (b)  are  provided  to  ensure  the  time  of  exit  map  is  upper 
semi -conti nuous  and  lower  semi -conti nuous  respectively. 

q.e.d. 


2.7  Chain  Recurrence. 

In  this  section  the  notion  of  chain  recurrence,  first  used  by 
C.  Conley  and  R.  Bowen  will  be  introduced.  (See  pp.  36-38  of  Conley 
(1978)).  This  concept  turns  out  to  be  a  very  useful  kind  of  recurrence 
from  the  application  standpoint. 

Let  (X,4>)  be  a  flow  and  let  S  be  a  compact  (Hausdorff) 
invariant  subset  of  X  . 

Defi ni tion:  Let  U  be  an  open  cover  of  S  .  Let  x,y  e  S  and  t  >  0  . 

A  (U, t)  -  chai n  from  x  to  y  means  a  sequence 

7 

x  =  xi>-">Vi  =  *!tr— -‘n  (ti  - 

for  all  i  =  1  ,*••  ,n  and  such  that  for  each  pair  (<f>(x.  »t)  ,x..+-| )  (i  =  1  ,***,n) 
there  is  an  element  of  the  cover  U  containing  both  members  of  the  pair. 
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Definition:  The  chain  recurrent  set  R(S)  of  S  is  defined  to  be  the 
set  of  all  x  e  S  such  that  for  any  choice  of  cover  U  of  S  and  any 
t  >  0  there  exists  a  (U,t)  -  chain  from  x  to  x  . 

Definition:  S  is  called  chain  recurrent  if  R(S)  =  S  .  It  is  called 
strong  gradient  like  if  R(S)  is  totally  disconnected  (and  consequently 
is  equal  to  its  rest  point  set). 

Proposition  2.5:  Given  x  e  S  .  Then  w(x)  ,  the  omega-limit  set  of 
x  ,  is  chain  recurrent. 

Proof:  (See  p.  38  of  Conley  (1978).) 

q.e.d. 


2.8.  Invariant  Manifolds. 

In  this  section  the  invariant  manifold  theory  for  rest  points 
will  be  discussed.  Some  general  references  for  this  topic  are:  Chapter 
IX  of  Hartman  (1982)  and  Section  9.2  of  Chow-Hale  (1982). 

Definition:  Let  (X,<j>)  be  a  flow  on  a  smooth  manifold  X  with  metric 

n  * 

d  .  (Think  of  X  as  an  open  subset  of  R  .)  Given  a  rest  point  x 

of  <j>  ,  the  stable  set,  Ws(x  )  ,  of  x  is  defined  as  the  set 

{x  e  X  :  d(<J>(x,t),x  )  -*■  0  as  t  +°°}  . 

The  unstable  set,  Wu(x  )  ,  of  x  is  defined  analogously  by  reversing 


time  t  . 
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If  the  rest  point  x  is  hyperbolic  (i.e.  all  eigenvalues  of 
the  linearization  (Jacobian  matrix)  have  non-zero  real  parts),  we  have 
the  following  stable  manifold  theorem. 

Theorem  2.6  (global  stable  manifold  theorem):  If  the  flow  (X,<J>)  i_s 

k  *  s  * 

C  (k  _>  1)  ,  and  x  is  a  hyperbolic  rest  point,  then  W  (x  )  (resp. 

u  *  k  * 

W  ( x  ) )  are  C '  immersed  submanifolds  of  X  ,  tangent  at  x  to  the 

generalized  eigenspace  corresponding  to  those  eigenvalues  of  the  lineariza¬ 

tion  which  have  negative  real  parts  (resp.  positive  real  parts) . 

Proof:  (See  Theorem  6.17  on  p.  152  of  Irwin  (1980).) 

q.e.d. 

*  s  *  u  * 

Remark:  If  x  is  hyperbolic,  W  (x  )  (resp.  W  (x  )  )  is  also  known  as 

* 

the  stable  manifold  (resp.  unstable  manifold)  of  x  . 

★ 

When  some  of  the  eigenvalues  of  the  linearization  at  x  have 

zero  real  parts,  there  are  a  number  of  locally  invariant  manifolds  through 
* 

x  . 

Theorem  2.7  (strong  stable  manifold  theorem):  Consi der 

x  =  Ax  +  u(x)  (x£Rn) 

where  u  e  C1  ,  u(0)  =  0  and  Du(0)  =  0  .  If  the  n  by  n  matrix  A 
possesses  s  (s  >  0)  eigenvalues  having  negative  real  parts,  then  the 
set 
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WSS(0)  =  {x  e  X  :  lim  j  log  d(<J>(x ,t)  ,0) )  <  +  °°) 

t->  +O0 

is  a  locally  invariant  manifold  of  dimension  s  . 

Proof:  (See  Theorem  6.1  on  p.  242  of  Hartman  (1982).) 

q.e.d. 

Definition:  The  set  Wss(0)  is  called  the  strong  stable  manifold  of 
x  =  0  .  The  strong  unstable  manifold,  Wsu(0),  can  be  defined  in  a  simi¬ 
lar  manner. 

Theorem  2.8  (center  manifold  theorem):  Consider 

x  =  Ax  +  u(x,y)  (x  c  Rm) 

y  =  By  +  v(x,y)  (y  e  Rn) 

near  (x,y)  =  (0,0)  ,  where  all  the  eigenvalues  of  A  have  zero  real 

k 

parts  and  all  the  eigenvalues  of  B  have  non-zero  real  parts,  u,v  e  C 
(k  >  1)  ,  u(0 ,0)  =  0  ,  v (0 ,0)  =  0  ,  Du(0 ,0)  =  0  ,  and,  Dv(0,0)  =  0  . 

Then  there  exists  a  function  h  :  (Rm,0)  ->  (Rn,0)  such  that 

(i )  the  set 

^  Wc(0)  =  { (x ,y)  e  (Rm+n , (0 ,0) )  :  y  =  h(x) >  , 

is  a  locally  invariant  manifold  in  Rm+n  ,  and , 

(ii)  Wc(0)  consists  exactly  of  those  solutions  (x(t),y(t))  for 
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which  | y ( t )  |  (teR)  is  smal  1 . 

Proof:  (See  Theorem  2.1  on  p.  313  of  Chow-Hale  (1982.) 

q.e.d. 


Definition:  Wc(0)  is  called  the  center  manifold  of  (x,y)  =  (0,0)  . 

Theorem  2.9  (center  stable  manifold  theorem):  Consider 

x  =  Ax  +  u(x,y)  (x  e  Rm) 

y  =  By  +  v(x,y)  (y  e  Rn) 

near  (x,y)  =  (0,0)  ,  where  all  the  eigenvalues  of  A  have  non-positive 
real  parts,  all  the  eigenvalues  of  B  have  positive  real  parts, 
u,v  e  C^  (k  _>  1)  ,  u(0,0)  =  0  ,  v(0,0)  =  0  ,  Du(0,0)  =  0  ,  and , 

Dv(0,0)  =  0  .  Then  there  exists  a  C^  function  g  :  (Rm,0)  ->  (Rn,0) 
such  that 


(i )  the  set 


WCS(0)  =  { (x ,y)  £  (Rlimi ,  (0 ,0) )  :  y  =  g(x)}  , 


rni+n 


is  a  locally  invariant  C  manifold,  and. 


(ii)  Wcs(0)  consists  of  exactly  those  solutions  (x(t),y(t))  for 
whi ch  | y ( t ) |  ( t >0 )  is  smal 1 . 


Proof:  (See  Theorem  2.11  on  p.  319  of  Chow-Hale  (1982).) 


q.e.d 
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Defi ni tion:  Wcs(0)  is  called  the  center  stable  mani fold  of  (x,y)  =  0  . 
Similarly,  one  can  define  the  center  unstable  manifold,  Wcu(0),  of 
(x,y)  =  (0,0)  . 

Remark:  It  should  be  pointed  out  that  these  invariant  manifolds  are 
tangent  to  the  respective  generalized  eigenspaces  of  the  linearization. 


2.9.  Linearization  Theorems. 

Linearization  theorems  allow  us  to  say  something  about  the 
local  structure  of  a  flow  near  a  point.  To  start  with  there  is  the  so- 
called  "straightening-out"  theorem  which  states  that  near  a  regular  point 
(i.e.  not  a  rest  point)  of  a  C  flow,  C  change  of  coordinates  near 
the  point  can  be  made  so  that  locally  the  flow  becomes  a  parallel  flow 
(see  Theorem  21.6  on  p.  56  of  Abraham-Robbi n  (1967)).  In  the  case  of  a 
hyperbolic  rest  point  of  a  flow,  we  have  the  well-known  theorem  of 
Hartman  and  Grobman  which  says  that  locally  (i.e.  near  the  rest  point), 
there  is  a  flow  equivalence  between  the  original  system  and  the  linear 
system  that  arises  from  linearizing  the  original  system  at  the  rest  point 
(see  Theorem  7.1  on  p.  244  of  Hartman  (1982)  or  Theorem  7.2  on  p.  110  of 
Chow-Hale  (1982)).  In  the  general  case,  we  have  the  following  theorem 
due  to  Palis  and  Takens: 

Theorem  2.10  (Pal i s-Takens ) :  Consider 


x  =  Ax  +  u(x,y) 
y  =  By  +  v(x,y) 


(x  e  Rm) 

(y  e  Rn) 


(2.4) 
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near  (x,y)  =  (0,0)  ,  where  u,v  e  (k  _>  1)  ,  u(0,0)  =  0  ,  v(0,0)  =0  , 
Du(0,0)  =  0  ,  and ,  Dv(0,0)  =  0  .  Assume  that  all  eigenvalues  of  A  have 

zero  real  parts  and  all  eigenvalues  of  B  have  non-zero  real  parts.  Then 
system  (2.4)  is  locally  C°  -  equivalent  to  the  decoupled  system 


w  =  Aw  +  u(w,h(w))  (w  e  Rm) 

z  =  Bz  (z  £  Rn) 


(2.5) 


near  (w,z)  =  (0,0)  ,  where  h  defines  the  center  manifold  of  (0,0)  . 


Proof:  See  the  second  theorem  on  p.  341  of  Palis-Takens  (1977).  This 
fact  is  also  mentioned  in  Remark  2.16  on  pp.  322-323  of  Chow-Hale  (1982). 

q.e.d. 

Using  Hartman-Grobman ' s  theorem,  we  can  prove  the  following 
proposition  which  tells  us  when  the  stable  manifold  of  a  boundary  hyper¬ 
bolic  rest  point  will  have  empty  intersection  with  the  positive  cones  R+  . 

Proposition  2.11 :  Consider 

x  =  f ( x )  (x  e  Rn) 


and  assume  that  cl ( Rn )  is  positively  invariant.  Let  xQ  e  b(R")  be  a 

/ 

hyperbolic  rest  point.  If  Wu ( xQ )  n  (Rn  \  cl(R^))  t  (1  ,  then 
WS(xQ)  n  r"  =  0  . 


9 
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c  n 

Proof:  Suppose  not,  that  is  suppose  W  (x  )  n  R+  i  0  .  If  we  can  show 
that  for  fixed  x-j  e  W  (xQ)  ,  X£  e  W  (x  )  where  x-j  ,X£  t  xQ  ,  there  is  a 
point  x  which  is  arbitrarily  close  to  x-|  and  a  time  t  >  0  so  that 
cf)(x ,t )  is  arbitrarily  close  to  x^  ,  then  this  will  contradict  the 
positive  invariance  of  cl(R^)  .  Since  this  can  essentially  be  reduced 
to  a  local  problem,  it  suffices  to  establish  the  following. 

Consi der 


x  =  Ax 

(x 

6  Rn) 

y  =  By 

(y 

e  Rm) 

where  all  the  eigenvalues  of  A  have  negative  real  parts  and  all  the 
eigenvalues  of  B  have  positive  real  parts.  Fix 

★  *  m+n  *  * 

(x  ,0)  and  (0 ,y  )  e  R  with  x  f  0  ,  y  f  0  .  Then  there  exists 

(xN,yN)  e  Rm+n  and  tN  >  0  such  that  (xN,yN)  +  (x  ,0)  and 

*  * 
4»(  (xN  ,yN)  ,tN)  -  (0,y  )  as  N  +  +°°.  For  a  proof,  let  tN  =  N  ,  xN  =  x 

and,  =  e  y  .  Then  clearly  x^  •*  x  and  yN  ->  0  ,  since  all 
eigenvalues  of  B  have  positive  real  parts.  Furthermore  we  have, 


t.,A  tNB  1mA  * 

e  yN)  =  (e  x 


x 


N 


* 

y  ) 


-y 


(0  ,y*) 


as  N  ■>  +00  as  desi red. 


q.e.d. 
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2.10.  Saddle  Nodes. 

In  this  section  a  sufficient  condition  for  the  existence  of 
saddle  nodes  and  some  related  calculations  will  be  discussed.  (See  pp. 
323-329  of  Chow-Hale  (1982).) 


Consider  the  system 

y-l  =  p(y  ) 

y  =  Ay  +  Q(Y) 


(2.6) 


where  y  =  (y2»***»yn)  e  »  Y  =  (y-j  ,y)  e  Rn  ,  P  and  Q  are  of  order 

two  or  higher  in  Y  and  all  the  eigenvalues  of  the  (n-1)  by  (n-1) 
matrix  A  have  negative  real  parts.  Clearly  Y  =  0  is  a  rest  point  of 
system  (2.6)  and  the  variational  matrix  of  (2.6)  at  Y  =  0  has  0  as  a 
simple  eigenvalue  and  the  remaining  (n-1)  eigenvalues  have  negative 
real  parts. 


Theorem  2.12  (saddle-node):  Let  y  =  ^-|(y-|)  with  ^(0)  =  0  be  the 
solution  of 


Ay  +  Q(y-|  ,y)  =  0 
near  (y-|  ,y)  =  (0,0)  and  let 

7 

'P2(yi  )  =  p(yn  .^(yp) 


(2.7) 


(2.8) 


If 
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^2(y-|)  =  n y i  +  °( |y-|  I3) 

and  n  f  0  then  Y  =  0  is  a  saddle  node:  the  strong  stable  manifold  of 
Y  =  0  is_  (n-1)  -dimensional  and  is  tangent  to  the  hyperplane  y-j  =  0  , 
the  1  -  dimensional  center  manifold  of  Y  =  0  is  tangent  to  the  y-|  - 
axis,  and,  Y  =  0  is  stable  on  one  side  of  its  strong  stable  manifold 
and  unstable  on  the  other  side. 

Proof:  Refer  to  Theorem  65  on  pp.  340-346  of  Andronov-Leontovich-Gordon- 
Maier  (1975)  for  the  case  when  n  =  2  and  pp.  326-327  of  Chow-Hale  (1982) 
for  the  general  case. 

q.e.d. 


Remarks :  (i)  It  turns  out  that 

n=2^~7(0)  (2>9) 

sy-i 

(ii)  If  n  >  0  ,  then  the  domain  of  stability  of  Y  =  0  is  on 
the  side  of  the  strong  stable  manifold  for  which  y-j  <  0  ;  if  n  <  0  ,  it 
is  just  the  reverse. 

The  above  theorem  is  stated  in  "canonical"  form.  In  applica¬ 
tions,  it  is  often  necessary  to  transform  the  problem  into  this  form. 

Some  of  that  calculation  will  be  carried  out  in  the  following. 

Consider  the  system 


Z  =  G(Z) 


(2.10) 
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where  Z  e  Rn  ,  G(0)  =  0  ,  and  the  variational  matrix  of  system  (2.10) 
at  Z  =  0  has  0  as  a  simple  eigenvalue  and  the  remaining  (n-1)  eigen¬ 
values  have  negative  real  parts.  Let  v-j  be  an  eigenvector  correspondi ng 
to  the  eigenvalue  0  and  let  v2»***>vn  be  a  basis  of  the  generalized 
eigenspace  correspondi  ng  to  the  remaining  (n-1)  eigenvalues.  Let 

J-  L_ 

N  =  E v-j » •  •  •  »v  1  be  the  n  by  n  matrix  whose  l  n  column  is  v^  and 
let  Y  =  N  "'z  .  Under  this  transformation,  (2.10)  takes  the  form  (2.6). 
The  y-|  =  0  hyperplane  is  transformed  back  to  the  hyperplane,  in  the  Z- 
coordinates,  spanned  by  the  vectors  v2s***»vn  (denoted  by  <v2 >* * ’ »vn>  ) , 
the  y-j  -  axis  is  transformed  back  to  the  line  spanned  by  the  vector  v-| 
(denoted  by  <v-j>  ),  and  the  positive  y-|  axis  is  transformed  back  to 
the  half  line 


Ov,  €  Rn  :  X  >  0} 


Let  N  =  [n 


G(Z)  =  (G1(Z),---.Gn(Z))  ,  and 


G,(Z) 


l 


n 


n 


+  0(|z|3) 


(i  =  1 ,  •  -  -  ,n)  .  Also  let  1^  =  (0 , -  - 
column  and  the  rest  are  all  0 1  s 


0's  ,  and, 


0,2,0,***,0)  with  a  "2"  in  the 


,0)  with  1 1 s  in  the  hth  and  kth 


columns  and  the  rest  are  all  0's  .  Then 
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Hence,  if  n  4  0  ,  Z  =  0  will  be  a  saddle  node.  The  strong  stable 

r  q 

manifold  W  (0)  of  1=0  is  tangent  to  <V2»***,vn>  and  the  center 

Q 

manifold  W  (0)  is  tangent  to  <v-|>  . 


2.11.  Omega-Limit  Sets. 

In  this  section  a  number  of  simple  but  useful  technical  results 
relating  to  omega-limit  sets  will  be  presented.  These  results  will  be 
used  repeatedly  in  the  main  body  of  the  thesis. 

Proposition  2.13:  Consider 


X  =  f ( X )  (x  €  Rn)  (2.11) 

n  ^ 

where  f  e  Lip  .  Assume  that  cl(R+)  is  positively  invariant,  let  x 

be  a  point  on  b(R+)  wi th  x-  =  0  and  such  that  f. (x  )  >  0  .  Then 

*  n 

x  4  w(x)  for  any  x  e  cl(R+)  . 

*  2 
Proof:  Since  x  is  a  regular  point,  if  f  is  C  we  can  use  the 

"straightening-out  theorem"  to  construct  a  change  of  coordinates 
* 

near  x  so  that  locally  the  flow  becomes  a  parallel  flow.  From  this 

and  the  positive  invariance  of  cl(R^)  ,  we  can  easily  see  that 

x  4  w(x)y  for  any  x  e  cl(R^)  .  In  the  general  case  when  f  e  Lip  , 

★ 

f. (x  )  >  0  implies  that  x^(t)  is  strictly  increasing  near  t  =  0  , 

★ 

where  x(t)  denotes  the  solution  of  system  (2.11)  such  that  x(0)=x  . 
Since  x.,-  =  0  ,  for  sufficiently  small  t  <  0  ,  x(t)  4  cl(R+)  .  On  the 
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other  hand,  if  x  €  w(x)  for  some  x  €  cl(R^)  then  x(t)  e  w(x)  also. 
By  choosing  t  <  0  sufficiently  small,  this  implies  w(x)  £  cl(R^)  , 
contradicting  the  positive  invariance  of  cl(R+)  . 

q.e.d. 

Lemma  2.14  (Butler-McGehee) :  Consider 

x  =  f(x)  (x  e  Rn) 

1  *  n 

where  f  £  C  .  Let  x  £  R  be  a  hyperbolic  rest  point.  If 

*  n  * 

x  £  w(x)  for  some  x  e  R  ,  then  either  w(x)  =  {x  }  or  both 

w(x)  n  Ws(x  )  and  w(x)  n  Wu(x  )  contain  points  other  than  x  . 

Proof:  (Refer  to  Appendix  1  in  Freedman-Wal tman  (1984).) 

q.e.d. 

Remark:  This  lemma  together  with  Proposition  2.11  provide  us  with  a 
useful  tool  for  proving  persistence. 

The  following  proposition  shows  that  under  certain  conditions 
0  is  not  in  the  omega-limit  set  of  any  orbit.  The  point  here  is  that 
the  flow  is  only  continuous  and  we  may  not  be  able  to  use  the  Butler- 
McGehee  lemma  as  stated. 

Proposition  2.15:  Consider 


x.j  =  f  i(x1,**«,xn,y)  (i  =!,•••»  n) 
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y  =  g(x] ,••• ,xn,y) 

where  f..,g  e  Lip,  cl  )  is  positively  invariant,  and  0  is  a  rest 
point.  Let  x  =  x-j  +  •  •  •  +  xp  and  suppose  that  there  is  a  neighbourhood 
U  of  0  i_n  cl(R^)  such  that: 

x  £  k-jX  and  y  £  -k^y  (2.12) 

in  U  ,  for  some  k-j , k^  >  0  .  Moreover,  assume  that  the  non-negative 
y  axis  is  invariant  and  0  is  globally  stable  over  it.  Then 
0  i  w(x1  ,•••  ,x  ,y)  for  all  (x-|,***,x  ,y)  e  rJJ  for  which 
0+  ( x-j , •  •  •  ,x  ,y)  is  bounded. 

Proof:  Let  e  >  0  be  such  that  (2.12)  is  valid  on 

V  =  {(x-j  ,•••  ,xn,y)  e  cl(<+1 )  :  x  <  c  ,  y  <  e)  . 
k,t  -k„t 

Since  x(t)  >  x(0)  e  and  y(t)  <  y(0)  e  c  ,  the  eventual  exit  set 
V°  and  the  immediate  exit  set  V  are  given  by: 

V°  =  {(x-j  ,•••  ,xn»y)  e  V  :  x  >  0} 

V"  =  { ( X-j , •  •  •  ,xp ,y )  e  V  :  x  =  e}  . 

Therefore  ,  V  is  a  Wazewski  set  and  hence  the  time  of  exit  map  is  contin¬ 
uous.  Based  on  this,  we  can  complete  the  proof  using  an  argument  similar 

to  the  proof  of  the  Butler-McGehee  lemma. 


q.e.d. 
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The  following  is  an  extension  of  the  Butler-McGehee  lemma  to 
general  (not  necessarily  hyperbolic)  rest  points. 

Proposition  2.16:  Consider 


X  =  f(x)  (x  e  Rn) 

9 *  *  * 

where  f  e  CT  .  Let  x  be  a  rest  point.  ]f  x  s  w(x)  for  some 

x  e  Rnand_  x  Wcs(x*) ,  then  w(x)  n  Wcs(x  )  contains  a  point  other  than 

* 

x  . 


Proof:  Use  Palis-Takens  linearization  theorem  to  construct  a  closed 
neighbourhood  U  of  x  such  that  V  =  U  \  W  (x  )  is  a  Wazewski  set 
with  V°  =  V  .  We  can  then  repeat  the  argument  as  in  the  proof  of  the 
Butler-McGehee  lemma  to  complete  the  proof. 

q.e.d. 


CHAPTER  3 


A  ONE-LOCUS,  TWO-ALLELE  MODEL  WITH 

FERTILITY  /  VIABILITY  DIFFERENCES 


In  this  chapter,  the  global  convergence  to  rest  points  of  a 
continuous,  one-locus,  two-allele  genetic  model  that  incorporates  non¬ 
linear  birth  and  death  rates  as  well  as  unequal  genotypic  fertilities/ 
viabilities  will  be  analyzed.  This  model  will  be  derived  in  Appendix  1. 


3.1.  Introduction  and  the  Model. 


This  chapter  is  concerned  with  a  dynamical  system  of  the  form 


X  =  — j  B(x)  -  -T  D(x) 
(u+v) 


where 


-  ~x~  D(x) 


(u+v)' 
2 


X3 


v 


x. 


(u+v) 


?  B(x)  -  -f  D(x) 


(3.1) 


V  = 


flXl  +  2  f2x2 


f0X0  +  7T  fnx 
3  3  2  2  2 


(3.2) 


x  =  x-,  +  x  +  x 


(3.3) 
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fl  ’  f2  »  f3  >  0  (3.4) 

and 

=  dt”  =  c*en" vative  with  respect  to  time  t. 

Here  x-j  ,  ,  and  x^  denote  the  number  of  organisms  of  geno¬ 

type  AA  ,  Aa  ,  and  aa  respectively  of  a  one-locus,  two-allele  problem 
where  A  and  a  stand  for  the  two  allele  types.  System  (3.1)  is  a  gener¬ 
alization  of  the  model  (without  predator)  studied  in  Freedman-Wal tman 
(1978,  1982).  In  these  two  papers,  the  f. 's  are  taken  to  be  unity. 

In  order  to  preserve  some  consistency,  the  same  notations  as 
in  these  two  papers  will  be  adopted.  Models  similar  to  system  (3.1) 
have  been  used  by  K.  Beck  in  her  study  of  the  population  genetics  of 
cystic  fibrosis.  Beck  (1982).  (See  also  Beck-Keener-Ricciardi  (1982, 

1984)  for  related  models.)  In  Section  3.1  of  Waltman  (1984),  there  is  a 
nice  discussion  on  the  derivation  of  these  types  of  models.  The  effect 
of  the  introduction  of  a  predator  population  to  the  persistence  of  the 
system  for  which  the  prey  population  (in  the  absence  of  the  predator)  is 
modelled  by  system  (3.1)  will  be  studied  in  Chapter  4.  Other  related 

works  of  interest,  besides  the  ones  mentioned  above  are:  Akin  (1979, 

1983),  Aronson-Weinberger  (1975,1978),  Beck  (1984),  Butler-Freedman- 
Waltman  (1982),  Christiansen-Fenchel  (1977),  Conley-Fife  (1982),  Hadeler- 
Glas  (1983),  Hadeler-Liberman  (1975),  Losert-Akin  (1983),  Nagylaki  (1977), 
Nagylaki-Crow  (1974),  and,  Shahshahani  (1979). 
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In  system  (3.1),  B(x)/x  and  D(x)/x  are  the  natural,  intrin¬ 
sic  birth  and  death  rates  respectively  of  the  entire  population  x  . 

The  f.  's  can  be  thought  of  as  corresponding  to  fertility  or  viability 
of  the  genotypes  x^.  's.  (Refer  to  Appendix  1  for  the  derivation 
of  the  coefficients  of  B(x)  .)  We  will  make  the  following  assumptions 
on  the  functions  B  and  D  : 


(HI)  B,D  :  cl  (Rj)  -+  cl  (r|)  are  C2  , 

(H2)  B  (0 )  =  D  ( 0 )  =  0  , 

(H3)  B ’ (0)  >  D'(0)  , 

(H4)  B 1 (x)  ,  D'(x)  >0  for  x  e  c1(r|)  , 

(H5)  there  exists  a  unique  K  >  0  ("carrying  capacity")  such  that 


B(K)  =  D(K)  ,  and. 


(H6)  B ' ( K)  <  D'(K) 


where  '  =  =  derivative  with  respect  to  x 


Assumption  (HI)  is  sufficient  to  ensure  existence  and  uniqueness 
of  initial  value  problems  for  all  t  0  and  for  the  invariant  manifold 
and  saddle  node  calculations  to  follow.  Assumption  (H2)  implies  that 
there  can  be  no  birth  nor  death  when  there  is  zero  population.  (H4)  states 
that  the  birth  and  death  functions  are  increasing  functions  of  the  entire  popu¬ 
lation.  (H3)  implies  that  for  small  populations,  the  birth  rate  increases 
more  rapidly  than  the  natural  death  rate.  (H5)  states  that  there  is  a 
carrying  capacity  of  the  environment,  at  which  point,  birth  rate  equals 
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natural  death  rate.  (H6)  says  that  at  the  carrying  capacity  the  birth 
rate  increases  less  rapidly  than  the  natural  death  rate. 

In  Section  3.2,  we  will  show  that  system  (3.1)  defines  a  contin- 

3  2 

uous  semi-flow  <j>  on  cl  (R  )  and  <j>  is  a  smooth  (C  )  local  flow  on 

4“ 

3 

R+  .  Moreover,  all  solutions  <f>(x  ,t)  with  initial  conditions 

x0  =  (x1o,x2o,x3o)  €  cl(R^)  \  {(0,0,0)} 
approach  the  closed,  positively  invariant  simplex 

cl (U)  =  (X  e  cl(R+)  :  x  =  K}  (3.5) 

3 

where  U  =  (X  £  R+  :  x  =  K)  .  In  Section  3.3,  we  discuss  the  rest  points 
of  the  flow  4>  .  A  complete  global  picture  of  the  flow  <{>  on  cl(U)  is 
described  in  Section  3.4  by  means  of  the  "projective"  flow  $  of  4 
on 

cl  (S)  =  {(x-|,x3)  e  cl(R^)  :  X-J  +  x3  <  K}  (3.6) 

where 

2 

S  =  { (x-j  ,x3)  £  R+  :  x-|  +  x^  <  K}  . 

Using  these  results,  we  show  in  Section  3.5  the  main  result  of  the 
chapter:  for  all  X  £  cl  (R^)  ,  4>(X,t)  converges  to  some  rest  point  of 
(p  ,  as  t'  tends  to  +°°.  It  turns  out  that  which  rest  point  <J>(X,t) 
converges  to  is  predictable  and  these  results  are  directly  analogous  to 
the  discrete  case  (see  Chapter  3  of  Roughgarden  (1979)).  Finally, 

Section  3.6  is  a  short  discussion  of  the  results. 
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3.2.  The  Associated  Flow. 


Let  us  define  the  right  hand  side  of  system  (3.1)  to  be 
(0,0,0)  when  X  =  (x.|,x2,x3)  =  (0,0,0)  .  That  is,  define 


by 


where 


F  =  (FrF2,F3)  :  cl(R3)  -  R3 


F(X)  =  (F1(x1 ,x2,x3),F2(x1 ,x  ,x3),F3(x1 ,x2,x3))  , 


F-j  (X) 


=  < 


( u+v ) 


0 


B(x)  - 


X1 

—  D(x) 
x 


if  X  e  cl(R^)  \  {(0,0,0)} 


if  X  =  (0,0,0) 


2uv 


x. 


(u+v)' 


B(x)  -  -f  D(x) 


X 


if  X  e  cl(JT)  \  {(0,0,0)} 


F2(X)  =  { 


0 


if  X  =  (0,0,0) 


f3(x) 


(u+v) 


B(x)  - 


—  D(x) 
x 


if  X  £  cl(r)  \  {(0,0,0)} 


=  { 


0 


if  X  =  (0,0,0) 


where  u,v  and  x  are  given  by  (3,2)  and  (3.3).  The  system  (3.1) 


becomes 


(3.7) 


X  =  F  ( X ) 


(3.8) 
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3 

Clearly,  F  is  locally  Lipschitz.  For  example,  if  X  £  cl(R+)  with 
0  <  |X|  <  r  , 

| F ( X )  -  F(0) |  <  3 1 B ( x )  +  D(x)]  <  C  |x|  <  C'  |X|  , 

for  some  constants  C,C'  >  0  depending  on  r  only.  Also,  F(0)  =  0  . 

We  will  summarize  some  of  the  properties  of  system  (3.8)  in  the 
following  proposition. 

Proposition  3.1 :  System  (3.8)  defines  a  continuous  positive  semi -flow  <p 

on  cl(R+)  which  becomes  a  smooth  local  flow  when  restricted  to  R+  . 

The  simplex  cl (U)  defined  in  (3.5)  is  closed,  positively  invariant 

under  4>  and  is  globally  stable  on  cl(R+)  \  {(0,0,0)}  ,  Furthermore, 

the  vector  field  F  points  into  R?  on  b(R^)  \  cl(H  UH  )  ,  H 

X1  x3  X1 

and  H  are  positively  invariant  and  (0,0,0)  is  a  repelling  rest 
x3 

point  of  <J)  . 

Proof:  The  definition  of  the  flow  <f>  by  (3.8)  and  its  smoothness  prop- 

3 

erties  are  clear.  It  is  also  easy  to  see  that  F  points  into  R+  on 

b(R^)  \  cl(H  UH  )  .  For  example  for  X  £  H  ,  we  have,  x,  = x_ =  0  , 

+  X|  Xj  Xp  o 

so  that  F^(X)  =  F^(X)  =  1  B(X)  >  0  .  Consequently  F(X)  points  into 

R  .  From  this  and  from  the  fact  that  H  and  H  are  invariant,  it 
+  x-,  x3 

3 

follows  from  Proposition  2.1  that  cl(R  )  is  positively  invariant. 
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On  the  other  hand,  since 

x  =  B(x)  -  D(x)  (3.9) 

the  assumptions  (HI)  -  (H6)  in  Section  3.1  imply  that  x  =  K  is  a  sink 
for  system  (3.9)  and  for  all  x(0)  =  xQ  £  rJ  ,  we  have,  0  <  x(t)  +  K  , 

as  t  -*  +  °°.  This  shows  that  cl(U)  is  globally  stable  over 

cl(R^)  \  {(0,0,0)}  . 

q.e.d. 

Remark:  Let  V  :  cl(R2)  \  {(0,0,0)}  ->  be  defined  by  V(X)  =  x  . 

Then  for  all  K'  e  r|  ,  V'^K')  is  a  closed  simplex  and  for  all  K'  f  K  , 

the  vector  field  F  on  V~^(K')  points  towards  V~"*(K)  =  cl  ( U )  . 


3.3.  The  Rest  Point  Set. 

By  Proposition  3.1,  all  rest  points  on  <J>  other  than  (0,0,0) 
must  lie  on  cl(U)  .  Hence  (K,0,0)  and  (0,0, K)  are  the  only  rest 
points  in  b(R2)  besides  (0,0,0)  . 


Proposition  3.2:  Define  the  curve  H  :  r|  +  U  by 


H(c)  = 


c2K 


2cK 


K 


,(l+c)  (1+c)  (1+c)1 


(3.10) 


for  all  c  c  R 


,1 


1 

+ 


Then  all  rest  points  of  <j>  jn_  R+  must  lie  on 


h(r;)  . 
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★  ★  ★  ★  ^ 

Proof:  Let  X  =  (x^x^jX^)  e  R+  be  a  rest  point  for  <J>  and  let 
* 

*  *  *  * 

c  =  — Sr  .  Then  c  >  0  and  x  =  x,  +  x0  +  x0  =  K  .  Let 
2  x!  123 


*  *  ]  * 
u  =  f1x1  +  £  f2x2 


*  *  ]  * 
V  "  ^3X3  +  ~1  ^2X2 


(3.11) 


According  to  (3.1),  we  have, 


u 


*2 


*  *  p 

(u  +v  r 


B(K)  - 


x 


-l  D(K) 


=  0 


*k  'k 

2*  V*  2  B(K) 
(u  +v  )d 

*o 

v 

*  *9 

(u  +v  r 


*jL 

K 


D(K)  =  0 


B  ( K)  - 


■Y  D(K) 


=  0  , 


so  that. 


*9  *  *  *  9 

u  K  =  x1  (u  +v  ) 


(3.12) 


*  *  *  ,  *  *  2 
2u  v  K  =  x2  (u  +  v  ) 


(3.13) 


★  9  *  ,  *  *2 

v  K  =  x^  (u  +  v  ) 


(3.14) 


since  B(K)  =  D(K)  .  Therefore 


(3.13)  (3.13) 

(3.12)  *  (3.14) 
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gives 


4  = 


X1  X3 


or 


x-j  =  c 


2  * 
x3 


•k  2  * 

Now,  x  =  K  implies  (c  +  2c  +  l)x^  =  K  and  hence 


★ 

X 


c2K  ZcK 
(1+c)2  ’  (1+c)2 


as  desired. 


q.e.d. 


Remarks:  (i)  c1(H(r]))  =  H(r])  U{(K,0,0) ,(0,0,K)}  (3.15) 

-  "t"  ~r 

It  will  be  referred  to  as  the  Hardy-Wei nberg  manifold.  It  is  a  closed 
segment  of  a  parabolic  curve  in  cl ( U)  . 


( i  i )  1  i  m  H  ( c ) 

c+0 


(0,0, K)  and  lim  H(c)  =  (K,0,0) 
c->  +°° 


Let 


a 


1 


and 


(3.16) 


Proposition  3.3: 

(i)  If  (a-|,a3)  =  (0,0)  ,  then  H(p|)  is  the  set  of  rest 

3 

points  of  (p  in  R  . 

(ii )  If.  a^  ,  a^  >  0  or  if  a^  ,  a^  <  0  ,  then  <j>  has  exactly 

3 

one  rest  point  in  R+  ,  namely,  H(c)  wi th  c  =  a^/a-j  . 
(iii)  If  a-j  >_  0  ,  a3  £  0  a_nd_  (a^  ,a3)  f  (0,0)  or  if  a^  <  0  , 


40 


a3  —  0  anc*  ( a i  >33)  t  (0,0)  ,  then  <p  has  no  rest  points  in 

^  O 

Proof:  By  the  previous  proposition,  all  rest  points  X  of  c|>  in  R+ 
are  of  the  form  H(c)  given  by  (3.10)  for  some  c  e  ,  According  to 

T 

(3.11), 


u  =  -C-K- ■  ~2  (eg  +  f2) 


(1+c) 


(3.17) 


v*  '  -H  +  f3> 


(1+c) 


Substituting  this  into  (3.12)  gives,  after  some  simplifications. 


(l+c)2(cf-j+f2)2  =  (c2f1+2cf2+f3)2 
or 

cf-j  +  f2  =  cf2  +  f3  (3.18) 

or 

c  =  a^  /  a-|  (3.19) 


provided  a-j  f  0  . 

Thus,  (3.12)  imposes  no  restrictions  on  c  in  case  (i),  it 
requires  c  to  be  a^  /  a-|  in  case  (ii),  and  in  case  (iii)  no  c  e 
will  satisfy  (3.12).  To  complete  the  proof,  all  we  need  is  to  verify 
that  any  c  e  satisfying  (3.12)  also  satisfies  (3.13)  and  (3.14). 

But  this  follows  from  a  straightforward  calculation.  For  example,  (3.13) 
simplifies  to  (3.12)  using  (3.18). 


q.e.d. 
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Remark:  In  case  (i)  we  have  a  1 -dimensional  manifold  of  rest  points.  In 
general,  a  manifold  of  rest  points  (critical  points)  is  referred  to  as  a 
critical  manifold  (see  Conley-Fife  (1982)). 

3.4.  Projective  Flow. 


3 

In  order  to  study  the  dynamics  of  <j>  on  cl  (R  )  ?  we  will 


first  study  the  global  dynamics  of  cp  on  the  globally  stable  positively 
invariant  simplex  cl ( U )  .  This  is  equivalent  to  investigate  the  phase 
portrait  of  the  "projective"  semi-flow  $  of  $  on  cl (S)  (defined 
in  (3.6)).  Here  $  is  defined  by  (3.1)  with  being  replaced  by 
K  -  x-j  -  x3  throughout.  That  is,  we  have. 


x 


(3.20) 


where,  by  (3.2), 


u  =  fK  +  (a-j  +  f)x^  -  fx3 


(3.21) 


v  =  fK  -  fx-j  +  (a^  +  f)x3 


and 


L  =  B(K)  =  D(K) 


(3.22) 


f  =  f2  /  f 


(3.23) 
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By  lifting  the  flow  <j>  back  to  cl  (U)  (that  is,  if  (x-j  (t)  ,Xg(t) ) 
is  a  solution  of  (3.20)  in  cl  ( S )  ,  then  (x-|  (t)  ^(t)  ^(t) )  with 
XgU)  =  K  -  x-j  ( t )  -  x^(t)  will  be  a  solution  of  (3.1)  in  cl  ( U ) )  ,  we 
obtain  the  flow  4  on  cl (U)  .  Clearly,  by  Proposition  3.1  or  by  a 
direct  verification,  the  vector  field  defined  by  system  (3.20)  points 
into  S  on  b(S)  \  { (K,0)  ,(0,K)}  where  b(S)  denotes  the  boundary  of 
S  in  the  x-j-x^  plane  and  that  (K,0)  and  (0,K)  are  the  only  rest 
points  of  $  in  b(S)  .  Hence  S  is  positively  invariant  under  $  . 
According  to  Proposition  3.3, 


(3.24) 


with  0=8^/  a 


is  the  only  rest  point  of  4  in  S  ,  provided  a-j  ^  >0  or  a-j  ,a^  <  0  ; 
and  when  a-|  _>  0  ,  a^  £  0  ,  (a^  ,8^)  ¥  (0,0)  or  when  a-j  £  0  ,  a3  _>  0  , 


(a-jja^)  ¥  (0,0)  ,  i>  has  no  rest  points  in  S  . 


For  the  remainder  of  this  section,  we  will  restrict  our  atten¬ 
tion  to  the  case  when  (a^  ,a^)  ¥  (0,0)  .  The  case  when  (a-pa^)  =  (0,0) 


for  the  "full"  system  (3.1)  will  be  discussed  in  Section  3.5. 

The  variational  matrix  M(x-|  ,Xg)  of  (3.20)  at  an  arbitrary 
point  (x-jjX^)  in  cl(S)  is  given  by: 


2u(fu+a-|V+fv )  1 


2u(a3U+fu+fv ) 


M( X-j  »x^ )  L 


(3.25) 


2v(fu+a-|  v+fv) 


2v(a2U+fu+fv)  1 


(u+v)3 
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Thus,  in  particular. 


M(K,0)  =  t 


f- 


0 


2(a3+f)' 


1 


-1 


(3.26) 


M  ( 0 ,  K )  = 


L 


-1 


2(a1+f) 

f— 


0 


(3.27) 


*  ★ 

and,  if  (x^,x3)  exists, 


"k  k 

M(x-j  ,x3) 


L 

K 


2a3el 


RW 


-  1 


al6l 


a3e3 

RW 


2ale3 


-  1 


(3.28) 


where 


R  =  a-|  +  a^  ,  W  =  f •) f ^  -  ft 


(3.29) 


e,  =  a?  +  fR  ,  e  =  a2  +  fR 
11  ’33 


(3.30) 


Note  that  under  these  notations, 


*  a3  K 

X1  = 

1  R 


T 


X3 


a2  K 


R 


T 


a3KW  *  KW 
’  v  = 


*  * 

and,  u  +  v 


R 


R 


u 


3 


3 


KW 

R 
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Proposition  3.4:  The  linearized  stability  properties  of  the  rest  points 
of  4  i_n  cl(S)  are  given  by  the  following  diagram: 


Parameter  Region 

Rest  Points 

ai 

a3 

(K,0) 

(0,K) 

/  *  *  \ 

(x-j  ,x3) 

+ 

+ 

"  j" 

“  5  ” 

-»+ 

- 

- 

-  5  + 

-  5  + 

”  5  " 

+ 

- 

-  }  + 

* 

- 

+ 

—  9  + 

“  5* 

★ 

+ 

0 

“  5  “ 

-  ,0 

* 

- 

0 

-  9  + 

-  ,0 

★ 

0 

+ 

-  ,0 

"  j" 

* 

0 

- 

-  ,0 

- ,+ 

★ 

where  a^  ( i  =  1,3)  js_  +  (resp.  -,0)  means  that  a.  is  positive  (resp. 
negative,  zero),  under  a  rest  point  means  that  both  eigenvalues  of 

the  variational  matrix  M  at  that  rest  point  are  negative,  means 

that  one  eigenvalue  is  negative  and  the  other  is  positive,  "-,0"  means 
that  one  eigenvalue  is  negative  and  the  other  is  zero,  and,  indi - 

cates  that  the  rest  point  does  not  exist  in  that  parameter  region. 


a  -j  L  - 

Proof:  The  eigenvalues  of  M(K,0)  are:  -  and  -  y-  ,  those  of 


M(0,K)  are: 

a3*- 

f3K 

and 

-  and  those  of 

M(x.|  ,x3)  ,  if  (x1?x3) 

exists,  are: 

ala3L 

and 

-  .  To  complete 

the  proof,  all  we  need 

WK 

to  note  that  if  a-j  ,a3  >  0  ,  then  W  >  0  and  if  a-j  ,a^  <  0  »  then 


i  s 

W  <  0  . 


q.e.d. 


45 


Next,  we  will  investigate  the  phase  portrait  of  $  in  a  neigh¬ 
bourhood  of  each  rest  point. 

Proposition  3.5:  We  have  the  following  cases,  depending  on  the  parameter 
regions : 

(i )  a-j  >  0  ,  a3  >  0  : 

*  *)^ 

(K,0)  and  (0,K)  are  sinks,  (x^ ,Xg)  is  a  saddle; 

(i i )  a^  <  0  ,  a^  <  0  : 

(K,0)  and  (0,K)  are  saddles  whose  stable  manifolds  have  triv¬ 
ial  intersections  with  cl(S)  ,  namely  at  (K,0)  and  (0,K) 

respectively,  and  whose  unstable  manifolds  intersect  cl(S) 

★  ★ 

at  a  whole  separatrix  (out  of  four  for  a  saddle),  (x^  jX^)  j_s 
a  sink; 


(i ii )  a-j  >  0  ,  a  3  <  0  : 

(K,0)  is  a  sink,  (0,K)  is  a  saddle  whose  local  phase  portrait 


is  similar  to  that 

of  case  (i i ) ; 

(iv) 

a-j  <  0  ,  a^  >  0  : 

(0,K)  is  a  sink. 

(K,0)  is  a  saddle 

whose 

local  phase  portrait 

is  similar  to  that 

of  case  (i i ) ; 

(v) 

a-j  >  0  ,  a^  =  0  : 

(K,0)  is  a  sink, 

(0,K)  is  a  saddle 

node: 

it  has  a  neighbour- 

o 

hood  in  R  which  is  the  union  of  one  parabolic  sector,  two 
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hyperbolic  sectors  and  three  separatrices ,  (0,K)  is  stable 

within  the  parabolic  sector,  the  union  of  the  two  hyperbolic 

sectors  covers  a  neighbourhood  of  (0,K)  i_n  cl  (S)  with  the 
separtrix  dividing  the  hyperbolic  sectors  lying  wholly  inside 

cl (S)  ; 

(vi )  a-j  <  0  ,  =  0  : 

( K ,0)  is  a  saddle  whose  local  phase  portrait  is  similar  to 
that  of  case  (i i ) ,  (0,K)  is  a  saddle  node:  it  is  stable  with¬ 

in  the  parabolic  sector  which  covers  a  neighbourhood  of  (0,K) 
in  cl (S)  ; 

( v i  i )  a-j  =  0  ,  a^  >  0  : 

(0,K)  is  a  sink,  (K,0)  is  a  saddle  node  whose  local  phase 
portrait  is  similar  to  that  of  (0,K)  i n  case  (v) ; 

(viii )  a-|  =  0  ,  a^  <  0  : 

(0,K)  is  a  saddle  whose  local  phase  portrait  is  similar  to  that 
of  case  (i i ) ,  (K,0)  is  a  saddle  node  whose  local  phase  por- 

trait  is  similar  to  that  of  (0,K)  in  case  (vi). 


^3”^2 


•  1 


Proof:  When  a-j  <  0  (consequently  (K,0)  is  a  saddle),  ^ 

is  an  eigenvector  correspondi ng  to  the  negative  eigenvalue  -  jr  of 
M ( K , 0 )  .  Thus  the  corresponding  1  -  dimensional  eigenspace  has  slope 

L ^  .  (It  is  vertical  if  2f^  -  f2  =  0  )•  Now  since  the  line 

3  2  jc 

'  ? 

x-,+x3  =  K  has  slope  -1  and  2f  t"1 .0]  »  therefore  the  1  -dimen- 

3  2 
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r 

sional  stable  manifold  W  (K,0)  of  (K,0)  has  trivial  intersection  with 
cl(S)  ,  namely  at  (K,0)  itself  (because  of  the  stable  manifold  theorem. 
Theorem  2.6).  On  the  other  hand  (1,0)  is  an  eigenvector  corresponding 

a-,  L 

to  the  positive  eigenvalue  -  of  M(K,0)  .  Therefore  the  unstable 

manifold  Wu(K,0)  of  (K,0)  has  slope  0  and  has  non-trival  intersection 
with  cl(S),  namely  in  a  whole  separatrix  (out  of  four  for  a  saddle  point), 
due  to  the  positive  invariance  of  cl (S)  and  the  Hartman-Grobman  theorem 
(see  Section  2.9).  Similar  arguments  apply  to  the  case  when  a^  <  0  . 

This  discussion  and  Proposition  3.4  settle  the  cases  (i)-(iv). 

Next  consider  the  case  when  a-j  =  0  .  We  will  use  the  procedure 
described  in  Section  2.10  to  show  that  (K,0)  is  a  saddle  node.  Let 


1 


N  = 


2f  -f 
3  2 


1 


and  Y  =  N'^Z  .  Then  (3.20)  takes  the  form  (2.6).  Moreover,  by  a 

a3L 

strai ghtforward  calculation,  the  n  in  (2.9)  is  equal  to  - 7  .  Since, 

*  2f1r 

fo 

as  noted  above,  -r  i  [-1,0]  ,  the  assertions  made  under  the  cases 

zt3”r2 

(vii)  and  (viii)  are  easily  verified.  The  case  when  a^  =  0  is  similar. 

q.e.d. 
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Knowing  the  local  behavior  of  $  near  each  rest  point,  we  can 
now  describe  the  global  dynamics  of  $  on  cl  (S)  . 

Theorem  3.6:  For  all  X  e  cl(S)  ,  J(X,t)  converges  to  a  rest  point  of 
4>  in  cl (S )  .  More  precisely,  we  have  the  following  cases  depending  on 
the  parameter  regions: 

(i)  a1  >  0  ,  a3  >  0  : 

cl (S)  is  the  disjoint  union  of  the  three  stable  manifolds: 

Ws (K,0)  ,  W(i) * * * S(0,K)  ,  and  Ws(x*,x*)  .  WS(K,0)  and  Ws(0,K) 

are  2-dimensional  ,  and  Ws(x|,xp  is  1 -dimensional  which  acts 

r 

as  a  separatrix  dividing  the  two  regions  of  stability:  W  ( K ,0 ) 
and  Ws(0,K)  of  the  rest  points  (K,0)  and  (0,K)  ; 

(i i )  a-j  <  0  ,  a^  <  0  : 

(x-j  ,Xg)  is  globally  stable  on  cl (S)  \  { ( K ,0 ) , (0 ,K) }  ; 

(iii )  a-j  >  0  ,  a^  <  0  : 

( K ,0)  is  globally  stable  on  cl (S)  \  { (0 ,K) }  ; 

(iv)  a-|  <  0  ,  a^  >  0  : 

(0,K)  is  globally  stable  on  cl (S)  \  {(K,0)}  ; 

•7 

(v)  a-j  >  0  ,  a^  =  0  : 

(K,0)  is  globally  stable  on  cl (S)  \  {(0,K)}  ; 

(vi )  a-|  <  0  ,  a^  =  0  : 
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(0,K)  is  globally  stable  on  cl(S)  \  { ( K ,0) >  ; 

(vi  i )  a-j  =  0  ,  >  0  ; 

(0,K)  is  globally  stable  on  cl(S)  \  { ( K ,0) }  ; 

(vii i )  a-j  =  0  ,  a^  <  0  ; 

(K,0)  is  globally  stable  on  cl(S)\  { (0 ,K) }  . 

Proof:  In  order  to  show  that  for  all  X  £  cl (S)  ,  $(X,t)  converges 

to  a  rest  point  of  $  as  t  tends  to  +°°  ,  it  suffices,  according  to 
the  Poincare-Bendixson  Theorem,  to  show  that  there  is  no  periodic  orbit 
nor  rest  points  connecting  orbit  forming  a  loop  in  cl (S)  .  Using 
Proposition  3.5  and  simple  geometric  arguments,  one  could  easily  show 
that  there  are  no  rest  points  connecting  orbits  in  cl (S)  . 

On  the  other  hand,  since  every  periodic  orbit  must  enclose  a 
rest  point,  there  can  be  no  periodic  orbit  in  cases  (iii )  -  (viii ) . 
(Otherwise,  the  periodic  orbit  will  have  to  enclose  ( K ,0 )  or  (0,K) 
contradicting  the  positive  invariance  of  S  .) 

To  show  that  there  can  be  no  periodic  orbit  enclosing  the  rest 

•j* 

point  (x-j  ,  one  could  use  an  index  argument  in  case  (ii)  (since 

x-j  ,Xg)  is  a  saddle  in  the  case).  As  for  case  (i),  let  L-j  be  the 

★  * 

open  line  segment  joining  ( K ,0 )  and  (x-j  ,Xg)  (that  is,  end  points  are 
excluded)  and  let  l_2  be  the  open  line  segment  joining  (0,K)  and 
(x-j  ,Xg)  .  Also  let  T  denote  the  open  triangle  with  vertices  (K,0)  , 
(x-|  ,x^)  »  and  (K,0)  .  Then  a  strai ghtforward  calculation  shows  that 
the  vector  field  defined  by  system  (3.20)  points  into  S  \  cl(T)  on 


50 


L]  U  I2  •  For  example,  if  (x^  e  L-j  ,  the  dot  product  between  the 
vector  field  at  ( x-j  >x^ )  and  the  outward  normal  (l,l+2c)  of  L-j  is 
given  by: 


u 


.(u+v) 


7 


L  + 


.(u+v) 


7  -  L(l+2c)  =  .2v(cv^)L 

2  V  (u+v)2 


which  is  negative.  Thus,  there  can  be  no  periodic  orbit  enclosing 
*  ★ 

(x^ jX^)  .  Hence,  there  is  no  periodic  orbit  in  cl (S)  . 


The  rest  of  the  assertions  follow  directly  from  the  above  and 
Proposition  3.5. 


q.e.d. 


3.5.  Global  Dynamics. 

We  will  first  describe  the  linearized  stability  properties  and 
the  local  stability  properties  of  the  rest  points  of  (p  . 

Proposition  3.7:  The  linearized  stability  properties  of  the  rest  points 
of  <p  are  given  by  the  following  diagram: 
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Rest  Point 

Parameter  Region 

(K.0,0) 

(0.0, K) 

*  *  ★ 

( X-J  ,*2  >^2 ^ 

al 

a3 

+ 

+ 

"  j"  5“ 

»+ 

- 

- 

-  }  J  + 

—  j  j  + 

~  “  9  "" 

+ 

- 

~  9  “  5  “ 

-  j  j  + 

★ 

- 

+ 

—  ^  5  + 

★ 

+ 

0 

- ,-  ,0 

* 

- 

0 

—  5  5  + 

-,-,o 

* 

0 

+ 

- ,0 

★ 

0 

- 

- ,-  ,0 

—  5  —  5  + 

* 

where  means  that  all  three  eigenvalues  of  the  variational 

matrix  at  that  rest  point  are  negative,  etc.,  as  in  Proposition  3.4. 

Proof:  This  follows  from  Proposition  3.4  and  the  assumption  (H6)  in 
Section  3.1 . 

q.e.d. 

Proposition  3.8:  We  have  the  following  cases,  depending  on  the  parameter 
regions : 

(i )  a-j  >  0  ,  a3  >  0  : 

(K,0,0)  and  (0,0, K)  are  sinks  (i.e.  they  are  point  attractors) 
(x^,x*,x*)  is  a  saddle  with  a  2-dimensional  stable  manifold  and 
a  1-dimensional  unstable  manifold; 
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(ii )  a1  <  0  ,  a3  <  0  : 

★  *  * 

(x^,X2,X2)  is  a  sink,  (K,0,0)  is  a  saddle:  its  2-dimensional 

c 

stable  manifold  W  (K,0,0)  has  trivial  intersection  with 

3 

cl (R , )  ,  namely  at  H  .  ,  and  its  1 -dimensional  unstable  mani- 

xi 

fold  Wu(K,0,0)  intersects  cl(R^)  in  a  whole  branch  (out  of 
two  for  the  1-dimensional  unstable  manifold),  (0,0, K)  is  also 
a  saddle  whose  local  phase  portrait  is  similar  to  that  of 

(K,0,0)  ; 

(iii )  a-j  >  0  ,  a3  <  0  : 

(K,0,0)  is  a  sink,  (0,0, K)  is  a  saddle  whose  local  phase  por¬ 
trait  is  similar  to  that  in  case  (ii); 

( i v )  a-j  <  0  ,  a^  >  0  : 

(0,0, K)  is  a  sink,  (K ,0 ,0)  is  a  saddle  whose  local  phase  por¬ 
trait  is  similar  to  that  in  case  (ii); 

(v)  a-j  >  0  ,  a^  =  0  : 

(K,0,0)  is  a  sink,  (0,0, K)  is  a  saddle  node:  its  strong 

stable  manifold  Wss(0,0,K)  ,  has  trivial  intersection  with 
3  3 

cl  (Rj  ,  namely  at  H  ,  cl(R_J  is  on  the  side  of  the  strong 

x3 

stable  manifold  for  which  (0,0, K)  is  unstable; 

(vi )  a-j  <  0  ,  a^  =  0  : 

(K,0,0)  is  a  saddle  whose  local  phase  portrait  is  similar  to 
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that  i n  case  ( 1  i ) ,  (0,0, K)  is  a  saddle  node:  its  strong 
stable  manifold  has  trivial  intersection  with  cl(R+)  ,  namely 

3 

at  H  ,  cl(R,)  is  on  the  side  of  the  strong  stable  manifold 
x3  + 

for  which  (0,0, K)  is  stable; 

(vii )  a-|  =  0  ,  a^  >  0  : 

(0,0, K)  is  a  sink,  (K,0,0)  is  a  saddle  node  whose  local  phase 
portrait  is  similar  to  that  of  (0,0, K)  in  case  (v) ; 

(viii )  a^  =  0  ,  a^  <  0  : 

(0,0, K)  is  a  saddle  whose  local  phase  portrait  is  similar  to 
that  in  case  (ii ) ,  (K,0,0)  is  a  saddle  node  whose  local  phase 
portrait  is  similar  to  that  of  (0,0, K)  i n  case  (vi ) . 


Proof:  The  variational  matrix  of  (3.1)  at  Xecl(R3)  \ {(0,0,0)}  is  given 

by 


Mn 

M12 

i 

<y-> 

s! — 

M ( x i  ,x3  »x3 )  - 

M21 

M22 

M23 

-M31 

M32 

M33- 

(3.31) 


where 


2f  uv  2 

M,  =  — ! — T  B(x)  +  — 

(u+v)3  (u+v)' 


B'(x) 


X 


M 


12 


f2u(v-u) 

(u+v)J 


X- 


— “  B 1 ( x )  +  -I  D(x)  -  -1  D‘(X) 

( u+v ) ^  x 


B(X)  + 


9 
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Thus , 


2f3u  2  x,  x, 

M13  =  -  —  -  -  B(x)  +  — u-  -7  B '  (x)  +  -i  D(x)  -  4-  D'(x) 
3  (u+v)3  (u+v  )2  x2  x 


M 


2f^v(v-u) 


2uv 


x. 


x 


PI  ’J -  B(x)  +  nr  B  (x)  +  — n-  D(x)  -  — —  D  (x) 

2  (u+v)3  (u+v)2  X7  x 


fp(u-v)  p  i  jv  xl+x^ 

Moo  =  — — -y-  B(x)  +  4^-?  B'(x)  -  -^r  D(x)  -  4  D'(x) 


'22 


(u+v)' 


( u+v )  ‘ 


X 


2f-u(u-v)  p  x?  x? 

Mo.  =  3  3  B(x)  +  ^2UV  2  B 1 (x)  +  4  D(x)  -  4  D‘(x) 


23 


(u+v) 


(u+v )‘ 


2flV  v2 

M^=  *  — ■ — p  B ( x )  +  — 

31  (u+v)3  (u+v) 


7  B 1 (x)  +  -f  D ( x )  -  -1  D'(x) 

X 


x. 


M„„  =  fzV(U  ^  B ( x )  +  vZ  2  B'(x)  +  ~3  D(x)  '  X  D'(x) 


32 


(u+v) 


(u+v)' 


2f9uv  w2  x,+x9  x. 

=  3  --  B ( x )  +  7  B 1  (x)  -  -U4  D ( x )  -  -±  D‘(x) 

(u+v)3  x3  x 


33  '  (u+v)3 


M(K,0,0)  = 


-e  -e  + 


0 


0 


aiL 


aiL 


0 


-e  + 


(f1-2f3)L 

-rpr— 


2f3L 

fiK 


L 

K 


C 


.32) 


where 


e  =  D ' (K)  -  B ' ( K)  >  0 


(3.33) 
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by  (H6)  of  Section  3.1.  Hence  the  eigenvalues  of  M(K,0,0)  are:  -e  , 


a  i  L  I 

-  yy  »  anc*»  ■  X  *  W^en  al  >  0  (consequently,  (K,0,0)  is  a  saddle) 


the  2-dimensional  eigenspace  corresponding  to  the  two  negative  eigenvalues: 

2f. 

i  c  c  nannorl  hw  +  \/cw~  +  r>y'C  ( 1  0  anrl  _ 

2f. 

-  <  D  _  thp  niKitivp  i nvari anrp  nf 


L  ■" 1  3 

-e  and  -  F  is  spanned  by  the  vectors  (1,0,0)  and  (0  ,  — y-  ,  1)  . 


3  3 

Since  — <  0  ,  the  positive  invariance  of  cl(R  )  and  the  stable  man¬ 


ifold  theorem  (Theorem  2.6)  imply  that  the  2-dimensional  stable  manifold 

3 

of  (K,0,0)  intersects  cl(R.)  at  H  only.  Similarly,  one  can  show 

X1 

that  when  a^  <  0  ,  the  2-dimensional  stable  manifold  of  (0,0, K)  inter- 

3 

sects  cl(R  )  at  H  only.  This  and  Proposition  3.7  settle  the  cases 

x3 

(i)  -  (iv). 


Next,  consider  the  case  when  a-j  =  0  .  We  will  use  the  proce¬ 
dure  described  in  Section  2.10  to  show  that  (K,0,0)  is  a  saddle  node. 

Let 


N  = 


-1  1 


0 


0  0 


0 


2f. 


and 
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Y  -  N  .  Then  (3.1)  takes  the  form  (2.6).  Moreover,  a  straightforward 

a^L 

calculation  shows  that  the  n  in  (2.9)  is  equal  to  — .  This 

2f]r 

settles  the  cases  (vii)  and  (viii).  The  case  when  a3  =  0  is  similar. 

q.e.d. 

We  can  now  describe  the  global  dynamics  of  the  flow  <f>  on 

cl(R+)  . 

Theorem  3.9:  For  all  X  e  cl(R  )  ,  <|>(X,t)  converges  to  a  rest  point  of 
<|>  as  t  tends  to  +°°  .  More  precisely,  we  have  the  following  cases 
depending  on  the  parameters  regions: 

(i )  a-j  =  0  ,  a^  =  0  : 

c1(H(r|))  is  the  set  of  rest  points  of  4>  in 
cl(R^)  \  {(0,0,0)}  .  cl(R^)  \  {(0,0,0)}  is  foliated  by  the 
strong  stable  manifolds,  Wss(H(c))  (c  e  [0,+»])  ,  of  these 
rest  points,  where 


WSS(H(c))n  cl(R^)  =  {(x-,,x2,x3)  e  cl(R+)  \  {(0,0,0)} 


1  _  c 

X-|  "t-  — 2 —  ^2  _  cx3  =  ^ 


(3.34) 


is  2-dimensional  for  all  c  e  (0,+°°)  , 


WSS (H(0) )  n  cl(R^)  =  WSS(0,0,K)n  cl(R})  =  H  ,  and. 


Wss(H(+~))  ncl(R^)  =  Wss(K,0,0) n  cI(rJ)  =  Hv  are 


+  ' 


1-dimensional . 
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(ii )  a.|  >  0  ,  >  0  : 

cl(R^)  \  {(0,0,0)}  is  the  disjoint  union  of  Ws(K,0,0)  , 
WS(0,0,K)  and  Ws (x^ ,x3 ,x3)  .  Ws(K,0,0)  and  WS(0,0,K)  are 

c 

3-dimensional  ,  W  (x-|,x2,x3)  is  a  2-dimensional  immersed  sub- 

3 

manifold  of  cl (R  )  \  {(0,0,0)}  forming  a  separatrix  surface 
which  divides  cl (R+)  \  {(0,0,0)}  into  the  two  regions  of  stab- 
i 1 iy:  Ws(K,0,0)  and  Ws(0,0,K)  of  the  rest  points  (K,0,0) 
and  (0,0, K)  ; 


(i i i )  a-|  <  0  ,  a^  <  0  : 


*  ★  ★ 

(x-|,X2>x3)  is  globally  stable  on 
(K,0,0)  is  globally  stable  on  H 


cl(R})  \  cl (H 

X1 

and  (0,0,  K) 


U  H  ) 


x. 


is  globally 


stable  on  H ,  ; 

x3 

(iv)  a-j  >  0  ,  a3  <  0  : 

(K,0,0)  is  globally  stable  on  cl(R^)  \  cl(H%  )  and  (0,0, K) 

+  x3 

is  globally  stable  on  Hw  ; 

x3 


(v)  a-j  <  0  ,  a3  >  0  : 

(0,0, K)  is  globally  stable  on  cl(R^)  \  cl(Hw  )  and  (K,0,0) 
-  +  x-, 

is  globally  stable  on  H  ; 


(vi )  a-|  >  0  ,  a3  =  0  : 

(K,0 ,0)  is  globally  stable  on  cl(R^)  \  cl ( H  )  and  (0,0, K) 

-  i-  x3 

is  globally  stable  on  Hv/  ; 

x3 
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(vii)  a-|  <  0  ,  a3  =  0  : 

(0,0, K)  is  globally  stable  on  cl(R^) 

is  globally  stable  on  H  ; 

X1 


cl(H  )  and  (K,0,0) 

X1 


(vii i ) 


a-j  =  0  ,  a^  >  0  : 

(0,0, K)  is  globally  stable  on 

is  globally  stable  on  H  ; 

X1 


cl(R^)  \  cl  (H  )  and.  (K,0,0) 
X1 


(ix)  a-j  =  0  ,  a^  <  0  : 

(K,0,0)  is  globally  stable  on  cl(R^) 
is  globally  stable  on  H 

x3 


cl (H  )  and  (0,0, K) 
X3 


Proof:  Case  (i)  is  discussed  in  Freedman-Wal tman  (1978).  It  was  shown 

there  that  if  we  let  u  =  u(0)  ,  v  =  v(0)  and  c  =  u  /v  then 

o  '  o  v  '  o  o 

u(t)  =  cv(t)  for  all  t  e  R  .  Furthermore,  x-j  (t)  -  c  x^( t)  ->  0  ,  and 
x^(t)  -  2cx3(t)  ->  0  as  t  +  °° .  Hence, 


(x-j  (t)  ,x2(t)  ,x3(t))  + 


2cK  K  \ 

(l+c)2  ’  (1+c)2/ 


as  t  ->  +°°.  (Note  that  the  assumption  "death  is  certain": 


.+<30 


D(x) 


0 


x 


dx  =  +°°  made  in  Freedman-Wal tman  (1978)  is  automatically  satis¬ 


fied  under  the  hypotheses  (HI )  -  (H6)  in  Section  3.1  because 

lim  D^x)  =  —  >  0  .)  The  assertions  made  under  case  (i)  are  just  a 

X++-  x  K 

rephrase  of  these  results  in  geometric  language. 
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We  will  now  restrict  ourselves  to  the  case  when  (a-j,a3)  i 
(0,0)  .  Recall  the  following  properties  of  the  omega-limit  set,  w(X)  , 
of  a  precompact  non-negative  orbit,  0+(X)  ,  of  a  flow  <j>  : 

(a)  w(X)  is  compact,  connected,  invariant  and  chain  recurrent 
(Proposition  2.5)  and 

(b)  cf ( X , t )  approaches  w(X)  as  t  -*  +°°  . 

For  all  X  e  cl(^)  \  {(0,0,0)}  ,  w(X)  c  cl(U)  by  the  Proposition  3.1. 
According  to  the  Theorem  3.6,  the  only  compact,  connected,  invariant,  chain 
recurrent  subsets  of  cl(U)  are  the  singleton  rest  point  sets.  Thus 
4>(X,t)  converges  to  a  rest  point  in  cl  ( U )  . 


The  assertions  made  under  the  cases  (ii)  -  (ix)  follow  from  the 
above  observation  and  Proposition  3.8. 


q.e.d. 


Remarks :  (i)  It  is  clear  that  one  could  weaken  the  assumption  (H3)  in 

Section  3.1  of  strictly  increasing  B  and  D  somewhat  without  changing 
the  statement  of  the  Theorem  3.9  nor  its  proof. 


(ii)  One  could  also  add  in  an  equation  for  the  carrying  capac¬ 
ity  K  of  the  form:  K  =  G(K)  ,  as  in  Beck  (1982).  If  we  impose  enough 
conditions  on  G  so  that  Kq  >  0  is  globally  stable  over  the  positive 
K  axis  then  the  asymptotic  behavior  described  in  Theorem  3.9  remains 
valid  with  K  replaced  by  Kq  . 
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3.6  Discussion. 

The  biological  interpretation  of  the  parameter  regions  in 
Theorem  3.9  are  as  follows: 

(i)  a-j  =  0  ,  a^  =  0  :  selection  values  are  all  equal, 

(ii)  a-|  >  0  ,  a^  >  0  :  heterozygote  inferiority, 

(iii)  a-j  <  0  ,  a^  <  0  :  heterozygote  superiority, 

(iv)  a-j  >  0  ,  a^  <  0  :  incomplete  dominance, 

(v)  a-j  <  0  ,  a^  >  0  :  incomplete  dominance, 

(vi)  a-j  >  0  ,  a^  =  0  :  selection  against  the  dominant  allele 

(vii)  a-j  <  0  ,  a^  =  0  :  selection  against  the  recessive  allele 

(viii)  a-j  =  0  ,  a^  >  0  :  selection  against  the  dominant  allele 

and , 


( i x )  a-j  =  0  ,  a^  <  0  :  selection  against  the  recessive  allele 


a  , 

A  , 
A  , 


a  . 


Thus  the  results  stated  in  Theorem  3.9  are  directly  analogous  to  the  dis¬ 
crete  case  (see  Chapter  3  of  Roughgarden  (1979)).  The  idea  of  conver¬ 
gence  to  rest  points  for  a  one-locus,  n-allele  genetic  model  is  not  new. 
One  could  refer  to  Butler-Freedman-Wal tman  (1982)  and  Hadeler-Glas 
(1983)  for  models  along  similar  lines,  and  to  Losert-Akin  (1984)  for  the 
Fisher-Wright-Haldane  model. 
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The  usefulness  of  using  a  continuous  model  like  system  (3.1) 
lies  in  the  possibility  of  adding  other  populations  into  the  system  in 
the  form  of  differential  equations.  In  the  next  chapter,  the  genotypic 
effects  of  predator-prey  interactions  will  be  studied  using  these  ideas. 


CHAPTER  4 


A  PREDATOR-PREY  MODEL  CONSISTING  OF  THREE  PREY 

GENOTYPES  WITH  FERTILITY/VIABILITY  DIFFERENCES 


In  this  chapter  a  model  of  a  predator-prey  interaction,  where 
the  prey  population  consists  of  three  genotypes  and  is  modelled  by  system 
(3.1),  is  proposed.  Sufficiency  conditions  leading  to  the  evolution  of 
pure  strains  are  given.  This  extends  results  of  Freedman-Wal tman  (1978, 
1982).  As  well,  conditions  are  derived  which  guarantee  the  persistence 
of  all  components  of  the  system. 


4.1  Introduction  and  the  Model. 

In  this  chapter,  the  model  of  the  previous  chapter  is  extended 
to  include  a  predator.  More  specifically,  the  model  is  given  by: 


2  x 

x,  =  — 1 2  B(x)  -  -1  E D ( x )  +  yP-j  ( x ) ] 
(u+v) 

L  =  ■  -2 3--v-o  B(x)  -  [D(x)  +  yP?(x)] 

2  (u+v)2  x  2 

2  x 

X,  =  — '—o  B(x)  -  -r  ID(x)  +  yP,(x)] 

3  (u+v)3  x  3 


y  =  y [-s  +  k  l  -j-  Pi (x)] 
i=l 


(4.1) 


where  x  ,  given  by  (3.3),  is  the  total  prey  population  and  y  is  the 
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predator  population,  u  and  v  are  given  by  (3.2)  and  s  and  k  are 
assumed  to  be  positive. 

The  genetic  component  of  this  model  arises  by  considering  a 
one-locus,  two-allele  problem  as  given  by  system  (3.1).  The  ecological 
component  arises  by  imbedding  an  intermediate  predator-prey  model  with 
continuous  natural  birth  and  death  processes  in  a  model  which  accounts 
for  different  predator  functional  responses  (given  by  P. (x)  (i=l,2,3)) 
by  the  predator  on  the  different  prey  genotypes  x.  (i=l,2,3)  .  For  a 

discussion  of  intermediate  predator-prey  models,  see  Chapter  4  of  Freed¬ 
man  (1980)  and  the  literature  cited  there. 

Besides  the  assumptions  (HI)  -  (H6)  made  in  Section  3.1  on  the 
birth  and  death  rates,  we  make  the  following  assumptions  on  the  predator 
functional  responses: 

(H7 )  Pi  :  cl(Rj)  ->•  c1(r[.)  are  C2  (i=l,2,3)  , 

(H8)  P.(O)  =  0  (i=l  ,2,3) 

(H9)  P'.(x)  >  0  (i=l  ,2,3)  for  all  x  e  cl(Rj)  • 

Assumption  ( H 7 ) is  sufficient  to  ensure  existence  and  unique¬ 
ness  of  initial  value  problems  for  t  _>  0  .  (H8)  implies  that  there  is 
no  predation  in  the  absence  of  the  prey  population.  (H9)  states  that 
the  predator  functional  responses  are  strictly  increasing  as  functions  of 
the  number  of  prey. 

The  remainder  of  this  chapter  is  organized  as  follows:  In 
Section  4.2,  we  construct  the  flow  on  cl(R^)  defined  by  system  (4.1) 


64 


and  list  some  of  its  properties.  The  boundary  rest  points  of  this  flow 
and  their  stability  are  considered  in  Section  4.3.  In  Section  4.4,  condi¬ 
tions  leading  to  the  evolution  of  pure  strains  and  the  existence  of  a 
globally  stable  boundary  rest  point  are  given.  The  classification  of  the 
parameter  values  into  persistence  and  non-persistence  cases  is  discussed 
in  Section  4.5.  We  finish  with  a  short  discussion  of  the  results  in 
Section  4.6. 


4.2  The  Associated  Flow. 

In  this  section  we  will  show  that  system  (4.1)  defines  a 
continuous  semi-flow  on  cl(R+)  which  becomes  a  smooth  (c"*)  local  flow 
when  restricted  to  rJ  .  We  will  also  show  that  it  is  A-dissipative  and 
discuss  its  boundary  invariant  sets. 

As  in  system  (3.1),  system  (4.1)  is  not  defined  when  x  =  0  . 
Let  us  redefine  system  (4.1)  as  follows.  We  define  system  (4.2)  to  be 


X1  = 


x, 

y— 2  B(x)  -  -+  [ D( x )  +  yP-|  (x)J 


(u+v) 


Xo  = 


Xo  = 


2uv-  2  B(x)  -  -J-  [D(x)  +  yP2(x)]  ,  x  >  0 
(u+v) 

2  X-3 

-  V---2-  B(x)  -  ~Y  [D(x)  +  yP3 ( x) ] 


(u+v) 


3  x. 

y  =  y[-s  +  k  l  —  p. (x)] 
i=l  x  1 


(4.2) 


Xi  =0  (i  =  1,2,3) 


,  x  =  0 
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y  =  -sy 


or 


x,  =  F.(x,,x9Jxvy)  (i  =  1,2,3) 


ilAl,A2,A3 


Let 


F  =  ( F, ,F  ,F,,F  )  . 


1  ’  2’  3 5  4 


(4.3) 


Clearly,  F.  e  Lip  (i  = 1 ,••• ,4)  on  cl(R^)  .  For  example. 


consider  the  two  points  (x^x^x^y)  and  (0,0, 0,y)  where 
x  =  X1  +  x2  +  x3  >  0  .  Then 

3 

|G(x-,  ,x  ,x  ,y)  -  G(0,0,0,y)  |  <  s  |y-y|  +  ky  l  P.(x) 

i=l  1 

3 

1  s  I y-y I  +  ky  I  P!(z.)x 

i=l  1  1 

for  some  0  £  z.  £  x  (i  =  1,2,3)  and  hence  G  e  Lip  . 

Proposi tion  4.1:  System  (4.2)  defines  a  continuous  semi-flow  on  cl(R^) 

1  4 

which  becomes  a  C  local  flow  when  restricted  to  R+  .  It  is  A  -  d i s s i - 
pative  where 

A  =  { ( x-j , x2 , x^ ,y )  e  cI(rJ)  :  x  £  K  and_  kx  +  y  £  j  (M+sK)}  (4.4) 


and  M  =  max  (B(x)  -  D(x)  :  x  €  [0 , K] } .  The  edges  and  faces  H 

- — ~  ■■  mmf.  —  ■  ■■  -  ■  -  A 


5 
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H 


x,  >y  ’  xvy  ’  Hx,  ,x9,x  are  1nvari‘ant  under  the  flow  and  the  vector  field 
o  12  3 

F  »  defined  in  (4.3),  points  into  R4  on  b(R4)  \  cH H  II  H  n  h  i 

- -  +  —  +  x-,,y  x2,y  xrx2,x3;* 


Proof:  Using  Proposition  (2.1)  it  is  easily  seen  that  cl(R4)  is  posi¬ 
tively  invariant.  System  (4.2)  is  A  -  dissipative  follows  from 

3  x. 

x  =  B(x)  -  D(x)  -  y  l  ~~~  P..(x)  <  B(x)  -  D(x) 

i=l  x  1 

k x  +  y  =  k [ B ( x )  -  D(x)]  -  sy 


and  Proposition  2.2.  The  rest  of  the  assertions  can  be  easily  verified, 
similar  to  the  proof  of  Proposition  3.1. 

q.e.d. 


Remarks:  (i)  The  invariance  of  H  can  be  interpreted  as:  if  there 

X1 

were  no  "a"  gamete  nor  predator  to  start  with,  then  there  would  not  be  any 
"a"  gamete  nor  predator  for  any  other  time.  The  invariance  of  the  other 
boundary  invariant  sets  can  be  interpreted  similarly. 

(ii)  The  global  dynamics  of  system  (4.1)  on  the  invariant  set 

cl(H  )  is  given  by  Theorem  3.9.  Its  dynamics  on  cl(H  ) 

x1’x2,x3  x.,y 

(i  =  l,3)  is  given  by  the  intermediate  predator-prey  model: 

xi  =  B(xi)  -  D(xi)  -  yPi(xi) 


y  =  y[-s  +  kP^  (xi )] 


(4.5) 
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4.3  Boundary  Rest  Points  and  Their  Stability. 


In  this  section  those  rest  points  of  system  (4.1)  lying  on 
b(R^)  and  their  stability  properties  will  be  considered.  Clearly, 

<  0  and 


E  (0, 0,0,0)  is  an  unstable  rest  point  (because  y 


x=0 

y>0 


0<x«l 

y=0 


>  0  ).  In  fact,  we  have  the  following  proposition. 


Proposition  4.2:  Eq  7  w(x-|  ,x2,x3,y)  for  any  (x-|,x2,x3,y)  e  . 

4 

Proof:  Since  in  a  neighbourhood  of  EQ  in  cl(R+)  ,  we  have, 
x  £  B(x)  -  D(x)  -yP(x)  £  k^x 

y  £  y[-s+kP(x)l  £  -k2y  for  some  k-j  ,  k^  >  0  , 

and  the  flow  on  cl(Hy)  is  given  by  y  =  -sy  ,  where 

P(x)  =  max  (Pi(x)  :  1  £  i  £  3}  ,  the  result  follows  immediately  from 

Proposition  2.15. 

q.e.d. 

Corollary  4.3:  Let  (0,0, 0,y)  e  cl(H  )  .  Then  (0,0, 0,y)  i  w(x] ,x2#x3,y) 

A 

for  any  (x-|,x2,x3,y)  e  R+  . 

Proof:  Suppose  not,  since  w(x] ,x2 ,x3,y)  is  compact,  connected  and  invar¬ 
iant,  the  global  stability  of  Eq  on  cl (H  )  implies  that 
Eq  e  w(x1  ,x2,x3,y)  contradicting  Proposition  4.2. 


q.e.d. 
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Clearly,  E^(K, 0,0,0)  (resp.  EgCOjOjKjO))  is  the  only  rest 

point  in  H  (resp.  H  )  . 

X1  x3 

Let 

b.  =  -s  +  kPi(K)  (i  =  1,2,3)  .  (4.6) 

A  necessary  and  sufficient  condition  for  H  (resp.  H  )  to  contain 

x-j  »y  x ^  ?y 

a  rest  point  is  b-j  (resp.  b^)  >  0  .  Moreover,  (H9)  in  Section  4.1 
implies  that  these  rest  points,  if  they  exist,  are  unique.  Let  us  denote 
them  by  E^x^ ,0,0, y-j )  and  E4(0,0,x3,y3)  . 

Rest  points  lying  in  H  are  of  the  form 

X-j  5X2? x^ 

★  ★  ★ 

Eg(x^ jX^jX  ,0)  =  (H(c),0)  where  H(c)  is  given  by  (3.10).  According  to 
Proposition  3.3,  in  the  case  when  a-j  =  a3  =  0  ,  (H(c),0)  is  a  rest 
point  for  all  c  e  (0,+°°)  .  Hence  there  is  a  curve  of  rest  points.  In 
the  other  cases,  E^  exists  if  and  only  if  a-j a3  >  0  ;  and  if  it  exists, 
then  it  is  unique  and  is  given  by  (H(c),0)  with  c  =  a3/a^  . 

The  linearized  stability  of  E^  (h=l,***,5)  is  governed  by 
the  sign  of  the  real  parts  of  the  eigenvalues  of  the  variational  matrix 
M  =  [M-jjl  1  <1-  j<4  of  system  (4.1)  evaluated  at  Eh  ,  where 

9F. 

Mij  =  Sz1  ’  Zi  =  xi  ^  =  lj2’3)  »  and’  z4  =  y  • 

Recall  the  definition  of  a.  (i  =  l,3)  ,  L  ,  and,  e  from 
(3.16),  (3.22),  and,  (3.33)  in  Chapter  3. 


Clearly, 


69 


M(E1)  = 


-e  -e  + 


0 


0 


0 


V 


f  1  K 


0 


0 


-e  + 


(fl"2f3)L 

“V 


2f3L 

flK 


K 


0 


■P-j  ( K) 


0 


0 


(4.7) 


alL 

M(E-j)  has  eigenvalues:  -e  ,  - 

L  1 


,  and,  b-j  .  Since  two  of 


these,  -e  and  -  are  negative,  E-|  has  a  stable  manifold  of  at 
least  two  dimensions.  The  actual  dimensions  of  W  (E^)  and  W  (E]) 
depend  on  the  signs  of  a-j  and  b-j  .  Analogous  statements  can  be  made 


for  E2  . 


M(E3)  is  given  by 


M-j-j  =  B'  (x-j )  -  D'U-j)  -  y-|  P-,  (x-j ) 


M12  =  B '  (X] )  -  D '  ( x-j )  - 


f0B(x, )  D(x, ) 


2“^!  7  ,  ^1'  ~  n./77  N  .  1  1  1 


fixi  xi 


+  — —  -  y^fxp  + 


X- 


H13  =  B*  (x-j )  -  D'fxp  - 


2f,B(x.)  D(x.)  y/^x,) 

-^-  +  -^--7^)  +  J=LJ- 

flXl  X1  X1 


M14  =  -P-j  (x-j )  ,  M21  =  0  ,  M22 


f2B(x1)  D(x1)  y-,  P2 (x-j ) 

flxl  X1  X1 


2f.B(x1) 

M23  7=  ’  M24  "  0  ’  M31  "  °  ’  M32  '  0  ’ 

Vl 
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M 


DU-,)  +  y-|p3(xi ) 


33 


,  M34  =  0  ,  M41  =  ky-j  P-j  (x] ) 


—  r  —  Pi (xi )  po(xi ^ i 

M42  =  kyl  [P^V  -  +  “"=  ]  ’ 

X1  X1 


M 


43 


-  r  -  W  p3(xl)l 
=  ky1  [p](x-,)  -  L  +  4-J-j  > 


M44  =  0 


(4.8) 


After  interchaning  the  second  and  fourth  row  and  columns  respectively 
of  M(E3)  and  calling  the  resulting  matrix  M  ,  we  see  that  : 

^41  =  ^42  =  ^43  =  °*  Analysis  this  matrix  shows  that  the  eigenvalues 


corresponding  to  vectors  lying  in  H 


x-,  ,y 


are  such  that  the  sign  of  their 


real  parts  is  the  same  as  the  sign  of  .  The  other  two  eigenvalues  of 
M(E3)  are  given  by  and  .  Clearly  <  0  but  has 

undetermined  sign.  Analogous  statements  can  be  made  regarding  E4  .  For 

future  reference,  we  denote  1^22 ( ^3 )  anc*  ^22^4^  by  ^l  anc* 
respectively.  That  is,  we  define 


d.  =  J- 

i  — 

x . 


f„B(x.) 


v  -  -  y.P2(x.) 


(i  =  1.3) 


(4.9) 


^  ^ 

The  stability  of  E5(x^  jX^x^O)  on  Hx  ^x  ^x  is  given  by 

*12  3 

Proposition  3.7.  To  compute  the  stability  of  E,-  in  the  y- direction 
we  note  that 


,  k  v 
=  "s  K  L 
K  i=l 


x. 


Pi(K) 


(4.10) 


71 


is  an  eigenvalue  of  M(E^)  corresponding  to  eigenvectors  in  y-  direc¬ 
tion,  and  hence  is  asymptoti cal ly  stable  or  unstable  in  the  y- 
direction  as  d  is  negative  or  positive.  The  case  when  f-j  =  =  f^ 

is  of  special  interest.  For  each  E^  =  (H(c),0)  (c  e  (0,+°°))  ,  M(E^) 
has  one  zero  and  two  negative  eigenvalues  corresponding  to  eigenvectors 
on  H  .As  for  the  y- direction,  we  have  the  following  propo- 

X-|  » X£  J  x^ 

si ti on. 

Proposition  4.4:  Let  f^  =  f^  =  fg  and  b^ ,bg  >  0  .  Jjf 

P2(K)  >  |  -  /  (P-j  (K)  -  |)(P3(K)  -  f)  (4.11) 

then  for  all  (xi0>x2o’y3o’-yo^  e  R+  ’  y(t)  >  0  ' 

t-H-00 

Proof:  Let  hi  =  P..(K)  -  |  (i  =  1,2,3)  .-  bi  >0  ( i  =  1,3)  implies  that 

h .  >  0  ( i  =  1,3)  .  Moreover,  h2  >  -  /h.|h3  by  (4.11).  According  to 

Lemma  4.5  (see  below),  we  have, 

d  =  -s  +  -k-^  [c2Pn(K)  +  2cP9(K)  +  P,(K)j 

(1+c)2  1  2  2 

=  — — 9-  [c2h-,  +  2ch9  +  h»]  >  0 
(1+c)2  1  2  2 

for  all  c  e  [0,+°°]  ,  where  d  is  defined  in  (4.10).  Hence  the  Hardy- 

Weinberg  manifold,  c!(H(r|))  is  normally  hyperbolic.  Suppose  that 

lim  y(t)  =  0  .  Using  an  argument  similar  to  the  Butler-McGehee  lemma 
t-H"°° 
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(Lemma  2.14)  and  its  extension  (Proposition  2.16),  one  shows  that 
w  =  w(xn  ,x9  ,x~  ,y  )  must  intersect  cl(H  )  at  a  point  other 

10  l.0  oO  0  X-| 


than  c1(H(r}))  .  The  global  dynamics  of  system  (4.1)  on  cl(H 


Xi » x2 ,  x3 


) 


as  given  in  case  (i)  of  Theorem  3.9  implies  that  Eq  e  w  or  w  is 
unbounded.  In  either  case  we  have  a  contradiction. 


q.e.d. 


Remark:  For  b-j , b^  >0  ,  (4.11)  is  satisfied  if  b2  >  0  . 

Lemma  4.5:  Let  h^h^sh^  e  R  .  Then 

1(c)  =  — - — o  [c2h,  +  2ch0  +  hJ  >0  for  all  c  £  [0 ,-h=°3 

(1+c)2  1  2  3 

if  and  only  if  h-|,h3  >  0  and  h2  >  -  /  h^h^  . 

Proof:  Since  1(0)  =  h3  and  I(+°°)  =  h-|  ,  the  conditions  h-j  >  0  and 
h^  >  0  are  clearly  necessary.  Now  for  c  £  (0,+°°)  ,  1(c)  >0  if  and 
only  if  c2h-|  +  2ch2  +  h3  >  0  ,  and, 

c2h-!  +  2ch2  +  h3  =  (c  /TTj”  -  /hp2  +  2c(h2  +  /  h-^  ). 

This  completes  the  proof  of  the  lemma. 


q.e.d. 
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4.4  Evolution  to  Pure  Strains. 

In  this  section,  criteria  will  be  given  which  guarantee  the 

disappearance  of  one  of  the  gamete  types.  As  well,  sufficiency  condi- 

4 

tions  ensuring  or  E^  to  be  globally  stable  over  R+  will  be 
considered. 


From  the  definition  of  u  and  v  in  (3.2), 


.f^u  +  f0v 


u  =  u  [-j - -y~  B(x)  -  1  D ( x )]  -  ^  [f^fx)  +  ^  f2x2P2(x)] 


(4.12) 


v  =  v  [ -  2 

L  (u+v)2 


3  2  B(x)  -  ^  D(x)]  -  [f3x3P3(x)  +  1  f2x2P2(x)]  . 


Therefore, 


•  • 


a, u-a^v 


v  u_  J 1 u  “  3 v  D  /  \  y  r  1 


v  u 


(u+v) 


2  B(x)  -  ^V[iflf2(P2(x)_Pl(x))xlX2 


+  f1f3(P3(x)-P1(x))x1x3  +  1  f2f3(P3(x)  -  P2(x))x2x31 


Theorem  4, 

.6:  Assume  that 

(i) 

fl  f2  f 3 

9 

(ii) 

P-|  (x)  £  P2(x) 

<  P3(x) 

for  all  x  e  [0 , K]  , 

(iii) 

-s  +  kP-j  (K)  > 

0  ,  and 

9 

(iv) 

there  exist 

kl ,k3  -  0 

not  both  zero  such 

that 

P2(x)  -  P,(x) 

>  k-j  x  and 

p3(x)  -  P2(x)  £  k3x 

for  all 

.  (4.13) 


x  €  [0,K]  . 
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Then  for  all 


(xlo’x2o’x3o’yo)£R+ 


x2(t) 


x3(t) 


-*■  0  as  t  -*■  +00 


Proof:  The  proof  proceeds  in  four  steps. 

Step  1 .  1 im  y(t)  >  0  .  (ii)  and  (iii)  implies  b.  >  0 

t*+x>  1 

(i  -  1,2,3).  Hence  this  follows  from  Proposition  4.4. 


Step  2.  One  can  assume  that  x(t)  £  K  for  all  t  £  0  .  This 

is  because  system  (4.1)  is  A  -  dissipative,  where  A  is  given 

in  (4.4),  and  that  there  are  no  invariant  sets  in 

{ (xi ,x2,x3,y)  e  cl(R^)  :  x  =  K}  other  than  those  that  lie  in 

cl  (H  )  . 

x-|,x2,x3' 

A  A  A  A  /\ 

Step  3.  Points  of  the  form  (x-j  jX^x^y)  where  x^  >  0  ,  and, 
x-j  =  0  or  x^  =  0  are  not  in  w(x]0’x20»x30,x40^  *  Consider 

/\  /\  A  A  A  A 

the  case  when  x-j  =  0  and  x2  >  0  .  Then  F-|  (x-j  ,x-|  jX^y^)  >  0 
and  we  can  apply  Proposition  2.13.  Similarly  we  can  deal  with 
the  case  when  x^  =  0  and  x^  >  0  . 

Step  4.  Since  a-j  =  a^  =  0  ,  according  to  (4.13),  we  have. 


v 

v 


u  1  "  uv  [\  klxlx2  +  (ki+k3)xix3  +  2  k3X2x3] 
i-  lyV  +  T  x3y]  1  !klXly  +  k3X3yl 


5 


since  0  <  u,v  <  K  . 
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Case  1:  >  0  .  Then,  by  integrating. 


v 

v ( t)  <  ~  u(t)  exp 
uo 


x-|  (T)y(i)  dij 


(4.14) 


If  the  last  integral  diverges,  then  lim  v(t)  =  0  .  Otherwise, 

t-Ho° 

( Xi y )  ( t )  is  in  iJ  [0 ,+°°)  .  Since  ^  [x-|y(t)l  is  bounded,  (x^y)(t) 

is  uniformly  continuous.  As  uniformly  continuous  lJ[0,+<»)  functions 

have  limit  zero  when  t  +  +°°  ,  this  shows  that  lim  (x,y)(t)  =  0  . 

t-*+°° 

lim  y(t)  >  0  implies  lim  x,  (t)  =  0  and  hence  lim  x?(t)  =  0  ,  by 

t->+oo  t">+00  t^+°°  ^ 

Step  3.  Thus  lim  u(t)  =  0  and  consequently  lim  v(t)  =  0  ,  by  (4.14). 

t->+°°  t~++°° 

But  this  contradicts  Corollary  4.3. 


Case  2:  >  0  .  This  is  similar  to  Case  1. 


q.e.d. 


Remark:  The  case  when  P^  =  P^  find  k^  >  0  was  studied  in  Freedman- 
Waltman  (1978)  and  the  case  when  P-j  =  P^  and  k^  >  0  was  studied  in 
Freedman-Wal tman  (1982). 


Corollary  4.7:  If,  in  addition,  we  assume  that  exists  and  is  glob- 
ally  stable  (resp.  globally  exponentially  stable)  on  Hw  ,  then  it  is 

xi  >y 

globally  attracting  (resp.  globally  stable)  on  R+  . 

Proof:  The  global  attractivity  of  E^  follows  from  Theorem  4.6, 

1 im  y(t)  >  0  and  a  chain  recurrence  argument.  In  the  case  when  E^ 

t^-fo° 
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is  exponentially  stable  on  H  (i.e.  the  two  eigenvalues  of  M(E_) 

x-|  »y  o 

corresponding  to  eigenvectors  on  H  have  negative  real  parts),  the 

X-j 

s  s  c 

dimension  of  W  (Eg)  is  at  least  three.  Therefore  W  (Eg)  must  be  at 

most  1 -dimensional  and,  according  the  Theorem  4.6,  Eg  is  stable  on 

Wc(Eg)  n  rJ  .  The  stability  of  Eg  on  then  follows  from  the  Palis- 

Takens  linearization  theorem  (Theorem  2.10). 

q.e.d. 

Remark:  Conditions  for  E0  to  be  globally  stable  on  H  w  can  be 

o  x-j  >y 

found  in  Hsu  (1978)  and  Cheng-Hsu-Lin  (1981). 


Theorem  4.8:  Assume  that 

(i)  fi  £  fg  _>  fg  and  they  are  not  all  equal,  and, 

(ii)  P-j  (x)  £  Pg(x)  £  Pg(x)  for  all  x  e  [0 , K' ]  ,  for  some  K'  >  K  . 

Then  for  all  (xi0,x2o’x3o’yo^  €  ’  x2^)»x3^  ^  0  —  *  *  +0°* 


Proof:  Again,  since  system  (4.1)  is  A  -  dissipative,  we  can  assume  that 
x(t)<K'  for  all  t>0.  Since  a]  >  0  and  a3  <  0  ,  from  (4. 1 3) , 
we  have 


v 

V 


anu  -  a0v 
J - l-  B  ( x ) 

(u+vr 


(4.15) 


If  a-j  >  0  (and  a3  £  0),  then  from  (4.15), 


a-j  u 
(u+v)2 


v 

v 


B(x) 
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Integration  yields: 


v  r  rt  a,u(T) 

v(t)  < -p  u(t)  exp  [-  - - - ~ 

o  ^0  (u(t)+v(t)j 


B(x(t) )  dxj  . 


If  1 im  u(t)  -  0  ,  then  1 im  v(t)  =  0  which  in  turn  implies 

t-*-+°°  t^-+°° 

1 im  x(t)  =  0  contradicting  Corollary  4.3.  Therefore  1 im  u(t)  =  u  >  0 
■£->-+00  -£->+00 

This  implies  the  integral 


■+00  a-j  u(t) 

- j  B(x(t))  dx  =  -h»  and  hence  lim  v(t)  =  0 

0  (u(t)+v(t))^  t-H°° 


Suppose  now  that  a-j  =  0  and  hence  a^  <  0  .  Then  from  (4.15) 


*  S  ,  a3V  p/  i 

7  "  IT  -  ” - tB(x) 


(u+v)1 


and 


v(t)  <  -p  u(t)  exp  [  [ 
o  LJ0 


a  v(t) 

- P  B(x(t) )  dx 

0  ( u ( t )+v ( t ) ) ^ 


If  this  last  integral  diverges,  then  again  lim  v(t)  =  0  .  Suppose  not, 

t-H-co 

i.e.  suppose  z(t)  =  - - «  B(x(t))  is  in  L*1  10,+°°)  .  By 

(u(t)+v(t)T 

Corollary  4.3,  x(t)  is  bounded  away  from  zero.  Therefore,  p:  z(t) 
is  bounded  and  so  z(t)  is  uniformly  continuous.  Thus  lim  z(t)  =  0  . 

t-*+°o 

Since  1 im  B ( x ( t ) )  >  0  ,  we  have,  lim  v(t)  =  0  . 

t-H-oo  t->-+0° 


q.e.d. 
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Corollary  4.9:  If,  in  addition,  we  assume  E^  exists  and  is  globally 

stable  (resp.  globally  exponentially  stable)  on  H..  ..  ,  then  it  is 

x-j  ,y 

globally  attracting  (resp.  globally  stable)  on  R+  . 


Proof:  Let  (x,  ,x0  .x0  ,y  )  e  R,  .  Then  w  =  w(x,  ,x0  ,x~  ,y)  c 
-  '  lo  2o  3o  Jo'  +  v  lo  2o  3o  ' 

cl(Hv  w)  ,  by  Theorem  4.8.  We  first  show  that  i  w  .  Suppose  not, 


x-,  ,y 


C  Q 

i.e.  suppose  E-jew,  then  by  Proposition  2.16,  w  n  W  (E-j)  contains 
a  point  other  than  E,  .  Since  WCS(E-,)  n  cl(R^)  c  cl(H  )  ,  this 

I  I  ' 

means  that  w  contains  a  point  in  cl (H  )  other  than  E,  ,  which  is  a 

X1  1 

contradiction.  Therefore  E,  l  w  .  The  global  dynamics  on  cl(H  ) 

I  x-j  9  y 

and  the  chain  recurrence  of  w  imply  that  w  =  (E^l  .  Hence  E^  is 

4 

globally  attracting  on  R+  . 


The  proof  of  the  stability  of  E3  on  R+  when  it  is  exponen¬ 


tially  stable  on  H 
Corollary  4.7. 


x]  ,y 


is  similar  to  that  given  in  the  proof  of 


g.e.d, 


Remark:  In  Theorems  4.6  and  4.8,  one  could  interchange  the  roles  of 
the  gamete  types  A  and  a  and  obtain  conditions  for  which  x-j  ( t)  ,x2(t)  -+  0 
as  t  -*•  +«>  .  The  same  observation  applies  to  Corollaries  4.7  and  4.9. 
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4.5  Persistence  and  Non-Persistence  Results. 

In  this  section,  conditions  leading  to  the  persistence  and  non¬ 
persistence  of  system  (4.1)  will  be  considered. 

In  the  case  when  all  selective  values  are  equal,  we  have  the 
following  theorem  for  persistence. 


Theorem  4.10:  Let  f-|  =  f^  =  f^  .  Assume  that  both  and  exi st 

and  are  globally  stable  on  H.#  and  H,  w  respectively.  If 

X-|  ,y  x^  , y 

d-j  ,  d^  >  0  and  (4.11)  holds,  then  system  (4.1)  is  persistent. 


Proof:  Let  (xi0»x20’x3o,y^  6  R+  and  w  =  w^xlo,x2o,X3o,y^  *  The 
existence  of  E^  and  E^  implies  b-j  ,  b^  >  0  .  Proposition  4.5  shows 


that  w  n  cl  ( H  )  =  cf>  .  Next  we  show  that  E.  /  w  .  Since  d,  >  0  , 

X1 ,x2 ,x3  J 

dim  Wu ( E, )  =  1  and  WCS(EJ  n  cl(R^)  c  cl(H  )  .  By  proposition  2.16, 

i  1  +  xi  ,y 


n  W05^)  contains  a  point  other  than  E^  and  hence  w  contains  a 

int  in  cl (H  )  other  than  Eo  .  But  this  is  impossible  because  of 

x-j  ,y  J 


the  global  dynamics  of  system  (4.1)  on  cl(H  )  and  that  E  ,  E,  l  w  . 

X-|  9J  V  \ 


Similarly,  we  can  show  that  E^  it  w  .  Knowing  these  facts,  it  is  then 

.  4. 

easy  to  show  that  w  n  b(R+)  =  <J>  . 

q.e.d. 


Remark:  Since  B(x. )  -  D(x.)  -  y^P^x^)  -  0  ,  in  the  case  when 
f 1  =  f 2  =  f 3  »  the  condition  d.  >  0  ( i  =  1 ,3)  simplifies  to 
P2(xi)  >  P-j^.)  (i  =  1,3)  . 
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In  the  case  when  the  selective  values  are  different,  we  have 
the  following  theorem. 


Theorem  4.11 .  Assume  that  (resp.  E^)  is  globally  exponentially 


stable  on  H 


x-j  ,y 


(resp.  on  H  )  if  it  exists.  We  have  the  foil 


x3’y 


owing 


table  for  persistence  and  non-persistence  of  system  (4.1): 


ai 

a3 

bi 

b3 

dl 

d3 

d 

P  /  NP 

la 

- 

- 

- 

- 

- 

NP 

lb 

- 

- 

- 

- 

+ 

P 

2a 

- 

- 

- 

+ 

- 

NP 

2b 

- 

- 

- 

+ 

- 

NP 

2c 

- 

- 

- 

+ 

+ 

+ 

P 

3a 

- 

- 

+ 

- 

- 

NP 

3b 

- 

- 

+ 

- 

- 

NP 

3c 

- 

- 

+ 

- 

+ 

+ 

P 

4a 

- 

- 

+ 

+ 

- 

NP 

4b 

- 

- 

+ 

+ 

- 

NP 

4c 

- 

- 

+ 

+ 

.  - 

NP 

4d 

- 

- 

+ 

+ 

+ 

+ 

+ 

P 

5 

- 

+ 

- 

- 

NP 

6a 

- 

+ 

- 

+ 

- 

NP 

6b 

- 

+ 

- 

+ 

+ 

P 

7 

- 

+ 

+ 

- 

NP 

8a 

- 

+ 

+ 

+ 

- 

NP 
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8b 

- 

+ 

+ 

+ 

- 

NP 

8c 

- 

+ 

+ 

+ 

+ 

+ 

P 

9 

+ 

- 

- 

- 

NP 

10 

+ 

- 

- 

+ 

NP 

11a 

+ 

- 

+ 

- 

- 

NP 

lib 

+ 

- 

+ 

- 

+ 

P 

12a 

+ 

- 

+ 

+ 

- 

NP 

12b 

+ 

- 

+ 

+ 

- 

NP 

12c 

+ 

- 

+ 

+ 

+ 

+ 

P 

13 

+ 

+ 

- 

- 

NP 

14 

+ 

+ 

- 

+ 

NP 

15 

+ 

+ 

+ 

- 

NP 

16a 

+ 

+ 

+ 

+ 

- 

NP 

16b 

+ 

+ 

+ 

+ 

- 

NP 

16c 

+ 

+ 

+ 

+ 

+ 

+ 

P 

where  a  "+"  (resp.  a  means  that  the  corresponding  parameter  is 
posi ti ve  (resp.  negati ve) ,  a  blank  means  that  the  corresponding  parameter 
value  is  irrelevant,  and,  P  (resp.  NP)  stands  for  persistence  (resp. 
non-persi stence)  of  system  (4.1). 


Proof:  We  will  illustrate  the  proof  by  proving  the  non-persistence  case 
(la)  and  the  persistence  case  (lb).  Under  case  (1),  E^  and  do  not 
exist,  En  (resp.  E0)  is  globally  stable  on  H  (resp.  H  )  ,  E. 

I  c  x-j  5 y  5 j  d 
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is  unique  and  is  globally  stable  on  H  (case  (iii)  of  Theorem 

X-|  ,  A£  , 

3.9).  In  case  (la),  d  <  0  .  Therefore  is  asymptotical ly  stable 

4 

on  R+  and  hence  system  (4.1)  is  non-persistent.  In  case  (lb),  it  is 
easily  seen  that  E^.E^Eg  i  w(xlo »x2o ,x3o ,yo^  for  a11 

^xlo,x2o’x3o,yo'  €  R+  *  (Thls  type  °f  argument  has  been  used  in  the 
proof  of  Corollary  4.9  and  will  not  be  reproduced  here.)  Therefore  the 
global  dynamics  on  b(R^)  implies  that  w  n  b(R^)  =  <J>  as  desired. 

q.e.d. 


4.6  Discussion. 

In  this  chapter  we  have  considered  a  predator-prey  model  where 
the  prey  population  consists  of  three  genotypes  with  fertility  /  viability 
differences,  and  with  different  predator  functional  responses.  We  have 
given  conditions  for  persistence  and  non-persistence  of  the  prey  geno¬ 
types. 

The  results  in  Section  4  can  be  interpreted  as:  if  a  prey 
gamete  type  has  selective  disadvantages  as  well  as  predation  disadvantages, 
then  it  will  die  out  eventually.  As  the  results  in  Section  5  show,  the 
ordering  of  the  selective  values  and  of  the  predation  functional  responses 
are  necessary  for  the  results  in  Section  4  to  hold. 

As  part  of  the  conditions  for  persistence  in  Theorem  4.11,  we 
have  shown  that  in  some  cases  where  the  subsystem  modelling  the  prey  popu¬ 
lation  growth  exhibits  non-persistence,  the  total  system  exhibits  persis¬ 
tence.  This  may  be  viewed  as  predator  regulated  survival  among  prey 
genotypes,  which  otherwise,  due  to  fertility /  viability  differences, 
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would  become  extinct. 

In  Freedman-Wal tman  (1978,  1980),  conditions  were  given  under 
the  assumption  of  no  ferti 1 i ty  /  viabi 1 i ty  differences,  for  two  of  the 
genotypes  to  become  extinct.  Here  we  have  generalized  these  conditions 
to  include  the  case  of  fertility  differences.  We  have  also  given  condi¬ 
tions  changing  the  outcome  to  one  of  persistence,  this  time  due  to 
fertility  /  viability  differences. 


CHAPTER  5 


A  TWO-LOCUS,  TWO-ALLELE  MODEL 


In  this  chapter,  the  model  considered  in  Chapter  3  for  one 
locus,  two  allele  will  be  extended  to  two  loci.  The  fertility  /  viability 
of  the  genotypes  will  be  taken  to  be  equal.  It  is  shown  that  solutions 
converge  to  rest  points  on  the  Hardy-Weinberg  manifold. 


5.1  Introduction  and  the  Model. 

Consider  a  two-locus,  two-allele  problem  where  A,  a  denote 
the  two  allele  types  at  the  first  locus  and  B  ,  b  denote  the  two  allele 
types  at  the  second  locus.  There  are  four  gamete  types  and  ten  geno¬ 
types.  If  the  number  of  gametes  of  type  AB  ,  Ab  ,  aB  and  ab 
are  denoted  by  u-|  ,  u^  ,  u^  and  u^  respectively,  then  it  is  natural  to 
denote  the  number  of  organisms  of  genotype  AB/AB  ,  AB/Ab  ,  AB/aB  ,  AB/ab  , 
Ab/Ab  ,  Ab/aB  ,  Ab/ab  ,  aB/aB  ,  aB/ab  and  ab/ab  by  x^  ,  >  x-j  3  »  , 

x22  *  x23  *  x24  *  x33  *  x34  »  x44  nespecti vely.  In  the  two-locus,  two- 
allele  model  proposed  here,  the  genotypes  AB/ab  and  Ab/aB  will  not 

be  distinguished.  We  use  x^  to  denote  the  number  of 
organisms  of  the  genotype  AB/AB  (another  notation:  AABB)  and  use 

x12  »  x13  5  X21  ’  x22  5  x23  *  X31  5  x32  and  X33  resPecti vel^  for  that  of 
the  genotypes  AABb  ,  AAbb  ,  AaBB  ,  AaBb  (AB/a  b  and  Ab/aB)  ,  Aabb  ,  aaBB  , 

aaBb  and  aabb  . 


The  model  we  wish  to  consider  in  this  chapter  is  of  the  form: 
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where 


U1  X 

X11  =  — 2  B(x)  -  -11  D(x) 

X 


•  2u, x,  9 

x12  =  B(x)  -  -jL  D(x) 


•  ^9  X-i  o 

x13  =  -f  B(x)  -  -ii  D(x) 

X 


2u-,u„ 

:21  =  — p-  B(x)  - 
x 


'21 


D(x) 


x. 


.  2(u,u,.+u9uJ  ~99 

**  ~  - L£2_3_  B(x)  .  _22  D(x) 


22 


x 


(5.1) 


2u?u  x 

x23  =  B(x)  -  -f-  D(x) 

x 


U-5  X--, 

x-n  =  B(x)  -  D(x) 


'31 


x 


x 


•  ^9^4  X  oo 

x32 =  ~^~r  B(x)  •  -jr  D(x) 


X 


33 


A 

T 

x 


x 


33 


B ( x )  -  D(x) 


x 


U1  X11  +  2  x12  +  2  x21  +  4  x22 


U2  X13  +  2  X12  +  2  X23  +  4  X22 


u3  X31  +  2  X21  +  2  X32  +  4  X22 


(5.2) 
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and 


U4  X33  +  2  X23  +  2  X32  +  4  X22 


x  =  Xn  +  •**  +  x33  . 


(5.3) 


For  a  derivation  of  system  (5.1),  see  Appendix  2.  The  birth  and  death 
functions  (B  and  D)  are  assumed  to  satisfy  the  same  hypotheses  as  those 
for  system  (3.1),  namely,  (HI )  -  (H6)  of  Section  3.1.  System  (5.1)  is  a 
no  fertility  /  viability  difference  generalization  of  system  (3.1)  to  two 
loci.  It  can  also  be  considered  as  a  continuous  version  of  the  well- 
known  model  of  two  loci  for  non-overlapping  generations  (see  Chapter  8 
of  Roughgarden  (1979)).  For  the  discrete  model  of  two  loci  (with  random 
mating,  no  selection,  linked  or  unlinked  loci),  it  is  known  that  the 
Hardy-Weinberg  equilibrium  relation  is  attained  gradually  over  a  number 
of  generations  (see,  for  example.  Section  2.6  of  Crow-Kimura  (1970)). 

An  analogous  result  is  true  for  system  (5.1).  In  contrast  to  this,  the 
recent  work  of  Akin,  Akin  (1979,  1982,  1983),  shows  that  in  the  Kimura 
model  for  two  loci  with  selection  and  recombination  (see  p.  197  of  Crow- 
Kimura,  1970),  it  is  possible  for  stable  periodic  solutions  to  exist. 

Other  related  work  of  interest  are:  Nagylaki-Crow  (1974),  Karlin  (1975), 
Karl  in-Carmel  1 i  (1975a,  1975b),  and  Nagylaki  (1977). 

The  rest  of  this  chapter  is  organized  as  follows.  In  Section 

Q 

5.2,  the  semi  flow  on  cl(£+)  defined  by  system  (5.1)  is  described  and 
some  of  its  elementary  invariance  properties  are  considered.  The  rest 
points  for  system  (5.1)  are  described  in  Section  5.3.  It  is  shown  that 
the  rest  point  set  (excluding  the  origin),  forms  a  2-dimensional  surface 
(with  boundary)  and  has  a  global  parametrization.  In  Section  5.4,  we 
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prove  the  main  result  of  this  chapter:  every  solution  of  system  (5.1) 
converges  to  a  rest  point  as  t  tends  to  +°°  .  It  turns  out  that  which 
rest  point  a  solution  converge  to  is  predictable  and  is  dependent  only 
on  the  initial  condition.  Finally,  Section  5.5  is  a  short  discussion  of 
the  results. 


5.2  The  Associated  Flow. 

As  with  system  (3.1),  system  (5.1)  is  not  defined  when 
X  =  (x-ji  ,•••  ,Xg3)  =  (0, •••,())  .  Let  F  =  (F-j-^  ,•••  9^33)  denote  the  vector 
field  defined  by  system  (5.1)  (i.e.  the  right  hand  side  of  (5.1)).  Define 
F  =  ((),•••, 0)  when  X  =  (0, •••,())  .  Then  system  (5.1)  becomes 

X  =  F(X)  .  (5.4) 

Clearly  F  :  cl(R^)  ->  R9  is  locally  Lipschitz.  As  in  (3.5),  we  define 
the  simplex 


U={XeR^:x=K} 


and  let 

cl (U)  =  (X  e  cl(R^)  :  x  =  K)  (5.5) 

be  the  closure  of  U  in  R9  .  We  will  summarize  some  of  the  properties 
of  system  (5.4)  in  the  following  proposition. 

Proposition  5.1:  System  (5.4)  defines  a  continuous  semi-flow  on  cl(R+) 
which  becomes  a  smooth  ( C  )  local  flow  when  restricted  to  R+  .  The 
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closed  simplex  cl(U)  defined  in  (5.5)  is  positively  invariant  and  is 

globally  stable  on  cl(R^)  \  {(0, ••*,())}  .  (0,***,0)  is  a  repelling 

rest  point.  H  ,  H  ,  H  ,  H  ,  H  w  ,, 

X11  X1  3  X31  X33  X1 1  *X1 2 ,X1 3  5  Hx-,rx21,x31  ’ 


,  and,  H  are  invariant.  The  vector  field  F 

*  - 5  V  Y  V  -  - 

31’  32’x33 


X13’X23’X33  *  *  x01,x00,x 


points  into  R+  on  b(R^)  \  cl (H  U  H  U 

X11  ’  12’  13  11 ’21 ’31 


H  U  H  )  . 

X13’X23’X33  X31’X32’X33 


Q 

Proof:  The  positive  invariance  of  cl(R+)  can  be  easily  shown  using 
Proposition  2.1.  As  with  system  (3.1), 

x  =  B(x)  -  D(x)  .  (5.6) 


Hence,  for  all  x(0)  >  0  , 


x(t)  K  (5.7) 

as  t  +°°  .  The  positive  invariance  and  global  stability  properties 
of  cl (U)  follow  from  (H5),  (H6)  of  Section  3.1  and  (5.7).  The  rest 
of  the  assertions  can  be  proved  as  in  Proposition  3.1. 


g.e.d. 
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5.3  The  Rest  Point  Set. 

In  this  section  the  set  of  rest  points  for  system  (5.1)  will 
be  described.  By  Proposition  5.1,  all  rest  points  of  (5.1)  other  than 
(0, •••,())  must  lie  in  cl(U)  .  First,  we  consider  the  interior  rest 
poi nts . 

Proposition  5.2:  Define  H  :  R2  ->  U  by 

2 

H(c-|,c2)  =  (H-|1(c1,c2),‘-*,H33(c1,c2))  for  (c-j,^)  e  R+ 


where 


2  2 

^11  ^C1 ,C?)  ~  ^1  ^"2 


11 v  1 *2 


(1+c-j  )2(1+c2)2 


K 


(l+c-j  )2(1+c2)2 


K 


K 


(1+c-,  )2(1+c2)2 


H00(c,  ,c9)  =  4c-, c 


K 


22 v  1 *2 


(1+c-,  )2(1+c2)2 


Hoo ( Ci » Co )  2c 


K 


23 v  1 *2 


(1+c-,  )2(1+c2)2 


K 


(1+c-,  )2(1+c2)2 
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K 


Hoo(^'"i  *Cq)  9  9 

32  1  2  2  (1+c1)2(1+c2)2 


H33*c1,c2^ 


K 


(1+c1)2(1+c2)2 


(5.8) 


2  9 

Then  H(R+)  is  the  set  of  all  rest  points  of  system  (5.1)  j_n  R  . 


Proof:  First  we  show  that  all  points  of  the  form  H(c-j,c2)  where 

9  ★ 

(c-j,c2)  e  are  rest  points  of  system  (5.1).  Let  H(c-j,c2)  =  X  = 

*  *  *  q 

(x^  ,•••  ,x23)  •  Clearly  X  e  R+  .  Define 


U-,  = 


1  *11 


*  1  * 
X,,  +  fr  X 


1  *  1  * 

+  —  V  +  —  X 

2  "12  2  21  4  22 


Un  = 


IU  = 


★  1  *  1  *  1  * 

V  4-  _  X  4-  —  X  +  —  X 

13  2  12  2  23  4  22 


*  1  *  1  * 
X~  n  +7Xnl  +  7T  X 


'31  2  21 


1  * 

2  "32  +  4  X22 


(5.9) 


and 


Then 


*  *  ]  *  i  *  ]  * 

u4  x33  +  2  X23  +  2  X32  +  4  X22 


x  =  X11  +  •**  +  X33 


U1  = 


c1c2K 


(1+C-j  )  (1+C2) 


*  c]  ^ 
u„  = 


(1+C-J  )  (l+c2) 


u3  = 


c2K 


(1+c-j  )  (l+c2) 


U4 


K 


(1+c-j )  ( 1  +c2 ) 


and  x  =  K  .  Therefore 
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*2 


Fn'x*> 


.  o 


*2 


* 

B(x  ) 


4 «„*,  ■  ^  - ,> i®  ■  o 

X 


Similarly,  one  shows  that 


F12(X*)  =  •••  =  F„(X*)  =  0 


33 


and  hence  X  is  a  rest  point. 


Next,  we  show  that  every  interior  rest  point  of  system  (5.1) 

2  * 

is  of  the  form  H(c-|  ,0^)  for  some  (c-|,c2)  e  R+  .  Let  X  = 

★  ★ 
x^3)  be  a  rest  point  of  system  (5.1)  with  x^.j  >0  ( i , j  =  1,2,3)  . 

★  ★  *  *  *  * 

Define  u,  ,  u0  ,  u„  ,  u„  ,  and  x  as  in  (5.9).  Clearly  x  =  K  .  Also 


'1  ’  “2  ’  3  ’  4 


define 


and 


Then 


d .  . 

TJ 


* 

x.  . 
ij 


(i  »J  =  1,2,3) 


x 


33 


r  •  = 


* 

u . 

l 

★ 

x 


(i  =  1 , *  *  *  ,4) 


33 


(5.10) 


(5.11) 


II 

CO 
*  CO 

dn  K 

(i) 

* 

33  = 

d12  K 

(ii) 

* 

33  = 

d13K 

(iii) 

* 

33  = 

d21K 

(iv) 

* 

33  = 

d22K 

(v) 
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Define  >  0  by 


* 

33  =  d23K 

(vi) 

★ 

33  =  d31K 

(vii) 

★ 

33  =  d32K 

(viii) 

33  =  d33K  =  K 

(ix) 

x. 


C-,  = 


 23 


2x 


33 


and 


x. 


Co  = 


 32 


2x 


33 


Then  2c-j  =  d^  and  2c2  =  d^  • 


(ii)  .  2r2  .  d12 


ITT  •  ~ 


and 


1 


'll 


d12  4dlld13  ’  (vii )' 


9lve  W 

( i  v 

9"1  ve  -Vn 


2rl  .  d21 


r3  d31 

and  (n)  .  Hi  =  ^23  =  d 

and  Till  •  r4  d33  d23 


/  i 1 ).  .  . 
TiTiJ  • 

2ri 

r2 

.  dl  2 
d13 

give 

(ii).  (ii)  . 

ttt  imy  • 

.  2r4  _ 

d32 

and 

(viii) 

.  2r3  _  d32  _ 

r3 

d31 

(ix) 

r4  d33 

d32 

d  =  c2 

(iv)  . 

2r3  .  d 

21  and  (M. 

31  2 

• 

(i)  ‘ 

r-j  d 

n  (vii) 

(iv)  . 
(vii)  ‘ 

d21 

■  2C2 

/dll  ' 

(vi)  .  2r4 
(i i i )  r2 

_  d23 
dl  3 

»<»  a.s-‘? • 


Thus , 


di 2  2ci y/d l7 


Since 


r2  .  C1 


vd 


11 


13  =  _ 
rl  M 


4  _  1 


11 


v/d 


11 


and , 
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(v)  .  2r4  j.  2r2r3  .  d22 

(i)  '  r,  r2  '  dn 


therefore 


d22  "  2c1c2  +  2  y/dll 


Substituting  all  of  the  above  into  the  equation 


* 

x 


K  ,  we  have. 


(dH  +  2ci'/VT  +  c2  +  2c2/d^  +  2cic2  +  2^y\  +  2c  +  c2  +  2c2+1J  x*3  =  K 

which  in  turn  implies  that 

*  _  1 

X33  ZZ  J 

(y/dll  +  C1  +  c2  +  ^ 

Substituting  this  into  (i)  and  by  noting  r-j  =  +  j  +\  d2i  +\  d22  » 

2  2 

we  have,  d^  =  c^c2  .  From  this,  one  can  easily  show  that 
X*  =  H ( c-j  ,c2)  . 

q.e.d. 


Next,  we  observe  that  the  boundary  rest  points  of  system  (5.1) 

must  lie  on  cl(H  )  ,  cl(H  x  )  ,  cl (H  )  or 

11*  12*  13  11’  21’  31  31  *  32 *  33 

cl(H  )  .  From  this  we  have  the  following  proposition. 

x13*  23*  33 


2 

Proposition  5.3:  cl(H(R  ))  is  the  set  of  rest  points  of  system  (5.1) 


in  cl(R;)  \  { ( 0 , • • • ,0)}  . 


Proof:  This  follows  from  Proposition  5.2,  the  above  observation  and  case 

(i)  of  Proposition  3.3.  Note  that,  for  example,  cl(H(R+))  n 

cl(H  )  is  the  Hardy-Weinberg  parabolic  segment  given  in 

X1 1 ,x21  ,x31 
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Proposition  3.3. 


q.e.d. 


Remark:  As  in  Section  3.3,  cl(H(R^))  =  H(  [0,+o°K0,+°°] )  is  called 

the  Hard.y-Weinberg  manifold.  In  this  case  it  is  a  2-dimensional  mani¬ 
fold  (surface)  with  boundary. 


5.4  Global  Convergence  to  the  Hardy-Weinberg  Manifold. 

In  this  section  we  will  show  that  every  solution  X ( t )  of 
system  (5.1)  with  X (0 )  e  cl(R^)  \  {(0, •••,())}  converges  to  a  rest 
point  in  cl(H(R^))  as  t  •+  +°°. 

First  we  define  some  auxiliary  quantities.  Define 


r  XAA  =  X11  +  X1 2  +  X1 3 

<  xAa  =  X21  +  X22  +  x23 

xaa  =  X31  +  x32  +  x33 


(5.12) 


'  XBB  =  X11  +  X21  +  X31 

<  xBb  =  X12  +  X22  +  x32 

^  xbb  =  X13  +  X23  +  x33  5 

UA  XAA  +  \  xAa  ua  =  xaa  +  2  XAa 


(5.13) 


(5.14) 


and, 
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UB  XBB  +  2  XBb 


Ub  xbb  +  2  xBb 


(5.15) 


Clearly 


and 


UA  +  ua  =  XAA  +  xAa  +  xaa  =  x 


UB  +  ub  =  XBB  +  xBb  +  xbb  =  x  > 


UA  =  U1  +  u2  ua  =  u3  +  u4 


UB  =  U1  +  u3  ub  =  U2  +  U4  ' 


(5.16) 


(5.17) 


UA(0) 


uB(°) 


Proposition  5.4:  Let  c-|  =  -  ^  and  =  u  ^ 

a  b 


Then 


(i)  uA(t)  =  c1ua(t)  ,  uB(t)  =  c2ub(t) 


(5.18) 


for  all  t  >  0  , 


(ii)  (uA(t) ,ug(t)) 
and 


(uB(t) ,ub(t)) 


-> 


C1K  K 


) 


i+ci  ’  i+^/ 


K 


c2K 

T+^  ’  w2> 


(5.19) 


(5.20) 


as  t  +°°  ,  and. 


( i i i )  (xAA(t),xAa(t),xaa(t)) 


-y 


aa 


(Xbb^)  ,xBb^ ,xbb^  ^  ^ 


c^K 

2c,K 

K  ^ 

(5.21) 

(l+c1 )2 

(l+c,)2  ’ 

(1+cp2/ 

c2K 

2c2K 

K  ^ 

(5.22) 

(l+c2)2 

(l+c2)2  ’ 

(l+c2)2/ 

as  t  ■>  +°°  . 
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Proof:  Clearly 


and 


-  _  _A 

(AA  T 
x 


x 


B ( x )  -  D(x) 


Aa 


2unu  xn 

— g-5-  B(x)  -  D ( x ) 


aa 


u  x 

a  B(x)  -  D ( x )  , 


T 

X 


X 


^RR 

Xnn  =  4  B(x)  -  -2H.  D(X) 


"BB  2 

X 


x 


xBb  =  4  B<x>  -  *-f-  D(x) 


xbb  =  ^  B(x)  -  ^  D(x)  . 


X 


X 


(5.23) 


(5.24) 


Using  a  result  in  Freedman-Wal tman  (1978)  (mentioned  in  the  proof  of 
case  (i)  of  Theorem  3.9),  we  know  that  (i)  holds  and  that 

xAA(t)  -  clxaa(t)  -  0  •  xAa(t)  '  2clxaa(t)  -  0 

(5.25) 

xBB(t)  “  c2xbb(t)  ^  0  »  xBb^^  '  2c2Xbb^  *  0  » 

as  t  +«>  .  (ii)  and  (iii)  follow  from  (5.7),  (5.16)  and  (5.25). 

q.e.d. 
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Remark:  Proposition  5.4  and  (5.17)  show  that 


^(t)  =  cA(t)  -  ua(t)  +  u4(t)  =  Clua(t)  -  ub(t)  +  u4(t) 

u2(t)  =  ub(t)  -  u4(t)  (5.26) 

u3(t)  =  ua(t)  -  u4(t) 

for  all  t  _>  0  ,  where 

UA(0)  u  (0) 

ci  =  yoy  and  c2 =  u^oy  • 

uA(0)  u  (0) 

Proposition  5.5.  Let  c-j  =  and  c2  =  ^  ^  .  Then 

a  b 

(u1 (t)  ,u2(t) ,u3(t) ,u4(t))  + 

(  C1C2K  c1K  c2K  K  \ 

((1+Cl){l+C2)  ’  (1+C-,  )(l+c2)  ’  (l+Cl)(l+c2)  ’  (1+c, )  (l+c2))  (5-27) 

as  t  •+  +  00 . 


Proof: 
Step  1 : 


The  proof  is  divided  into  three  steps. 

u-j  ,  ,  u3  ,  u4  satisfy  the  following  system  of  differential 

equations. 


U1  (U1  +  u-!u2  +  u-|u3  +  2  uiu4  +  2  U2U3^  V  U1  x 


•  2  1  1  b(x) 

u2  -  ( u2  +  u-| U2  +  u2u4  +  2  u-] u4  +  2”  U2U3^  2  ~  ^ 

X 


D(x) 
2  x 
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“3  =  (U3  +  U1U3  +  U3U4+IU1U4+IU2U3)  ^  ' 

X 

u4  -  (u2  +  u2u4  +  u3li4+lUlu4  +  lu2u3) 

X 

where  x  =  u]  + u2 + u3 + u4  .  This  follows  directly  from  (5.1),  (5.2), 
and,  (5.3). 


DUi 
3  x 


.  D(x)_ 
4  x 


(5.28) 


Step  2: 


u1(t)-c2u2(t)  =  (ulo-c-,u2o)  exp 


ft  /(l+c1)ua(T)B(x(x))  d(x(t))- 


LJ0 


2x2(t) 


x(t) 


di 


u3(t)-c2u4(t)  =  ( u3o"c2u4o^  exP 


t /r(l+c1)ua(x)B(x(T))  d(x(t))Nj 

L'oV  ^ 


u1(t)-c]u3(t)  =  (ulo-Clu3o)  exp 


ft  /(l+C-,)Ua(T)B(x(T))  p(x(T)) 


LJ0  \  2x^(t) 


x(t) 


dx 


u2(t)-c1u4(t) 


(u2o-clu4o)  exp 


t  /(l+c1 )ua(x)B(x(x)) 


LJ0 


2x2(t) 


(5.29) 


This  can  be  proved  by  writing  down  the  linear  differential  equation 
that  each  of  these  functions  satisfies.  For  example,  by  (5.26)  and 
(5.28), 


U2~C1U4  ^  U2~C1 U4  ^ 


(l+c2)ubB(x) 


2x 


D(x) 

x 


(U2“C1 U^) 


r  (l+c-,)uaB(x) 


D(x) 

X 
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Step  3: 

u]  (t)  -  c2u2(t)  0  u3(t)  -  c2u4(t)  +  0 

u-j  (t)  -  u3( t)  0  u2(t)  -  c1u4(t)  ->  0  ,  (5.30) 

as  t  +  4co  .  By  (5.7),  (5.19),  and,  (5.20), 

(HCl)ua(t)B(x(t))  D(x(t))  B(K) 

2x2(t)  T  '  2K 

as  t  -►  +°°  .  Therefore,  the  integral  in  (5.29)  is  divergent  and  conse¬ 
quently  we  have  (5.30) . 

Using  (5.30),  it  is  then  easy  to  see  that  (5.27)  holds  because 

u-j  +  u2  +  u3  +  u^  =  x  .  (5.31 ) 

q.e.d. 

u«(0)  u  (0) 

Theorem  5.6:  Let  cn  =  — anc)  co  =  — nrv  •  Then 
-  1  ua(0)  2  ub(0)  - 

(x-ji  (t),**«,x33(t))^  H(c.j,c2)  (5.32) 

as  t  ->  +°°  ,  where  H(c-|,c2)  is  defined  in  (5.8). 

Proof:  To  show  (5.32),  it  suffices  to  show 

x 1 1  ( t )  -  c-| c2x33( t)  ,  x-j  2 ( t )  -  2Ci  c2x33(  t)  ,  x -|  3 ( t )  -  Ci ^23^ "^)  5 

O 

x21(t)  -  2c1c2x33(t)  ,  x22(t)  -  4c1c2x33(t)  ,  x23(t)  -  2c]x33(t)  , 

X31  (t)  -  c2x33(t)  ,  x32(t)  ~  2c2x33( t)  0 


(5.33) 
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as  t  -*  +°°  .  We  illustrate  this  by  showing  z(t)  =  "^)  ‘  ^C2X33^^  ^  ® 

as  t  -*  +»  .  Clearly  z(t)  satisfies  the  linear  differential  equation 


z(t)  +  Q-j  (t)z(t)  =  Q2( t) 


where  Q-j  (t)  =  ,  Q2(t)  =  2u4(t)  [u3( t)  -  c2u4( t)] 

Solving  this  equation  yields 


z(t)  =  z(0)  exp  [- 


0 


Q-|  (x)dT]  + 


-t  ft 

exp  [- 

0 


Q-j  (x)dT]  Q2(s)  ds  . 


We  will  show  that  the  first  and  second  terms  tend  to  0  as  t  tends  to 
+  00  . 


Since  lim  Q,  (t)  =  ,  by  (3.7),  therefore 

t++co  K 


-foo 


Q,  (t)dt  =  +°°  ,  and  hence,  lim  exp[- 

0  t->+°° 


0 


Q-j(T)dT]  =  0  . 


On  the  other  hand,  by  (5.29), 


0 


exp  [- 


Ql (x)dT]  Q2(s)ds  =  2(u3(0)-c2u4(0) )  exp[- 


[r  DWl)i  dT]  . 

0 


xT^T 


0 


1+c, 

u4(s)  exp  [— y 


S  dT]  BWlU  ds  .  (5.34) 

0  x  ^t)  x  (t) 


To  show  that  the  L.H.S.  of  (3.34)  tends  to  0  as  t  tends  to  -*»  , 
it  suffices,  by  L'Hospital's  rule,  to  show  that 
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u4(t)  exp[- 


1  im 

t->+°o 


1+c2  f*  ^4(^)B(xCx))  B(x(t)) 

— - tph - a  J  — ? — 

2  Jo  X{T)  x2(t) 


exp[ 


D(x(t) )  riTl  D(x(t)) 

~7(7rdT]  ~w 


=  0 


(5.35) 


0 


Since 


1  im 

t-*+°° 


u4(t)  =  (l+Cl)(l'+c2) 


5 


1  im 


B ( x ( t) )  .  B(K) 
x  ( t )  K 


>  0  ,  and. 


1  im 

t-v+°o 


D(x(t))  _  D(K) 
x  ( t )  K 


>  0 


9 


therefore,  it  suffices  to  show, 


lim 

t  -H-°° 


0 


■D(x(x 
'  x(t 


(1+c2)u4(t)B(x(t)) 

x2(t) 


•]  dx  =  +00 


But  this  follows  from 


lim 

t-*-+°° 


D  ( x  ( t ) ) 

x(t) 


(l+c2)u4(t)B(x(t))i  Q{K)  l+c2  K  B(K).B(K)^n 

7 M  _1  =  K  '  2  ^+C2  *  K2  '  2K 


q.e.d. 


Remark:  Theorem  5.6  implies  that  cl(R+)  \  ,0)}  is  foliated  by 

the  strong  stable  manifolds  Wss(H(c-|  ,c2) )  »  ((c-|,c2)  e  [0,+°°] x  [0,+°°] ) 

2 

of  the  rest  points  H( c-j  ,c2)  on  the  Hardy-Weinberg  manifold,  cl(H(R+))  . 
For  0  <  c-j  ,  c2  <  +00  » 


WSS  ( H  ( c-j  ,c2 ) )  n  cl(R+)  =  ' {(x^  ,*«*,x33)  6  cl(R^)  \  { (0,  •  •  •  ,0) }  : 
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1-c 


(X1 l+xl 2+xl 3)  +  2  1  ^X21+X22+X23^  "  C^X31+X32+X33^  "  0 


1  'Co 

(X1 1+X21+X31 )  +  ~~2  ^X12+X22+X32^  "  C^xl 3+X23+X33^  =  °* 


is  7-dimensional.  For  c-j  =  0  ,  0  <  <  +00  , 


WSS(H(0,c2))  n  cl(R^)  =  {(xir---,x33)  e  cl(Hx  >)(  >){  )  \  {(0,- •  •  ,0)} 

31  32  33 


1  “c2 

X31  +  2  X32  C2X33  " 


is  2-dimensional.  For  0  <  c-j  <  +°°  ,  =  0  , 


WSS(H(c,,0))  n  cl(R®)  =  {(x11(”*,x„)  ecl(H 

I  +  u  33  X1  3 9  23 9  33 


)  \  {(0.---.0)} 


]-cl 

X1 3  +  2  X23  "  C1X33  0} 


is  2-dimensional.  For  c-j  =  +°°  ,  0  <  <  +°°  , 


Wss(H(+co,c2))  n  cl(R+)  =  { (x-j -j  »,**>x33)  £  cl(Hx  ^x  ^ 

11  12  13 


) \{ (0» • • • ,0) } 


1  “c2 

X11  +  "T~  x12  '  C2x13  0} 


is  2-dimensional.  For  0  <  c-j  <  +°°  ,  c2  =  +°°  , 


w  S(H(cr+oo))  n  cl(R+)  =  {(x11,...,x33)  £  cl  (Hx  jX  jX  )\{(0, —  ,0)> 

1 1  21  31 


1_ci 

X11  +  “2  X21  ‘  C1X31  "  0} 
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is  2-dimensional.  Finally, 


WSS (H(0,0) )  n 

Cl (R®) 

=  H 

Y 

33 

WSS  (H(05+°°) )  n 

cl(R®) 

=  H 

X31 

WSS  (H(-+o°,0) )  n 

cl(R®) 

=  H 

X1 3 

WSS  (H(+°°,+00) )  n 

cl(R®) 

=  H 

X11 

are  1 -dimensional . 


5.5  Discussion. 

In  this  chapter  we  have  shown  that  if  (x-j -j  (t) ,  ••  •  ,x33(t) ) 
a  solution  of  system  (5.1)  other  than  (0,***,0)  ,  then 

Xy,  (t) :x12(t) :x]3(t) :x21 (t):x22(t):x23(t):x31 (t):x32(t) :x33(t) 

2  2  2  2  2  2 
+  c-j c2  :  2c^c2  :  c-j  :  2c-^ c^  :  4c-|C2  :  2c^  :  c^  :  2c^  :  1 

at  t  -*  + 00  ,  where 

x  1 1  ( 0 )  +  x-j  2(0)  +  x^  3(0)  +y  x21(0)  +y  x22(0)  +  j  ><23(0) 

Q  -  -  "I  =j  1 

X31  (0)  +  x3-|  (0)  +  x33(0)  +  2  x2]  (0)  +  2  x22^0^  +  2  x23^ 


x-j  1  (0)  +  x2i  (0)  +  x3-|  (0)  +],  x-j 2(0)  x22 ( 0 )  +ir  x32(0) 

Xi 3(0)  +X23(0)  +  x33(0)  x^2(0)  x22 ( 0 )  +\  x32^0^ 


is 


and 
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This  is  directly  analogous  to  the  discrete  case  as  discussed  in  Section 
2.6  of  Crow-Kimura  (1970).  It  should  be  mentioned  that  model  (5.1)  is  a 
first  attempt  to  extend  the  one-locus,  two-allele  model  (3.1)  to  two 
loci.  And  as  such,  many  of  the  standard  two  loci  features  (for  example, 
selection,  recombination,  etc.)  are  not  incorporated.  In  the  next  chap¬ 
ter,  the  effect  of  predator  mediated  selection  will  be  studied  using  a 
predator-prey  model  in  which  the  prey  population  is  modelled  by  system 
(5.1).  The  selection  comes  from  choosing  different  predator  functional 
responses  for  the  various  prey  genotypes. 


CHAPTER  6 


A  PREDATOR-PREY  MODEL  CONSISTING 


OF  NINE  PREY  GENOTYPES 


In  this  chapter  a  model  of  predator-prey  interaction,  for 


which  the  prey  population  consists  of  nine  genotypes  corresponding  to  a 
two-locus,  two-allele  problem  and  modelled  by  system  (5.1),  is  proposed. 
Some  sufficiency  conditions  leading  to  the  evol ution  of  pure  strains  as 
well  as  to  the  persistence  of  all  components  of  the  system  are  given. 
These  conditions  serve  to  illustrate  the  similarity  as  well  as  differ¬ 
ences  between  the  one-locus  and  the  two-locus  models. 

6.1  Introduction  and  the  Model. 

In  this  chapter,  the  model  of  the  previous  chapter  is  extended 
to  include  a  predator.  More  specifically,  the  model  to  be  studied  is 
given  by: 


x-ji  -  — ^  B(x)  -  [D(x)  +  yP ii  ( x ) 3 


x 


x 


~  ID(x)  +  yP12(x)] 


x  i  0  -  — ^  B(x)  -  ^  [D(x)  +yP-j3(x)l 
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X22 

X23 

X31 

x32 

X33 

y 


2(u1u.  +  u9u  )  x99 

— ~2~  ~~ B(x)  ■  ~jr  ID(x)  +  yp22(x)J 

x 


X?"3 

— B<x>  -  -jr  iD<x)  +  yP23(x)] 

X 


U-  x?, 

~Z  B ( x )  -  -f-  [D(x)  +  yP31(x)] 

X 


2U«U,  *oo 

— V1  BW  -  ~T  ID(x)  +  yp32(x)j 


X 


_  u4 

'  T 

x 


B  ( x )  - 


x 


f-  ED(x)  +  yP33(x): 


3  x.  . 

yl-s  +  k  l  -li  P..(x)] 
i  ,j-l  x  1J 


(6.1) 


where  u-j  ,  u^  ,  u^  ,  u^  are  given  by  (5.2),  x  is  given  by  (5.3),  and, 
s,k  >  0  .  As  in  the  case  of  systems  (4.1)  and  (5.1),  x .  •  ( i , j=l ,2,3) 

•  J 

are  the  nine  genotypes  of  the  prey  population  and  y  denotes  the  preda¬ 
tor  population.  The  birth  and  death  functions  (B  and  D)  are  assumed  to 
satisfy  the  assumptions  ( HI )  -  ( H6 )  Section  3.1.  The  predator  functional 

responses  P.  .  (i,j  =  1,2,3)  are  assumed  to  satisfy  (H7)  -  (H9)  of  Section 
3 

4.1 ,  namely. 


(H7)  P.j  :  cl(R|)  -  cI(rJ) 

(h8)  p.yo)  =  0 


(i,j  =  1,2,3)  are  C1 
(i,j  =  1,2,3) 

(i  ,j  =  1 ,2,3)  for  x  t  cl (rJ)  . 


(H9) 


P!  . ( x)  >  0 
1J 
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In  Section  6.2,  we  construct  the  flow  defined  by  (6.1)  and 
study  its  boundedness  and  invariance  properties.  Some  of  the  boundary 
rest  points  of  system  (6.1)  and  their  stability  properties  are  considered 
in  Section  6.3.  In  Section  6.4,  some  sufficiency  conditions  leading  to 
the  evolution  of  pure  strains  as  well  as  the  global  stability  of  a  boun¬ 
dary  rest  point  will  be  given.  In  Section  6.5,  some  sufficiency  condi¬ 
tions  which  guarantee  the  persistence  of  the  system  will  be  discussed. 

We  finish  with  a  short  discussion  of  the  results  in  Section  6.6. 


6.2  The  Associated  Flow. 


As  in  the  case  of  system  (4.1),  system  (6.1)  is  not  defined 
when  x  =  0  .  When  x  =  0  ,  we  define  system  (6.1)  to  be 


y  =  -sy 


(6.2) 


and  write  (6.1)  and  (6.2)  as 


xij  Fij(xll’*'*’x33’y) 
y  ~  G(x-|-| ,  *  *  *  ,x33»y)  • 


(6.3) 


We  have  the  following  analogue  of  Proposition  4.1: 

Proposition  6.1:  F .  .  (i,j  =  1,2,3)  ,  G  e  Lip  on  c1(r].^)  and 

*  J 

F..  (i,j  =  1,2,3)  ,  G  e  C1  on  r|°  .  Thus  system  (6.3)  defines  a  con- 
tinuous  semi -flow  on  cl(R^)  which  becomes  a  smooth  (C  )  local  flow 
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when  restricted  on  .  It  is  A  -  dissipative  with 


A  =  ( (x-|-|  ,•••  »x33>y)  e  cl  (r|°)  :  x  £  K  and 


kx  +  y  £  (M  +  sK)} 


(6.4) 


where  M  =  max  {B(x)-D(x):xe[0,K]}.  H  ,H  ,  H  ,  H  ,  H 

X11  X1 3  X31  x33 


H 

> 

H, 

H. 


iry  xi3,y 

X31 

,  H 

31 5  32  *  33 

xn* 

and 

31  ,X32,X33’y 

11  ’*“,F33,G) 

poi 

,  H  w  ,  H 


,  H 


,  H 


Xir*“,X33 


xn  ,x21  ,x31 »y  ’  x13,x23’x33,y  ’ 

are  invariant.  The  vector  field 


i[°  on  b(Rj°: 


V  =  cl (H 


U  H 


U  H 


X1 1  ,X1 2,X1 3,y  x-i^sX^ijX^-jsy  x^QjX^qjX^^jy 


'irfl2r 3i 


1 3  23  33 


UH  ..  U  H 

X31  ,X32,X33 


»y  xn 5  ’  *  *  ,x33 


)  . 


(6.5) 


Proof:  To  show  that  F, ]  6  Lip  on  c1(r|°)  ,  it  suffices  to  consider 
two  points  in  c1(r|^)  of  the  form  ( x-j i , •  •  •  jX^y)  and  (0#«**,0,y) 
where  x  =  x^  +  •••  +  >  0  .  Then 

<  U-J  B(x)  +  ^  £D(x)  +yP11  (x)] 

X 


Fn (xn ’*** ,x33,y^  -  Fn (°>*#* »°»y) 


<  B(x)  +  D(x)  +  yP-,1(x) 


<_  [B '  (z-j )  +  D*  (z2)  +  yP-J i  ( z3) ]  |x| 


33  5 
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for  some  0  £  z-|  ,  £  x  .  Therefore  e  Lip  on  cl(pj°)  and 

•  •  •  p 

12’  ’33 


similarly  we  show  that  ,  G  e  Lip. 


The  positive  invariance  of  c1(r|^)  can  be  easily  verified 
using  Proposition  2.1.  Hence  system  (6.3)  defines  a  continuous  local 
semi-flow  on  c1(r|^)  which  becomes  a  c"*  local  flow  when  restricted 


to  R 


10 

+ 


To  show  that  system  (6.3)  is  A  -  dissipative  for  A  defined  in 
(6.4),  it  suffices  to  observe  that 

x  £  B(x)  -  D(x) 

kx  +  y  <  k [B (x)  -  D(x)]  -  sy 


and  use  Proposition  2.2.  This  also  shows  that  system  (6.1)  defines  a 
continuous  semi-flow  on  c1(r|^)  . 

The  invariance  of  the  sets  listed  is  clear.  It  is  also  easy 
to  show  that  the  vector  field  (F, i , •  •  •  jF^jG)  points  into  on 

b(p|°)  \  V  .  For  example,  for  (X,y)  =  (x-| -| ,  •  •  •  ,x33,y)  e  b(pj°)  with 

x  i  =  0  ,  F-j  i  (X,y)  =  B(x)  >  0  and  Fr((X,y)  =  0  if  and  only  if 

X 

U-|  —  0  (i.e.  ^ii  — ^12_ ^21=  x22^  * 


q.e.d. 


no 


6.3  Some  Boundary  Rest  Points  and  Their  Stability. 

In  this  section  some  of  the  rest  points  of  system  (6.1)  which 

lie  on  b(R^)  will  be  considered.  According  to  Proposition  6.1,  these 

rest  points  must  lie  on  V  .  Clearly,  Eg(0, •••,())  is  a  rest  point  and 

by  Proposition  5.3,  E  =  { (H(c-j  ,c2)  ,0)  :  (c-j  ,c2)  e  [0,-k»]x[0,-k»]}  is 

the  set  of  rest  points  of  system  (6.1)  in  cl(H  )\  {En}.  Let 

xll’*'*’x33 


bij  =  -s  +  k  p i  j  ( K )  ( i ,  j  =  1 , 3 )  . 


(6.6) 


Then  H  (i,j  -  1,3)  contains  a  rest  point  if  and  only  if  b..>  0 

xi  j  *y  1J 

(i,j  =1,3)  .  If  they  exist  these  rest  points  are  unique.  Let  us  denote 

these  rest  points  by  E^  (i  =  1  ,2,3,4)  if  they  exist.  That  is,  let 

E1  =  (*11  »  E2  =  (0,0,x13,0,«  ••  ,0,y-l3)  , 

E3  =  (0,---,0,x31,0,0,y31)  ,  and  =  (0,* • •  ,0 ,^33,733) •  The  invariant  sets 

H  ,  H  ,  H  ,  ,  and,  H  w  „  may  or  may 

X1 1  ,X21  »X31  ,y  Xn,X12’X13’y  x1  3,x23,x33,y  X31,X32’  33,y 


not  contain  rest  points.  Some  conditions  are  known  which  guarantee  the 
existence  or  non-existence  of  these  rest  points  (see  Freedman-So-Wal tman 
(1984)).  In  this  thesis  we  will  not  use  these  conditions. 


The  rest  point  EQ  is  unstable.  Indeed  we  have  the  following 
Proposition. 

Proposition  6.2:  Eg  /  w(x-| -j  ,*  •  •  ,x33 ,y)  for  all  ( x-j  -j  ,•  •  •  ,><33 ,y) 6  . 


Proof:  This  follows  from  Proposition  2.15  because  in  a  small  enough 


neighbourhood  of  Eg  ,  we  have. 


Ill 


x  >_  B(x)  -  D(x)  -  yP(x)  >.  k^x  for  some  >  0 

y  £  yC-s  +  kP(x)]  <  -k2y  for  some  k2  >  0  , 

where  P(x)  =  max  (P.  . (x)  :  i  ,j  =  1,2,3}  ,  and  the  flow  on  cl(H  )  is 

'  j  y 

given  by  y  =  -sy  . 

q .e .d . 


Corollary  6.3:  (0,---,0,y)  /  w(x^  ,x33,y)  for  all 

,  ,  JO 

(Xi i , •  •  •  ,x33  ,y)  £  . 

Proof:  This  follows  from  Proposition  6.2  and  the  global  stability  of 
Eq  on  the  invariant  set  cl ( H  )  . 

q.e.d. 


The  stability  properties  of  each  (H(c-|,c2),0)  in 

cl ( X  )  \{En)  is  determined  by  Theorem  5.6.  To  compute  the 

xn,“',x33  u 

stability  of  (H(c^,c2),0)  in  the  y-direction  we  note  that 


d ( c-j  ,c2)  =  -s  + 


2 - 2  [cic2P11(K)  +  2cic2P12(K) 


(1+C-j)  (1+C2) 


+  c^P13(K)  +  2c1c|p2i  (K)  +  4c-,c2P22(K)+  2c-,P23(K) 


+  c3P31(K)  +  c2P32(K)  +  P33(K)] 


(6.6) 


is  the  eigenvalue  of  the  variational  matrix  of  system  (6.1)  at 
(H(c-|  ,c2)  ,0)  that  corresponds  to  eigenvectors  that  have  a  non-zero  y 
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component.  A  condition  for  the  persistence  of  y  is  given  by  the  follow¬ 
ing  proposition. 

Proposition  6.4:  If  d ( c-j  ,c2)  >  0  for  all  c-j  ,c2  e  [0,+°°]  ,  then 
Il^y(t)  >  0  for  all  (x11o>-"  .x33o.yo)  c  R™  . 

Proof:  The  proof  is  similar  to  that  of  Proposition  4.4.  Again  we  know 
that  E  is  repelling  in  the  y-di recti  on  and  therefore  in  order  for 
lim  y(t)  =  0  ,  the  omega-limit  set  w(x, -,  ,*••  ,x^  ,y  )  of 

t-H-oo  1  10  0 

c 

(x-j  lo  **  *  *  ,x33o  ,yo^  must  ’’ntersec‘t  the  stable  manifold,  W  (E)  ,  of  the 

s  10 

invariant  set  E  at  a  point  other  than  E  itself.  Now  W  (E)  n  cl (R  ) 

cd(H  )  .  The  global  dynamics  on  cl  (H  )  as  given 

X1 1 9  *  ’  *  ,x33  X11 *  *  33 

by  Theorem  5.6  implies  that  either  w  contains  Eq  or  w  is  unbounded. 

In  either  case,  we  have  a  contradiction. 

q.e.d. 

Remark.  Unlike  as  in  Lemma  4.5,  we  are  unable  to  obtain  a  necessary  and 
sufficient  condition  on  P-jj(K)  so  that  d(c-|,c2)  >  0  for  all 
Ci  ,  c2  e  [0,+°°]  .  However,  from  the  definition  of  d (c-j  ,c2)  given  in 
(4.7),  it  is  clear  that  d(c-|,c2)  >  0  for  all  c-|  ,  c2  e  [0,+°°]  if 
Pi  . (K)  >-|  for  all  i  ,j  =  1,2,3  . 

The  linear  stability  of  E.  (i  =  1,2, 3,4)  is  governed  by  the 
signs  of  the  real  parts  of  the  eigenvalues  of  the  variational  matrix 
M(Ei )  (i  =  1 ,2,3,4)  of  system  (6.1 )  evaluated  at  E^  (i  = 1  ,2,3,4)  . 

For  illustration,  we  will  present  M(E-j)  below.  To  simplify  notations. 


we  use 
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Mij,hk  to  denote  'I(E1)sMy>ij  to  denote  ^(E,) 

nK  i  j 

we  move  y  to  the  second  row. 


etc.  and 


xn 

y 

X1 2 

X13 

X21 

X22 

X23 

X31 

,  X 

ro 

X33 

X11 

? 

-ve 

? 

? 

? 

? 

? 

? 

? 

• 

? 

y 

+ve 

0 

? 

? 

? 

? 

? 

? 

? 

? 

X1 2 

0 

0 

? 

+ve 

0 

+ve 

+ve 

0 

0 

0 

X1 3 

0 

0 

0 

-ve 

0 

0 

0 

0 

0 

0 

X21 

0 

0 

0 

0 

? 

+ve 

0 

+ve 

+ve 

0 

X22 

0 

0 

0 

0 

0 

? 

+ve 

0 

+ve 

+ve 

X23 

0 

0 

0 

0 

0 

0 

-ve 

0 

0 

0 

X31 

0 

0 

0 

0 

0 

0 

0 

-ve 

0 

0 

x32 

0 

0 

0 

0 

0 

0 

0 

0 

-ve 

0 

X33 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-ve 

where  +ve  (resp.  -ve,0,?)  means  that  the  entry  is  positive  (resp.  is 
negative,  zero,  has  indeterminate  sign).  Furthermore,  the  non-zero 
entries  are  given  by: 


Mll,ll  =  B,(xll}  "  D 


(xll)  - 


yi  i  Pi i (xi i ) 


M 


n,y 


(x^)  <  o 


Mn,i2  Mii,2i  B'(xn)  "  D,(xn)  -  ynpii (xn } 

Mn  ,22  =  "  ~~  B ( x-j i )  +  B 1  ( xi i )  -  D'  (xn )  -  y-| i p-j i (xn }  • 

2xn 
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Mirij  ■ 

=  -  J  B(x,  i )  +  B '  ( X-j  -j )  -  D1  (x-j  -j )  -  y^P^x^) 

X11 

(for  ( i ,  j )  = 

(1,3)  ,  (2,3)  ,  (3,1)  ,  (3,2)  or  (3,3)) 

M  ,, 

y.n 

=  ky-j-jP-j-j  (X]i )  >  o 

M  •  . 
y.ij 

=  [xnP^(xn)  +  P  (xn)  -  P,  1  (x-,  n ) :  ,  (i.j)  t  (1,1) 

X11 

M12,12 

=  ~  CBtx^)  - D(xn) -YnP^tx-n):  =  cf  -  P12(xn)] 

X11  X11 

M12,13 

=  J-  B(x,,)  >  0 

X11 

M12,22 

=  1  B(x„)  >  0 

2x 

^Xll 

M12,23 

=  1  B(7 .)  >  0 

X11 

M1 3 ,1 3 

=  ~  [D(xn)  +  y11P13OTi1 )]  <  0 

X11 

M21  ,21 

i  _  _  _  y  ]  i  c  _ 

[B ( x -j  -j )  -  B ( x -j  -j )  y11p2l  (xH  ):|  —  ^1<  ^21  ^X1 1  ^ 

xn  xn 

M21  ,22 

=  B(x,,)  >  0 

2x 

^Xll 

M21 ,31 

=  1  B  ( x", , )  >  0 

X11 

M21 ,32 

=  1  B(x,,)  >  0 

X11 

115 


^22 ,22  -  B(x-j-|)  ~  DCx^-j)  -  yiiP22^xll^ 

X11 

=  C^ll(t  "  P22(xll^  ‘  \  B(xll)  ] 

X11 

M22 ,23  =  -  B^xl 1 ^  >  0 

xn 

^22 ,32  =  ~  B(xn )  >  0 

X11 

^22 ,33  =  -  B^XH  ^  >  0 
X11 

M23 ,23  =  ‘  ^D^xl  1 )  +  ^11P23^X11^  <  0 
X11 

M3i,3i  =-_t—  [D(xn)  +  ynp3i(xn)]  <  0 
X11 

m32,32  =  "  cD(xn)  +  yiip32(xn)]  <  0 

xn 

M33,33  =  '  =^~~  [D(xr|)  +  <  0 

X11 


Therefore  the  M(E-j)  has  five  negative  eigenvalues:  13  »  23  ’ 

^31  31  *^32  32  ’^33  33  ’  t*iree  rea"*  eigenvalues  with  indeterminate  signs: 
Mi 2  12  >^21  21  ’^22  22  ’  anc*»  S1*9n  rea^  Parts  °f  the  nemain- 


ing  two  eigenvalues  (corresponding  to  eigenvectors  lying  on  H 


x 


)  is 


11*- 


11 ,11  • 


the  same  as  that  of  the  sign  of  M 
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6.4  Some  Conditions  Which  Lead  to  Evolution  of  Pure  Strains. 


In  this  section,conditions  leading  to  the  global  stability  of 
boundary  rest  points  will  be  considered.  Due  to  the  complexity  of  the 
problem,  we  are  forced  to  focus  on  a  small  number  of  cases.  In  our  analy¬ 
sis,  we  will  need  to  know  the  global  dynamics  of  system  (6.1)  on  the 
boundary  invariant  set  V  .  The  assumptions  are  essentially  those  in 
Theorem  4.6  applied  to  each  of  the  four  4-dimensional  boundary  invariant 


sets  H 


,  H 


,  H 


X11  ,x12’X13,y  xirx21,X31,y  X13’X23,X33 


,  and  H 


X31  ,X32’X33 

The  added  assumptions  for  this  section  and  the  next  are  as  follows. 


(HT)  The  predator  y  persists.  Mathematically  speaking,  this 
means  that,  according  to  Proposition  6.4,  we  require 
d(c-,  ,c2)  >  0  for  all  c-,  ,  e  [0,+°°]  . 


(H2)  E-|  (resp.  E2,E3,E4)  is  globally  exponentially  stable 


on 


Hx  v  (resp*  Hx  y  •  H 
xll’y  x13’y 


,  H  w), 

x3i  »y  x 


33 ,v 


(Hi)  The  four  sets  of  predator  functional  responses:  P-,-,  /  P-j 2  /  p13 

Pll/P21/P31  ’  P13/P23/P33  ■  and’  P31  '  P32  '  P33  are 
ordered  with  the  heterozygote  one  between  the  two  homozygotes 

ones.  There  are  16  possibilities.  By  interchanging  the  roles 
of  the  allele  types  A  (resp.  B)  and  a  (resp.  b),  if  neces¬ 
sary,  we  are  left  with  3  cases  to  consider: 

(1)  P-ji  £  P12  <  P-,3  ,  P-,  i  1  P2i  1  p3i  »  P1 3  -  P23  -  P33  » 


P31  -  P32  -  P33  * 
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(2)  PH  -  P1 2  -  P13  9  P11  1  P21  -  P31  5  P33  -  P23  -  P13  * 


P33  -  P32  -  P31  *  and 


(3)  Pn  <  P]2  <  P]3  >  P-ji  1  P2i  1  P3i  »  Pl 3  1  P23  -  P33  9 


P33  -  P32  -  P31  * 


The  case  when  P]1  <  P]2  <  P]3  <  P23  1  P33  <  P23  <  P31  <  Pgl  1 
P-jl  is  excluded  because  we  require  at  least  one  of  the  inequalities  to  be 
strict. 


In  each  of  the  above  three  cases  (1)  -  (3),  we  require  the  dif¬ 
ferences  of  the  appropriate  P.  .  ' s  be  bounded  below  by  a  non-negative 

*  vJ 

linear  function  of  x  on  [0,K]  .  For  example  P^  £  P12  <  P13  in 
case  (1)  should  read:  there  exist  ^  »  k-j2  -j3  _>  0  ,  not  both  zero, 
such  that  P12(x)  -  P^x)  >  ]2x  and  P-, 3(x)  -  P12(x)  >  k]2  13x  for 

x  e  [0 ,K]  . 

The  following  is  a  theorem  which  provides  conditions  under 
which  the  rest  point  F-j  is  globally  stable. 


Theorem  6.5:  Suppose  (HI)  and  (H?)  hold ,  Furthermore  assume  that 
(l)  P^  <  P]2  <  P13  1  P21  1  P22  -  P23  -  P31  -  P32  -  P33  9 

(ii)  there  exist  e  -j  2  9^13  9^21  9  ^  2  2  9  ^  2  3  9^31  9  ^3 2  9  ^33  ^  ^  * 
such  that 


P1 2 ( ^ )  —  p "j  1  ( x )  _>  e-j 2x  ,  P] 3 ( x )  ~  P"j2(^)  ^  ei3x  * 

P2i  (x)  ~  p  1 3  ( x )  il  e2i x  »  P22^x)  —  P2"|(x)  e22x  , 

P23(x)  —  P 2 2 ( x )  e23x  ,  p 3 1  ( x )  —  P23(x)  £31  x  , 
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P32 ^  ~  P31  (x)  7_  e32x  *  -  P32 ( x )  —  e33x 


for  all  x  e  [0,K]  ,  and. 


(iii)  y13%  -  P22^x13^  ^  7  B^x13^ 


y31%"  P22^x31^  ^  7  B^x31  ^  »  — 


1  n  t~T7 


y33^k  '  P22^X33^  ^  2  B^X33^  * 


then  E-j  is  global  ly  stable  on  R 


10 

+ 


Proof:  (HZ)  ,  (ii),  and  (iii)  imply  that  the  rest  points  (i  =  1  ,2,3,4) 
are  hyperbolic.  In  fact  by  studying  the  eigenvalues  of  the  variational 
matrix  at  each  of  these  rest  points,  as  given  in  Section  6.3,  we  know, 
dim  WU(E-| )  =  0  ,  dim  WU(E2)  >  1  ,  dim  WU(E3)  >  1  ,  dim  WU(E4)  >  2  . 
Therefore  E-j  is  asymptotically  stable  and  WU(E^)  (i  =2,3,4)  must  have 
non-empty  intersection  with  \  c1(r|B)  .  By  Proposition  2.11, 


WS(E.)  n  c!(r|B)  must  lie  in  b ( )  .  Theorem  4.6  tells  us  that 


Is 

,s 


,10' 


Ws(E4)  n  cl(R;u)  c  cl (H 


,10' 


)  ,  WS(E  )  n  cl (Rj0)  =  cl(H  ), 

x33  ,y  x31  *x32  *  33»y 


and  W  (E«)  n  cl(R_!  )  c  cl  ( H  )  .  From  this  one  can  easily 

L  X1 3  ,X23  *x33  ,y 


show  that  E-j  is  the  only  point  in  b(R^)  that  can  lie  in 
w  =  w(x^ ,••• ,x33,y)  for  any  (x-j 1  ,• • •  ,x33 ,y)  in  r|°  .  Moreover,  if 
E-j  e  w  ,  then  the  orbit  converges  to  E-j  as  t  ->  +°°  .  Using  REDUCE  2 
(a  programming  language  for  algebraic  manipulations),  one  shows  that 
u  u. 

— - -  < - [e00x,0x00  +  other  non-negative  terms]  .  Therefore,  as 

ua  UA  ~  uAua  33  12  33 

in  the  proof  of  Theorem  4.6, 
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u  (0) 

a 


ua(t)  iT^TUT  uA(t)  exP  K 


33 


0 


( X-|  2^33^ )  ( ^ )  dx  1 


If  the  integral  diverges,  then  lim  u  (t)  =  0  so  that  w  n  b(R10)  f  <f> 

t-*+°°  a  + 

and  hence  the  orbit  converges  to  E]  .  Now  suppose  z(t)  =  (x-|2x33y)(t)  e 
lJ[0,+«>)  ,  then  lim  z(t)  =  0  and  hence  1 im  x,?(t)  =  0  or 

t->-Ko  t^‘+00 

1 im  x..(t)  =  0  (since  lim  y(t)  >  0  by  (HI)  ).  This  again  implies 
t-H<x> 

w  n  b(R^)  f  (J>  and  so  the  orbit  converges  to  E-j  as  before.  This  shows 
that  E-j  is  globally  attracting  and  it  completes  the  proof  of  the  theorem. 

q.e.d. 


Remark:  The  assumption  (iii)  in  the  above  theorem  is  technical  and  is 
necessary  in  order  for  the  rest  points  E.  (i  =  1,2, 3, 4)  to  be  hyperbol¬ 
ic.  This  assumption  is  "generical ly"  satisfied. 


6.5  Some  Persistence  Results. 

In  this  section  we  will  consider  conditions  which  lead  to  the 

persistence  of  system  (6.1).  Besides  the  conditions  (HT)  -  (Hi)  made  in 

the  previous  section,  we  will  also  assume  that  the  rest  points  E^ 

(i  =  1  ,2,3,4)  are  hyperbolic.  This  amounts  to  requiring  the  k.  .  hl's  in 

I  j  5  n  I 

(Hi)  be  positive  and  that  -  P22^ij^  "  I  B^ij^  (i»3  =  1>3)  be 

non-zero. 

Theorem  6.6:  Under  the  above  assumptions,  system  (2.1)  is  persistent  if 

(i)  y-|-j[-|  -  P22^x11^  "  7  B(xn)  >  0  holds  under  cases  (1)  and 
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(3)  in  (H3)  of  Section  6.4,  and, 

(ii)  y^C |  -  P22^i  1 )]  "  \  >  0  >  and 

y33^-T<  ‘  P22^x33^  "  1  B^x33^  >  0 

hold  under  case  (2)  of  (H3)  in  Section  6.4. 


Proof:  We  will  only  discuss  case  (1)  of  (Hi)  in  Section  6.4.  The  other 
two  cases  are  similar.  Let  (x^  ,*•• ,x.„,y)  e  and  let 


v33 


+ 


,10. 


w  =  w(x^  ,•••  ,x33,y)  .  Our  aim  is  to  show  w  n  b(R| )  =  4>  ,  or  equiva¬ 
lently  w  n  V  =  <p  ,  where  V  was  given  in  (6.5).  Recall  that 

lim  y(t)  >0  by  (HI)  .  First  we  will  show  that  Ed  i-  w  .  Suppose  not, 

t->+oc' 

that  is,  suppose  e  w  .  Since  dim  WU(E^)  _>  2  ,  therefore  WU(E^)  n 

(r|B  \  c1(r|B))  1  4  ,  and  hence,  by  Proposition  2.11,  WS(E^)  n  r|^  =  4. 

It  follows  from  the  global  dynamics  on  cl(H  )  and 

x13’  23’  33,y 

c  10 

cl(H  )  as  given  in  Theorem  6.1  that  Wb(EJ  n  cl(R! )  c 

x31’  32*  33’y 


cl  (H 


x33,y 


)  .  The  Butl er-McGehee  lemma  then  implies  that  w  contains  a 


point  in  cl(H  )  other  than  E.  which  is  a  contradiction.  Next  we 
x33 ,y  4 

show  that  E^  4  w  .  As  in  the  case  of  E^  , 

Wu(E3)  >  1  implies  WU(E3)  n  (r|°  \  c1(r|°))  1  4  and  therefore,  by 


<~  1 0 

Proposition  2.11,  W  (E3)  n  R+  =  <p  .  The  global  dynamics  on 

c  10 

cl  (H  )  and  cl  (H  )  imply  that  W  (Ej  n  cl  (R  ) 

XH  »x2i  »x3i  >y  31  *  32  ’  33  ,y 

c  cl(H  )  .  The  Butl er-McGehee  lemma  says  that  w  must  contain 

X31  ,x32’x33’y 

a  point  in  cl  (H  v  Y  ..)  other  than  Eo,  which  is  a  contradiction. 

X31  ,  32  ’  33  *y 
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Hence  E^  (  w  .  Similarly  4  w  .  One  could  repeat  the  above  argument 
for  E-j  because  >  0  by  assumption.  Hence  E..  t  w  (i  =  1 ,2,3,4) 

and  w  n  V  =  cf)  follows  from  the  global  dynamics  on  V  . 


6.6  Discussion. 

Due  to  the  complexity  of  (6.1),  our  analysis  for  this  system 
is  not  as  complete  as  the  other  models  (3.1),  (4.1)  and  (5.1),  that  we 
studied  in  the  previous  chapters.  However,  Theorem  6.6  points  to  a  way 
of  making  the  system  persistent  which  otherwise  would  not  be  by  adjust¬ 
ing  the  predator  functional  response  of  the  double  heterozygote  (AaBb)  . 
This  result  indicates  a  situation  for  which  a  two-locus  problem  cannot  be 
analyzed  as  the  "sum"  effect  of  two  one-locus  problems. 


CHAPTER  7 


CONCLUDING  DISCUSSION  &  REMARKS 


In  this  thesis  we  study  four  continuous  models  in  population 
genetics  /  ecology  both  for  the  one-locus,  two-allele  problem  and  for  the 
two-locus,  two-allele  problem. 

In  Chapter  3,  a  continuous  model  (system  (3.1))  for  a  one- 
locus,  two-allele  problem  with  non-linear  birth  (B)  and  death  (D)  proces¬ 
ses  as  well  as  ferti 1 i ty  /  viabi 1 i ty  differences  is  considered.  The  popu¬ 
lation  x  under  consideration  is  divided  into  three  subpopulations  x-j  , 
x^  ,  Xg  corresponding  to  the  three  different  genotypes  AA  ,  Aa  and  aa 
with  fertility  /  viabilities  denoted  by  f-|  ,  f^  and  fg  .  The  main 
result  for  this  model  is  Theorem  3.9  which  provides  a  complete  picture  of 
the  global  dynamics  of  system  (3.1).  This  theorem  can  be  summarized  as 
follows. 

(i)  When  f-j  =  f^  =  f^  (no  selection),  all  solutions  converge  to  some 

equilibrium.  More  precisely,  if  (x-|  (t)  ,x2(t)  ,x3(t))  is  a  solution, 

K  2 

then  it  converges  to  the  equilibrium:  - y  (c  ,2c, 1)  where  c  is 

(i+cr 

X, (0)  +  J  x2(0) 

determined  by  initial  conditions,  namely  c  =  - =j -  .  Hence, 

x3 ( ° )  +  j  x2(0) 

all  solutions  converge  to  a  polymorphism  (positive  interior  equilibrium). 
Which  equilibrium  it  converges  to  depends  on  initial  conditions.  Phrased 
in  another  way,  this  result  also  says  that  the  proportion  x-j(x)  :  x^(t)  : 
x3(t)  tends  to  c2  :  2c  :  1  as  t  tends  to  positive  infinity.  This  is 
the  analogue  of  the  Hardy-Weinberg  equilibrium  relation  for  discrete  one- 
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locus,  two-allele  models  with  random  mating  and  no  selection. 

(ii)  When  f-j  ,  f3  >  (homozygous  advantage),  the  positive 
cone  is  divided  into  two  regions  by  a  separatrix  (dividing)  surface.  Solu¬ 
tions  initiating  from  one  region  will  converge  to  the  boundary  equilibrium 
(x-jjX^jX^)  =  (K,0,0)  whereas  solutions  initiating  from  the  other  region 
will  converge  to  the  boundary  equilibrium  (0,0, K)  .  Solutions  initia¬ 
ting  on  the  separatrix  surface  converge  to  the  unique  unstable  polymer¬ 
ic  2 

phism:  - (c  ,2c,l)  .  Unlike  the  case  when  there  is  no  selection,  c 

(l+cr 

is  independent  of  initial  conditions  in  this  case.  Instead,  it  is 

determined  by  the  fertilities  f-j  ,  f^  and  f^  ,  namely, 

f  -  f 

3  2 

c  =  -r - j-  .  Hence,  under  "most"  initial  conditions,  one  of  the  gamete 

T1  ‘  2 

types  becomes  extinct.  Which  one  will  actually  become  extinct  depends 
on  initial  conditions. 


(iii)  When  f^  >  ^  ,  f^  (heterozygous  advantage),  all  (positive) 

K  2 

solutions  converge  to  the  unique  stable  polymorphism:  - 0  (c  ,2c,l)  . 

(He) 

Hence  all  three  genotypes  coexist.  As  in  the  case  of  homozygous  advan¬ 
tage,  c  is  independent  of  initial  conditions  and  is  dependent  solely 


on  fertility  differences,  namely,  c 


(iv)  When  f^  >  f  >  f3  (selection  for  allele  A  which  has  incomplete 
dominance  over  allele  a)  or  when  f-j  >  f3  =  fg  (selection  for  the  reces¬ 
sive  allele  A),  or  when  f-,  =  f2  >  f3  (selection  for  the  dominant  allele 
A),  all  solutions  converge  to  the  boundary  equilibrium  (K,0,0)  .  Hence, 
the  gamete  type  a  becomes  extinct. 
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(v)  When  f^  >  f^  >  ^1  (selection  for  allele  a  which  has  incomplete 
dominance  over  allele  A),  or  when  ^3  >  ^2  =  ^1  (selection  for  the  reces¬ 
sive  allele  a),  or  when  f^  =  f^  >  ^1  (selection  for  the  dominant  allele 
a),  all  solutions  converge  to  the  boundary  equilibrium  (0,0, K)  .  Hence, 
the  gamete  type  A  becomes  extinct. 

All  of  the  above  results  have  their  analogue  with  discrete 

models. 

In  Chapter  4,  the  model  in  Chapter  3  is  extended  to  include  a 
predator  y  (system  (4.1)).  The  prey  population  x  is  modelled  by 
(3.1).  For  the  prey  population,  there  are  two  types  of  selection  forces 
at  work.  First,  there  is  the  one  studied  in  Chapter  3,  namely,  selection 
due  to  fertil ity  /  viabil ity  differences.  The  other  type  is  selection  due 
to  differential  predator  functional  responses.  (Recall  that  the  predator 
functional  response  of  x.  (i  = 1,2,3)  was  denoted  by  P.  (i  = 1,2,3)  .) 

Of  particular  interest  are  the  following  two  questions. 

(1)  Under  what  conditions  will  one  of  the  prey  gamete  types  become 
extinct? 

(2)  When  will  all  the  prey  genotypes  coexist  with  the  predator? 

The  concept  of  coexistence  is  formulated  in  terms  of  the  definition  of 

persistence:  if  x.(0)  >  0  (i=l  ,2,3)  and  y(0)  >  0  ,  then 

lim  inf  x. (t)  >  0  (i  =  1,2,3)  and  lim  inf  y(t)  >  0  .  In  the  case  of 
^  t->+°° 

no  selection  due  to  fertil ity  /  viability  differences,  the  answer  to 
question  (1)  is  provided  by  Theorem  4.6.  This  theorem  (modulo  a  number 
of  technicalities)  can  be  summarized  as  follows. 
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If  the  predator  functional  responses  of  the  prey  genotypes 
are  ordered:  P-j  <  <  P^  (hence  x-j  can  be  thought  of  as  the  one 

which  is  the  most  difficult  to  catch)  and  the  predator  y  can  survive 
on  x-j  alone  (hence  it  will  also  be  able  to  survive  on  the  other  two 
prey  genotypes  alone),  then  gamete  type  a  becomes  extinct. 

This  result  can  be  interpretated  as  saying  that  the  predator 
exhausted  the  two  prey  genotypes  X2  and  x^  which  are  easier  to  catch 
and  lives  on  the  prey  genotype  x-j  alone.  In  the  case  when  there  is 
selection  due  to  fertility  /  viability  differences,  the  answer  to  question 
(1)  is  provided  by  Theorem  4.8.  In  summary  this  theorem  states  the  fol¬ 
lowing. 


If  the  fertilities  are  ordered:  f-|  £  f^  £  f^  and  not  all 
equal  (hence  x-j  has  fertility  advantage)  and  the  predator  functional 
responses  are  also  ordered:  P-j  £  P^  £  P^  (hence  x^  has  predation 
advantage),  then  gamete  a  becomes  extinct. 

Note  that  in  this  theorem,  it  is  not  required  that  the  predator 
y  must  survive.  In  fact  it  can  also  go  extinct  along  with  x^  and  x^  • 
This  result  is  intuitively  clear  in  view  of  the  previous  results  for  no 
selection  due  to  differential  predator  functional  responses  and  for  no 
selection  due  to  fertility  /  viability  differences. 

The  answer  to  question  (2)  in  the  case  of  no  selection  due  to 
fertil ity  /  viabil ity  differences  is  provided  by  Theorem  4.10.  In  sum¬ 
mary,  this  theorem  states  the  following. 
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We  assume  that 

(i)  the  subsystem  with  only  x -j  and  y  has  a  globally  stable 

positive  interior  equilibrium  Eg  =  ,0,0 ,y-j )  (consequently 

Pi  ( K)  >  -£-)  and  the  subsystem  with  only  Xg  and  y  has  a 
globally  stable  positive  interior  equilibrium  =  (0,0,Xg,yg) 

r 

(consequently  Pg(K)  >■£-), 

(ii)  P2(x-|)  <  |  ,  P g Cxg )  <  f  ,  and, 

(iii)  P2(K)  >  f  -/<P1(K)  -  f)(P3(K)  -  f) 

then  system  (4.1)  is  persistent. 

Recall  that  s  is  the  death  rate  of  the  predator  y  in  the 
absence  of  prey  and  k  is  the  conversion  factor  from  prey  biomass  to 
predator  biomass.  Condition  (iii)  is  a  condition  to  guarantee  y  sur¬ 
vives.  It  says  that  Pg(K)  should  be  large  enough  so  that  y  will  not 
go  extinct.  Condition  (ii)  is  a  condition  to  guarantee  neither  Eg  nor 
E^  are  asymptotically  stable.  It  says  that  Pg(x.j)  (resp.  Pg(Xg)  ) 
should  be  small  enough  so  that  for  solutions  that  start  near  Eg  (resp. 
E^),  the  gamete  type  a  (resp.  A)  will  not  become  extinct.  Note  that 
conditions  (ii)  and  (iii)  are  not  mutually  exclusive  because  x^  <  K 
( i  =1,3)  and  Pg  is  strictly  increasing. 

In  the  case  when  the  fertilities  /  viabilities  are  all  differ¬ 
ent,  the  answer  to  question  (2)  is  provided  by  Theorem  4.11.  In  order 
to  explain  the  results  of  this  theorem,  let  us  recall  the  following  assump¬ 


tions  and  notations. 


. 

t 
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It  is  assumed  that  if  the  subsystem  consisting  of  x-j  and  y 
only  contains  a  positive  interior  equilibrium  E3  -  (x^,0,0,y^)  ,  then  it 
is  globally  stable.  Similarly,  if  the  subsystem  consisting  of  x 3  and 
y  only  contains  a  positive  interior  equilibrium  =  (0,0, x^y^)  , 
then  it  is  globally  stable. 

ai  =  f.  -  f2  (i  =  1 ,3) 

bi  =  -s  +  kP.(K)  ( i  = 1,3) 

1  _ 

d.-  =  z-[-  f  -  D(x . )  -  y .P?(x. )]  (i  =  1,3) 
xi  Ti  1  w 

d  =  -s  +  — - — o[c2P-,  ( K)  +  2cP9(K)  +  P  JK)] 

(1+c)2  1  2  3 

a3 

where  c  =  —  .  Note  that  d-,  (resp.  d~)  is  defined  only  when  E9 
a  i  i  o  o 

(resp.  E^)  exists  and  d  is  defined  only  when  a-j  and  a^  are  of  the 

same  sign.  Moreover,  b.  >  0  ( i  = 1,3)  is  the  necessary  and  sufficient 

condition  for  the  existence  of  E^  ( i  =  1,3)  .  By  noting  B(x. )  ■  D(x.. ) - 

y. Pi (x. )  =  0  ,  the  meaning  of  d..  will  become  a  little  bit  more  trans- 

1  -  s  -  ai  - 

parent  if  we  rewrite  it  as  —  [y^ (-^  -  PgCx.))  -  —  B ( x ^ ) ]  . 

xi  1 

Since  there  are  thirty  two  persistence  and  non-persistence 
cases  studied,  we  will  only  highlight  the  results  of  the  theorem  by  discus¬ 
sing  two  of  these  cases  here.  (See  Theorem  4.11  for  the  labelling  of  the 
cases. ) 
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(I)  Persistence  due  to  predation. 

In  case  (8c):  a^  <  0  ,  a3  >  0  ,  >  0  ,  b3  >  0  ,  d^  >  0  , 

d3  >  0  ,  system  (4.1)  is  persistent.  Since  a-j  <  0  and  a3  >  0  ,  we 
know  that  gamete  A  (consequently  and  x2  )  will  go  extinct  in 
the  absence  of  the  predator  y  .  b-j  ,  b3  >  0  says  that  y  can  survive 
on  either  x-j  or  x3  alone,  d.  >  0  (i  =  l,3)  reduces  to 

_  s  a  •  _ 

Pp ( x • )  <  tt  - ~  B(x.)  .  Hence  the  persistence  of  the  system  is  due  to 

f.y.  1 
r  i 

Pi  (1  =  1,2,3)  . 

(II)  Persistence  due  to  fertil ity  /  viabil ity  differences. 

In  case  (lb):  a-j<0,a3<0,b^<0,b3<0,d>0, 

system  (4.1)  is  persistent.  In  the  absence  of  fertility  /  viability 
differences  (that  is,  assume  a-j  =  a3  =  0),  the  system  is  non-persistent, 
because  b.  <  0  ( i  =  1,3)  .  However,  if  the  heterozygote  x2  has  selec¬ 

tion  advantage  due  to  fertilities  /  viabilities  and  the  predator  functional 
response  P^  of  x2  is  large  enough  to  make  d  >  0  ,  then  the  system  is 
persistent. 

In  Chapter  5  the  model  studied  in  Chapter  3  is  extended  to  a 
two-locus,  two-allele  model.  Only  the  no  selection  case  is  considered. 

The  whole  population  is  again  denoted  by  x  and  the  subpopulations  of 
the  nine  genotypes:  AABB  ,  AABb  ,  AAbb  ,  AaBB  ,  AaBb  ,  Aabb  ,  aaBB  ,  aaBb  ,  and, 
aabb  are  denoted  by  x-^  ,  x^  ,  x^3  ,  x^  ,  x^^  ,  x^3  ,  x^  ,  x32  ,  and,  x33 
respectively.  The  main  result  of  this  chapter  is  Theorem  5.6  which 
provides  a  complete  picture  of  the  global  dynamics  of  system  (5.1).  It 
is  an  analogue  of  the  Hardy-Weinberg  equilibrium  relation  for  a 
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discrete  two-locus,  two-allele  with  random  mating  and  no  selection  model. 
In  summary,  this  theorem  states  the  following. 


All  solutions  of  system  (5.1)  converge  to  some  equilibrium. 
More  precisely,  if  (x^ (t) ,x12(t) ,x] 3(t) ,x21 (t) ,x22(t) ,x23(t) ,x3] (t) , 
x32(t) ,x33(t) )  is  a  solution  of  system  (5.1),  then  it  converges  to  the 

K 


equil ibri urn: 


(i+Clr(i+c2) 


2  2  2  2  2  2 
"2"  ( c-j c^ ,2c^ c2*Ci  ,2c-jC2,4c-jC2,2c-|  ,c2,2c2,l ) 


as  t  tends  to  positive  infinity.  As  in  the  one-locus  case,  c-j  ,  c2 
depend  on  initial  conditions,  namely 


[x1l(0)+x12(0)+xi3(0)l  +y  [x21  (0)+x22(0)+x23(0)l 

c  -  . 

tx31 (0)+x32 (0 )+x3s( 0) ]  +  2  [x23 (0)+x22(0)+x23(0)] 
and 

[X-j  i  (0)+x2i  ( 0 )  +x ^ -j  (0)1  +  ^  [x]2^+x22^+X32^°^ 

2  [x13(0)+x23(0)+x33(0)]  +  \  [x12(0)+x22(0)+x32(0)] 

Hence,  the  dynamics  of  system  (5.1)  is  trivial.  Rephrased  in 
terms  of  proportions,  this  theorem  also  says  that  x^x)  :  x^2(t)  : 
x]3(t)  :  x21(t)  :  x22(t)  :  x23(t)  :  x31(t)  :  x32(t)  :  x33(t)  tends  t0 

c]c\  :  2c^c2  :  c^  :  2c^c*  :  4^^  :  2c]  :  c|  :  2c2  :  1  as  t  tends  to  positive 

infinity,  which  is  the  more  usual  form  in  population  genetics  literature. 

In  Chapter  6,  the  model  in  Chapter  5  is  extended  to  include  a 
predator  y  .  The  dynamics  of  the  prey  population  is  modelled  by  system 
(5.1).  Again  the  two  questions  in  Chapter  4  are  raised.  A  partial  answer 
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to  question  (1)  is  provided  by  Theorem  6.5  which,  modulo  a  number  of 
technicalities,  can  be  summarized  as  follows. 


Assume  that  each  of  the  four  two-population  subsystems  consti¬ 
tuting  y  and  x.-  (i,j  = 1 ,3)  contains  a  globally  stable  positive 

I  vJ 

interior  equilibrium.  If  there  is  a  complete  ordering  of  the  predator 


functional  responses:  P-j -j  <  P <  P ^  <  P ^  <  P^  <  P^  <  P ^  < 
P33  >  and,  if  the  predator  y  can  survive  on  the  most  difficult  to 
catch  prey  genotype  x-j -j  ,  then  the  gamete  types  Ab  ,  aB  ,  and,  ab  go 
extinct.  Hence  only  x^  and  y  survive. 


This  result  is  intuitively  clear  in  view  of  the  similar  result 
for  the  one-locus  case. 


A  partial  answer  to  question  (2)  is  provided  by  Theorem  6.6 
which,  modulo  a  number  of  technicalities,  can  be  summarized  as  follows. 


Assume  that  the  predator  y  survives  and  that  each  of  the 
four  subsystems  consisting  of  y  and  x,-|  (resp.  xq 3  »  x3q  »  x33  ) 
contains  a  globally  stable  positive  interior  equilibrium 
E-|  =  (x,-|  ,0,*  •  •  ,0,y.|  -j )  (resp.  E2  =  (0,0,^ 3,0, •  •  •  ,0,y-j  3)  , 

E3  =  (0 ,  •  •  •  ,0  ,x"3-j  ,0 ,0  ,y3i )  ,  E33  =  (0,--*,0,x33,y33)  ).  Under  each  of 
the  following  three  cases, 


( i )  p  <p  <p  <p  <p 

11  12  13  23 


P  <P  <P  <P  <P 
33  5  11  21  31  32  33 


and  P22(xn)  <  f  -  ~  B<xn)  > 


2y 


11 


33 


<  P12  < 

P13 

’  P11  <  P21  <  P31 

<  P32  < 

P31 

*  P22^X11  ^  <  k  " 

’  P33  <  P23  <  P13  * 


2  y 


11 
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P??(x3?)  <  k  ~ —  B(Xgg)  »  anc*> 


2y 


33 


fiiil  P  <P  <P  <P  <P  <P  <P 
11  12  *13  23  33  32  31 

and  P22(xll)  <  |  '  8(7^)  , 


D  <  P  <  P 

5  11  21  31 


2y 


11 


system  (6.1)  is  persistent. 


One  could  interpret  these  results  as  saying  that  if  we 
consider  the  system  as  a  one-locus  problem,  then  the  system  is  non-persis¬ 
tent;  but  as  a  two-locus  problem,  the  predator  functional  response  P^ 
of  the  double  heterozygote  x ^  can  change  (by  making  more  diffi¬ 

cult  to  catch)  a  non-persistent  system  to  one  of  persistence. 
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APPENDIX  1 


DERIVATION  OF  THE  ONE-LOCUS  TWO-ALLELE  MODEL  (3.1) 


If  we  denote  the  two  alleles  at  the  locus  by  A  and  a  , 
then  there  are  three  genotypes:  AA  ,  Aa  ,  and,  aa  .  The  mating  table 
looks  like: 


Mating  Type 

AA 

Aa 

aa 

AAx  AA 

1 

- 

- 

AA  x  Aa 

1/2 

1/2 

- 

AA  x  aa 

- 

1 

- 

Aa  x  Aa 

1/4 

1/2 

1/4 

Aa  x  aa 

- 

1/2 

1/2 

aa  x  aa 

- 

- 

1 

There  are  two  ways  of  arriving  at  the  coefficients  of  B(x) 
in  system  (3.1). 

(i)  f .  '  s  as  fertil ities: 

f.  (i  =1,2,3)  is  to  be  thought  of  as  the  number  of  gametes 
produced  by  each  individual  in  the  x.  population.  For 
example,  there  will  be  j  f^x^  A  gametes  and  j  f^x^  a 
gametes  produced  by  the  x^  population.  Thus,  there  are 
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f-|xl  +  -Jr  ?2x2  =  u  ^at  many  A  gametes  and 

^3X3  +  \  ^2X2  =  V  that  man-y  a  9ametes  ln  the  gamete  pool. 

By  random  union  of  gametes,  the  proportion  of  AA  genotype 

2 

produced  is:  (77777)  (77777)  =  — - — 9  .  Similarly  one  can  arrive 

u  v  u  v  (u+vp 

at  the  other  coefficients  in  system  (3.1)  this  way. 

(ii)  f.  's  as  viabilities: 

f  .x. 

The  fraction  .  V1  ■■  (i  =  1,2,3)  is  to  be  thought 
Vl  T2X2  3X3 

of  as  the  proportion  of  x-  that  involves  in  random  mating. 
Since  the  AA  genotype  is  obtained  from  the  mating  types: 
AA*AA  ,  AAxAa  ,  Aa*AA  ,  and,  Aa*Aa  with  probabilities:  1  ,  1/2  , 
1/2  ,  and,  1/4  respectively,  we  have 


(u+v)2 

as  the  proportion  of  AA  genotypes  produced.  The  other  two 
coefficients  of  B(x)  in  system  (3.1)  are  arrived  at  in  a 
similar  way. 


( 


Vi 


-)(- 


Vl 


f  1  xi  +  f2x2  +  V^lx]  +  ^2X2  +  ^3X3 


) 


+  1( 


Vl 


-)(: 


f2X2 


2  vf^  x^  +  f2x2  +  f^x^  x-|  +  ^2x2  +  ^3X3 


•) 


+ 1  ( 


V2 


-)(: 


flXl 


2  VVf2X2  +  f3X3  Vl+f2X2  +  f3X3 


) 


+  1  ( 


V2 


-)(: 


f2x2 


4  1 -|x-|  +  f2x2  +  *3*3'  Vflxl  +  f2X2  +  f3X3 


)  = 


APPENDIX  2 


DERIVATION  OF  THE  TWO-LOCUS  TVJO- ALLELE  MODEL  (5.1) 


The  derivation  is  similar  to  the  one-locus  case  in  Appendix  1. 
We  denote  the  two  allele  types  at  the  first  locus  by  A  and  a  and  the 
two  allele  types  at  the  second  locus  by  B  and  b  .  There  are  nine  geno¬ 
types:  AABB  ,  AABb  ,  AAbb  ,  AaBB  ,  AaBb  ,  Aabb  ,  aaBB  ,  aaBb  ,  and,  aabb  , 
the  number  of  which  we  denote  by  ,  x^  >  »  X33  respectively. 

Therefore  the  number  of  AB  (resp.  Ab  ,  aB  ,  ab)  gametes  in  gamete  pool 
is  given  by  fu-j  (resp.  fu^  ,  fu^  ,  fu^)  where  the  u^  's  were  defined  in 
(5.2)  and  f  is  the  common  fertility  of  the  genotypes.  The  total  number 
of  gametes  in  the  gamete  pool  is  given  by  fu-j  +  fu^  +  fu^  +  fu^  =  fx  .  By 
random  union  of  gametes,  the  proportion  of  AABB  genotype  produced  is: 
fu,  fu,  u^ 

(-^-)  (-^-)  =  — ^  .  Similarly  one  can  arrive  at  the  other  coefficients  of 

X 

B(x)  in  system  (5.1)  this  way. 


The  alternative  derivation  given  in  Appendix  1  can  also  be 
carried  out.  Since  it  is  lengthy  but  similar,  we  will  not  write  it  down 
here. 
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